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A  test  apparatus  for  measuring  the  effects  of  waviness  is  described. 

A  computer  based  system  is  used  to  record  the  data  and  monitor  the  test. 

A  test  procedure  is  outlined. 

Results  from  wavy  seal  experiments  show  that  waviness  reduces  seal 
wear  compared  to  a  flat  face  seal  a  factor  of  four  to  forty  and  reduces 
torque  a  factor  of  four.  Long  term  tests  show  that  by  moving  the  wave  the 
effect  remains  indefinitely.  As  the  amount  of  waviness  is  increased,  fric¬ 
tion  torque  decreases  and  leakage  increases.  Wear  profiles  show  that  the 
seal  faces  have  a  tilt  in  addition  to  waviness.  This  alters  the  pressure 
distribution  significantly  when  compared  to  the  radially  parallel  wavy 
faces  analyzed  previously.  Tests  using  a  converging  radially  tapered  seal 
show  that  the  taper  wears  away  with  time.,. 

A  theory  is  developed  for  predicting  tnermoelastic  instability  in  seal 
rings.  Results  show  that  ring  deflection  has  a  predominant  influence  at 
low  numbers  of  waves.  Experimental  results  are  needed  to  verify  the  theory. 

The  hydrodynamic  theory  previously  developed  for  wavy  face  seals  is 
completely  revised  and  presented  as  a  more  comprehensi ve  theory.  This  was 
necessary  because  of  the  tilt  variation  on  the  seal  faces.  Using  this  theory 
it  is  shown  that  certain  combinations  of  waviness  and  tilt  lead  to  minimum 
torque  and  minimum  leakage.  This  optimum  design  has  faces  which  are  either 
radially  converging  or  parallel  with  contact  all  around  the  seal  to  provide 
low  leakage.  These  results  point  the  way  to  even  better  seal  designs  than 
the  experimental  seal. 

Using  the  model  a  comparison  between  theory  and  experiment  was  made. 

The  theoretical  solution  requires  that  a  complete  deflection  analysis  and 
simultaneous  solution  be  made.  Agreement  between  experiment  and  theory  is 
good  for  both  torque  and  leakage.  Agreement  of  theory  with  measured  deflec¬ 
tion  is  excellent. 

In  conclusion,  the  report  shows  both  theoretically  and  experimentally 
how  waviness  can  be  used  to  great  advantage  in  a  water  seal  to  reduce  fric¬ 
tion  and  wear. 
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APPLICATION  OF  RESEARCH  TO  THE  NEEDS 
OF  THE  U.S.  NAVY 


Mechanical  face  seals  are  used  in  numerous  applications  in  Naval 
machinery.  These  applications  range  from  propeller  shaft  seals  to 
boiler  feed  pump  seals.  In  such  equipment  the  mechanical  seal  plays 
a  vital  role.  When  such  seals  fail,  repair  is  costly  both  in  terms 
of  lost  time  and  direct  costs,  so  any  improvement  in  seal  life  and 
reliability  would  be  of  significant  benefit. 

As  more  advanced  equipment  is  designed,  it  is  sometimes  difficult 
to  achieve  desired  performance  in  more  severe  service  environments 
with  the  present  state  of  the  art  of  seal  design.  Thus,  an  improve¬ 
ment  in  seal  technology  would  serve  this  important  application. 

The  immediate  objective  of  the  research  herein  is  to  further  the 
understanding  of  mechanical  face  seal  lubrication  phenomena.  The 
ultimate  objective  is  to  develop  the  capability  of  designing  contact¬ 
ing  face  seals  having  a  longer  life,  greater  reliability,  and  for 
extreme  environments.  Thus,  the  objectives  of  this  research  are  com¬ 
patible  with  mechanical  face  seal  needs  for  Naval  machinery. 
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Chapter  1 
Introduction 


Mechanical  Face  Seals 

In  applications  where  a  rotating  shaft  must  pass  from  one  fluid 
region  to  another,  contacting  mechanical  face  seals*  play  the  essen¬ 
tial  role  of  minimizing  the  transfer  of  fluid  between  the  regions. 
Applications  of  face  seals  range  from  water  pump  seals  to  process 
pump  seals  to  propeller  shaft  seals. 

The  performance  and  reliability  of  contacting  mechanical  face 
seals  are  of  great  importance  for  any  type  of  equipment  where  minimal 
leakage,  high  reliability,  and  long  life  are  necessary.  Even  for 
equipment  where  these  factors  are  not  so  critical,  seal  failures  and 
short  seal  life  lead  to  high  operating  cost  due  to  down  time  and  main¬ 
tenance  cost. 

Even  though  mechanical  face  seal  technology  has  been  steadily 
improving  over  the  past  several  decades,  further  improvement  in  the 
state  of  the  art  of  seal  design  would  be  most  beneficial.  Although 
seals  having  an  acceptable  life  and  reliability  can  be  designed  for 
many  applications,  further  improvement  in  seal  life  and  reliability 
would  result  in  significant  cost  savings  to  the  user.  Also,  there 
are  numerous  mechanical  face  seal  applications  where  seal  loading, 
reliability,  life,  and  leakage  requirements  are  difficult  to  achieve 
within  the  present  state  of  the  art.  Examples  of  such  applications 
are  seals  for  pumps  for  nuclear  power  plants  and  seals  for  large 
diameter  submarine  propeller  shafts.  Additionally,  the  friction 
losses  in  face  seals  represent  a  significant  fraction  of  energy  con¬ 
sumed  for  pumping  purposes.  Within  the  present  state  of  the  art  of 
seal  design, it  is  very  difficult  to  design  a  low  leakage  seal  that 


*The  class  of  low  leakage  face  seals  where  there  is  definite  contact 
and  wear  of  the  faces  as  opposed  to  hydrostatic  or  hydrodynamic 
where  a  definite  clearance  is  maintained. 


1 


also  has  a  low  friction  loss. 

The  main  barrier  to  the  advancement  of  the  state  of  the  art  is 
that  the  mechanics  of  seal  operation  are  not  well  enough  understood 
to  be  able  to  reasonably  anticipate  seal  performance  as  a  function  of 
design  parameters.  There  are  no  well  established  fundamental  theore¬ 
tical  bases  that  can  be  used  to  indicate  the  type  of  seal  design  that 
will  give  improved  performance.  Improvements  that  have  been  made  have 
been  brought  about  largely  by  trial  and  error  combined  with  elementary 
sealing  theories. 

In  order  to  be  able  to  predict  the  performance  of  contacting  face 
seals  as  a  function  of  design  parameters,  it  is  essential  that  the  lub¬ 
rication  mechanisms  between  the  faces  be  well  understood.  At  present 
it  is  known  that  hydrodynamic  or  hydrostatic  lubrication  plays  some 
role  in  providing  load  support  for  oil  seals  as  well  as  water  seals. 

But,  the  precise  nature  of  this  lubrication  is  not  known.  Several 
theories  have  been  put  forth.  However,  these  theories  have  not  been 
verified  for  contacting  face  seals,  and  it  is  not  possible  as  yet  to 
use  these  theories  for  the  design  of  contacting  face  seals. 

In  this  work,  the  results  from  further  research  into  hydrodynamic 
and  hydrostatic  lubrication  mechanisms  in  face  seals  is  reported.  Much 
has  been  learned  about  these  mechanisms.  Based  upon  experimental  work 
there  is  a  strong  promise  that  effective  use  of  hydrodynamic  and  hydro¬ 
static  mechanisms  can  be  made  to  dramatically  improve  face  seal  per¬ 
formance. 

Seal  Lubrication 

The  mechanical  face  seal  consists  basically  of  two  annular  rings 
which  rotate  relative  to  each  other  and  which  are  pressed  together  by 
spring  and  fluid  pressures  (see  Figure  1-1).  The  surfaces  that  rub 
together  are  generally  manufactured  as  flat  as  possible  initially  so 
as  to  minimize  leakage.  The  effective  gap  between  the  faces  is  ideally 
quite  small  (order  of  1  ym)  so  that  leakage  flow  across  the  faces  will 
be  quite  small.  The  difficulty  in  designing  a  mechanical  seal  is  in 
maintaining  the  gap  at  a  very  low  value  while  at  the  same  time  providing 
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Figure  1-1.  Mechanical  Face  Seal. 
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a  definite  lubricant  film  between  the  faces. 

The  load  that  must  be  supported  at  the  faces  of  a  mechanical  seal 
is  due  primarily  to  loading  caused  by  the  sealed  pressure.  The  unit 
face  load  can  be  expressed  as  some  fraction  B  of  the  sealed  pressure 
where  B  can  be  made  greater  or  less  than  unity  by  geometry  selection. 

The  load  support  at  the  faces  is  derived  from  fluid  pressure  and  me¬ 
chanical  pressure.  If  the  fluid  pressure  at  the  faces  is  large  enough 
to  support  all  of  the  load,  then  there  will  be  no  contact  and  no  ad¬ 
hesive  wear.*  If  none  of  the  load  is  supported  by  fluid  pressure,  the 
load  must  be  carried  by  mechanical  contact  and  the  wear  rate  will  be 
large. 

In  practice,  seals  often  operate  at  one  of  two  extremes.  At  one 
extreme,  a  large  gap  will  be  created  by  hydrostatic  or  hydrodynamic 
pressure,  all  of  the  load  will  be  supported  by  fluid  pressure,  and  the 
seal  will  leak  a  large  amount.  At  the  opposite  extreme,  the  gap  will 
close  completely  and  leakage  will  be  very  low.  Only  a  small  fraction 
of  the  load  will  be  carried  by  fluid  pressure,  and  wear  and  heat  gen¬ 
eration  will  increase. 

Based  on  the  above,  it  can  be  concluded  that  an  effective  seal 
should  operate  between  these  two  extremes— having  both  adequate  fluid 
pressure  load  support  and  low  leakage.  The  seal  should  operate  so  that 
it  just  touches  to  minimize  leakage  but  such  that  the  load  is  carried 
by  fluid  pressure,  not  mechanical  pressure.  To  do  this  requires  that 
any  fluid  pressure  generation  mechanism  used  to  provide  load  support 
to  the  seal  must  be  very  carefully  controlled.  At  present,  in  commer¬ 
cial  contacting  face  seals  this  lubrication  is  left  primarily  to  chance, 
and  often  such  seals  operate  at  one  of  the  extremes  mentioned.  Quite 
commonly  such  seals  will  operate  in  the  low  leakage  condition  where 
wear  progresses  at  a  definite  rate.  Such  seals  are  quite  satisfactory 
for  many  applications.  However,  when  pressure  or  speed  requirements 
are  increased,  then  this  mode  of  operation  may  lead  to  rapid  failure 
through  wear  or  heat  checking. 

*There  may  still  be  abrasive  wear  even  if  the  surfaces  do  not  touch. 
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In  order  to  be  able  to  design  into  the  seal  the  proper  balance 
between  fluid  pressure  load  support  and  leakage,  the  origin  of  lubri¬ 
cation  mechanisms  must  first  be  fully  understood.  A  considerable  amount 
of  research  has  been  accomplished  in  this  area  over  the  years.  Various 
theories  have  emerged  and  have  been  described  in  a  previous  annual  re¬ 
port  [1].*  These  theories  fall  roughly  into  two  categories,  hydrostatic 
and  hydrodynamic.  In  this  research  program,  attention  has  been  directed 
toward  studying  the  effects  of  waviness  as  a  source  of  hydrodynamic  and 
hydrostatic  pressure  in  face  seals.  Particular  attention  has  been  fo¬ 
cused  on  hydrodynamic  effects  in  relation  to  roughness  and  low  viscosity 
applications  where  some  touching  is  expected.  Waviness  was  selected 
among  the  various  hydrodynamic  theories  because  it  is  controllable  in 
seals  and  a  better  understanding  could  lead  to  improved  seal  designs. 

Sources  of  waviness  in  seals  have  been  described  in  detail  pre¬ 
viously  [1].  Waviness  may  be  produced  accidentally  as  an  uncontrolled 
variable  or  intentionally  in  a  precisely  controlled  manner.  A  wavy 
seal  ring  is  shown  in  Figures  1-2  and  1-3. 

Wavy  Face  Seal 

The  concept  of  waviness  is  that  where  the  film  thickness  varies 
in  some  fashion  circumferentially  around  the  seal.  Generally  speaking, 
film  thickness  may  vary  radially  as  well  as  tangentially. 

h  =  h(r,e)  .  (1-1) 

In  the  present  work  interest  is  focused  upon  film  thickness  shapes  of 
the  following  functional  form 

h  =  hg  +  f^(r)  cos  ne  +  f£(r)  sin  ne  .  (1-2) 

At  any  particular  radius  r  the  film  shape  is  periodic  with  n  waves 
around  the  seal  and  is  therefore  wavy.  However,  film  shape  can  also 
vary  in  some  general  manner  with  r. 


♦Numbers  in  brackets  refer  to  List  of  References  at  the  end  of  report. 
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Figure  1-2.  Wavy  Seal  Geometry. 


If  f^(r)  =  const  and  f 2 ( r )  »  const,  then  the  faces  are  always  radi 
ally  parallel.  This  component  of  film  thickness  variation  is  commonly 
termed  waviness.  If  the  f^(r)  f  const,  then  the  faces  are  not  in  gen¬ 
eral  radially  parallel.  The  f^(r)  are  referred  to  as  tilt.  Thus,  the 
film  thickness  shapes  of  interest  are  combinations  of  waviness  and  tilt 
Since  at  any  radius  the  film  thickness  is  wavy,  the  combination  of 
waviness  and  tilt  defined  above  will  also  be  called  a  wavy  film  shape. 

The  reason  for  choosing  film  shapes  as  described  by  equation  (1-2) 
as  a  subject  for  study  is  that  these  shapes  can  conveniently  be  gener¬ 
ated  by  planned  mechanical  distortions  in  a  seal  ring,  and  the  shapes 
also  include  common  modes  of  unplanned  distortion  found  in  operating 
seals.  For  example,  generally  seals  undergo  a  uniform  tilt  due  to 
pressure  and  thermal  deformation.  When  rings  are  loaded  by  any  no n- 
axisymmetric  load  they  become  wavy  as  described  by  equation  (1-2). 

Now,  for  the  sake  of  illustration,  assume  a  seal  has  a  wavy  film 
thickness  shape  given  by  equation  (1-2)  where  n  =  2.  Assume  a  ring  is 
pressed  against  a  perfectly  flat,  moving,  mating  ring.  For  n  =  2  there 
will  generally  be  2  or  4  regions  where  film  thickness  is  smallest  (only 
2  are  shown)  and  where  contact  can  be  expected  to  occur  (see  Figures 
1-2  and  1-3).  In  general,  a  fluid  pressure  distribution  develops 
according  to  the  Reynolds  equation.  Both  hydrostatic  and  hydrodynamic 
components  are  included.  Angular  rotation  of  the  mating  face  causes 
pressure  to  build  up  in  the  converging  regions  as  the  minimum  film 
thickness  region  is  approached.  In  diverging  regions,  cavitation  will 
occur  as  the  pressure  attempts  to  become  negative.  For  an  outside 
pressurized  seal  the  cavity  may  appear  as  shown  in  Figure  1-2  depending 
upon  the  f^(r).  The  fluid  will  flow  in  streamers  across  the  cavity  and 
a  full  film  will  again  develop  at  the  downstream  cavity  boundary. 

The  problem  to  be  solved  is  to  determine  the  fluid  and  mechanical 
pressure  distribution  for  the  given  configuration.  Given  the  pressure 
distribution,  then  load  support  and  leakage  can  be  calculated,  and  the 
effect  of  various  parameters  can  be  studied.  Given  this  tool,  then 
waviness  shapes  which  provide  improved  performance  can  be  found. 
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During  the  first  portion  of  this  research  program  only  radially 
parallel  face  waviness  was  considered.  Much  effort  was  spent  in  devel¬ 
oping  solution  methods  and  studying  predicted  results. 

As  background,  the  pressure  distribution  for  wavy  parallel  faces 
has  been  solved  by  various  methods  by  Findlay  [2],  Pape  [3],  and 
Stanghan-Batch  and  Iny  [4]  for  perfectly  smooth  faces.  These  results 
show  that  even  a  small  waviness  produces  sufficient  load  support  for 
complete  liftoff.  However,  heavily  loaded  or  low  viscosity  seals  do 
not  operate  with  full  fluid  pressure  load  support  and  complete  separa¬ 
tion.  A  definite  wear  results.  Surfaces  of  such  seals  contact  during 
operation,  and  pressure  distribution  is  affected  by  interactions  with 
surface  roughness  (see  Figure  1-3).  Seals  of  this  type  operate  in  a 
mixed  friction  regime.  Fluid  pressure  load  support  due  to  waviness 
may  provide  a  significant  fraction  of  the  load  support,  but  asperity 
contact  must  provide  the  balance. 

As  a  starting  point  for  this  research  program,  these  important 
effects  were  included  in  a  wavy  seal  lubrication  model.  In  the  first 
annual  report  for  this  project,  reference  [5],  this  more  general  prob¬ 
lem  was  solved  using  a  one-dimensional  theory.  In  the  second  annual 
report,  [1],  the  much  more  complex  two-dimensional  solution  to  the 
above  problem  was  solved.  The  effects  of  waviness,  roughness,  asper¬ 
ity  contact,  wear,  cavitation,  and  elastic  deflection  were  included 
in  this  model.  Using  this  model,  predictions  were  made  for  the  rela¬ 
tive  wear  rate,  friction,  and  leakage  as  a  function  of  roughness, 
waviness,  speed,  size,  pressure,  viscosity,  and  material. 

A  number  of  conclusions  were  reached  based  on  these  first  two 
annual  reports: 

1)  The  effects  of  roughness  on  hydrodynamic  lubrication  are  not 
completely  understood.  Certain  fundamental  questions  remain 
concerning  the  roughness  model  used. 

2)  As  to  the  potential  of  utilizing  hydrodynamic  effects  caused 
by  parallel  face  waviness  to  advantage  by  design,  the  results 
show  that  wear  rate  and  friction  can  be  greatly  reduced  while 
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maintaining  leakage  at  acceptable  levels. 

3)  While  a  comparison  of  predicted  results  to  experimental  re¬ 
sults  given  in  the  literature  is  generally  good,  data  con¬ 
tained  in  the  literature  are  incomplete,  so  more  complete 
experimental  data  are  needed  for  comparison. 

4)  In  low  viscosity  or  heavily  loaded  applications  where  some 
touching  is  expected  to  occur,  waviness  will  wear  away  with 
time  and  any  benefit  derived  will  be  lost. 

5)  Based  upon  data  for  some  commercial  seals  and  using  the  model, 
it  was  determined  that  there  was  insufficient  accidentally 
caused  waviness  to  produce  significant  hydrodynamic  effects 

in  water.  One  cannot  generalize  to  say  that  such  effects  do 
not  occur  in  commercial  seals.  However,  the  model  is  provided 
so  that  the  question  can  be  answered  on  a  case  by  case  basis. 

Item  1)  was  treated  extensively  in  the  third  annual  report  [6]. 
Even  after  this  analysis  certain  fundamental  questions  remain  concern¬ 
ing  how  to  deal  with  roughness  in  lubrication  problems.  A  contribution 
to  the  literature  on  the  subject  was  made  [7].  This  thorough  analysis 
led  to  conclusions  allowing  certain  simplifying  assumptions  discussed 
in  Chapter  3  of  this  report. 

Item  2)  was  also  treated  extensively  in  the  third  annual  report 
[6],  A  methodology  for  the  design  of  a  wavy  face  seal  was  developed 
and  applied.  Theoretical  results  showed  large  reduction  in  friction 
and  wear  rate  compared  to  conventional  designs  whereas  leakage  could 
be  controlled. 

Concerning  item  3),  the  second  and  third  annual  reports  [1,6] 
describe  a  test  apparatus  designed  to  test  the  wavy  seal  theory. 

This  apparatus  has  been  in  operation  for  more  than  one  year,  and 
many  tests  have  been  conducted.  These  test  results  are  reported  in 
Chapter  2  and  compared  to  theory  in  Chapter  4  of  this  report. 

As  will  be  discussed,  the  test  results  do  verify  predicted  behav¬ 
ior.  However,  early  in  the  test  program  it  was  observed  that  film 
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thickness  was  not  of  the  radially  parallel  face  type,  but  consisted  of 
combined  waviness  and  tilt  as  previously  described.  This  necessitated 
modification  of  the  waviness  model  to  include  this  more  generalized 
type  of  film  thickness.  These  results  are  reported  in  Chapter  4  of 
this  report. 

Concerning  item  4)  above,  a  solution  to  this  problem  was  first 
proposed  in  the  first  annual  report  [5],  It  was  proposed  to  move  the 
waviness  slowly  around  the  seal  so  that  whatever  wear  occurred  would 
be  uniformly  distributed.  Then  the  shape  of  the  wave  would  be  pre¬ 
served  and  tests  using  a  constant  wave  could  be  made.  The  concept  is 
illustrated  in  Figure  1-4.  This  concept  was  incorporated  into  the 
test  apparatus  and  is  described  in  detail  in  references  [1]  and  [6]. 

In  summary,  in  this  report  experimental  results  using  waviness 
are  presented,  the  waviness  theory  is  modified  to  accommodate  certain 
unanticipated  film  thickness  shape  characteristics,  and  theory  and 
experiment  are  compared.  A  related  development,  that  of  the  role  of 
thermoelastic  instability  in  face  seals,  is  also  included  in  Chapter  3. 
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Distorted  Ring  Under  Pressure 
After  Wave  Has  Moved 


Figure  1-4.  Moving  Waviness  Concept 


Chapter  2 

Experimental  Results 


Test  Apparatus 

In  order  to  test  the  various  theories  of  seal  lubrication  devel¬ 
oped  in  this  research  program,  a  mechanical  seal  test  apparatus  was 
designed  and  fabricated.  The  details  of  the  design  are  included  in 
the  previous  summary  report  [6].  The  important  features  of  the 
apparatus  plus  changes  and  additions  since  the  previous  report  will 
be  discussed  here. 

The  design  of  the  test  apparatus  satisfies  several  criteria.  A 
wide  range  of  speed  and  sealed  pressure  is  available  in  order  to  ver¬ 
ify  modeled  behavior  over  a  wide  range  of  conditions.  Certain  other¬ 
wise  uncontrolled  variables  such  as  shaft  runout,  seal  waviness,  and 
seal  radial  taper  are  very  carefully  controlled.  For  testing  the  lu¬ 
brication  model  of  interest,  a  controllable  waviness  can  be  applied 
to  the  seal  faces,  and  in  order  to  reduce  the  effects  of  wear,  this 
waviness  is  moved  circumferentially  with  respect  to  the  primary  ring. 
Leakage  from  the  seal  is  measured.  Torque  is  measured  in  a  manner 
such  that  the  seal  torque  is  measured  independently  of  bearing  torque 
or  torque  arising  from  a  second  seal  used  in  a  tandem  arrangement. 

Many  of  the  distortions  commonly  found  in  cormercial  seals  are  mini¬ 
mi  zed . 

The  specifications  are  as  follows: 

Configuration:  Horizontal  shaft,  single  seal  with  torque  measured  on 
the  seal  housing.  Rotating  secondary  seal  ring. 

Size:  Accepts  up  to  92  mm  (3.63  in.)  nominal  shaft  size  seal. 

102  mm  (4  in.)  nominal  seal  face  diameter. 

Speed:  500  to  4000  RPM.  3  to  21  m/s  (9  to  70  ft/s)  for  a 

102  mm  (4  in. )  seal . 

Pressure:  0  to  6.9  MPa  (0  to  1000  psi)  continuous.  Vessel  is 

rated  at  13.8  MPa  (2000  psi). 
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Bearings: 


Five  ABEC-7  bearings  on  precision  spindle.  Bearings 
are  rated  at  a  maximum  of  53000  N  (12000  lb)  thrust 
at  4000  RPM  (corresponds  to  a  102  mm  (4  in.)  seal  at 
6.9  MPa  (1000  psi). 

3700  W  (5  Hp)  variable  speed  drive  (belt  type). 

Materials  are  selected  for  seawater  service. 

Balance  ratio  of  1.0  (lower  values  can  be  used). 

Zero  moment  design.  Designed  to  minimize  the  effects 
of  pressure  caused  radial  taper. 

Seal  waviness  is  controlled  by  variable  ring  distortion. 
It  can  be  moved  relative  to  the  primary  ring. 

Can  be  controlled  at  any  given  temperature  above  38°C 
(100°F) . 

Speed,  pressure,  temperature,  waviness,  torque,  and 
leakage. 

Figure  2-1  shows  the  general  layout  of  the  test  apparatus.  The 
bolted  vessel  on  the  right-hand  end  contains  the  pressurized  seal. 
Immediately  to  the  left  is  the  waviness  drive  plate  and  the  cage-like 
torque  sensing  element.  The  entire  test  vessel  is  supported  cantilever 
style  by  the  torque  sensing  element.  Further  to  the  left  is  the  bear¬ 
ing  housing  which  is  attached  to  the  frame.  The  variable  speed  motor 
drive  is  located  to  the  left  of  the  bearing  housing. 

On  the  level  below  the  test  vessel  in  Figure  2-1  is  located  all 
of  the  auxiliary  equipment.  To  the  left  end  is  the  oil  supply  system 
including  a  tank,  pump,  filter,  and  cooler.  On  the  right  end  is  the 
high  pressure  water  circulation  pump  and  heat  exchanger. 

Figure  2-2  shows  the  essential  features  of  the  seal  test  appa¬ 
ratus  in  detail.  The  shaft  turns  on  five  precision  angular  contact 
bearings  designed  to  operate  under  the  large  thrust  load.  Cooling 
water  is  circulated  through  the  test  chamber  through  openings  not 
shown.  The  torque  sensing  element  acts  as  the  support  for  the  vessel. 
The  pressure  vessel  is  designed  to  be  readily  taken  apart  by  removal 
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of  bolts.  The  secondary  ring  support  rotor  is  removed  by  removing  the 
capscrew.  To  minimize  axial  runout  at  the  seal  face,  the  end  of  the 
shaft  is  ground  true  with  the  bearings.  The  rotor  is  ground  both  where 
it  mates  with  the  shaft  and  the  secondary  seal  ring.  The  two  secondary 
seal  ring  faces  are  ground  parallel. 

The  primary  ring  is  made  of  Pure  Carbon  P658RC.  The  unusual  pro¬ 
portions  of  the  ring  result  from  adjusting  the  geometry  so  that  sealed 
pressure  causes  no  rotation  about  a  circumferential  axis.  The  primary 
ring  is  driven  by  lugs  which  engage  two  notches.  Two  as  opposed  to  a 

greater  number  of  notches  were  used  because  this  arrangement  gives  the 

primary  ring  maximum  freedom  to  float  radially  to  align  itself.  Second 
harmonic  waviness  is  minimized  by  driving  the  ring  through  the  centroid 
of  its  cross  section.  Any  remaining  second  harmonic  waviness  produced 
by  drive  forces  will  be  more  easily  flattened  than  the  third  harmonic 
waviness  built  into  the  ring,  so  it  is  not  considered  important.  The 
secondary  seal  is  located  at  the  left  end  on  the  inside  diameter.  The 

balance  ratio  for  the  design  shown  is  1.0.  Springs  to  the  left  of  the 

seal  provide  a  preload  of  about  0.2  MPa  (30  psi). 

The  secondary  seal  ring  is  fabricated  from  tungsten  carbide.  This 
ring  is  also  of  a  zero  pressure  moment  design.  An  unusual  feature  of 
the  particular  design  is  that  the  mechanical  force  at  the  right-hand 
side  of  the  seal  has  been  reduced  to  nearly  zero  by  placing  the  0-ring 
seal  on  the  right-hand  face  as  opposed  to  the  outside  diameter.  It  is 
useful  to  minimize  the  axial  force  to  reduce  coning  effects  caused  by 
friction  at  the  load  bearing  point  in  conjunction  with  radial  deforma¬ 
tion  due  to  pressure  and  temperature.  The  secondary  seal  ring  is  held 
in  place  by  spring  clips  (not  shown)  which  engage  notches  in  the  ring. 

Figures  2-2  and  2-3  show  how  the  moving  wave  concept  was  imple¬ 
mented  in  the  seal  test  apparatus.  The  drive  ring  is  designed  so  that 
the  primary  ring  can  float  and  align  itself  but  cannot  rotate  about  the 
shaft  axis.  The  primary  ring  thus  takes  its  alignment  from  the  face  of 
the  rotating  secondary  ring.  The  waviness  drive  cylinder  is  driven  by 
a  worm  at  low  speed  (one  revolution  per  9.47  h).  The  waviness  drive 
cylinder  seals  into  the  housing  with  two  0-rings  designed  so  that  gas 
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pressure  can  be  supplied  through  the  end  plate  and  into  the  waviness 
cylinder.  Each  of  the  six  holes  on  the  waviness  drive  cylinder  is 
connected  to  this  gas  pressure.  When  the  primary  ring  is  positioned 
over  these  six  holes  and  sealed  off  by  six  0-rings,  then  the  primary 
ring  becomes  subjected  to  six  small  regions  of  gas  pressure  on  its 
inside.  These  forces  produce  moments  about  the  centroid  of  the  ring. 
The  staggered  location  causes  alternating  moments  of  opposite  sign. 

Face  waviness  and  an  alternating  face  taper  are  produced  by  these 
moments.  As  the  waviness  cylinder  is  turned  by  the  worm  and  gear, 
the  waviness  pattern  moves  relative  to  the  seal  primary  ring. 

The  torque  sensing  element  is  shown  in  Figures  2-1  and  2-2.  The 
axial  load  that  must  be  carried  by  this  member  is  large,  53000  N  at 
6.9  MPa  (12000  lbs  at  1000  psi)  whereas  the  torque  to  be  measured  is 
quite  small  (down  to  a  few  N*m  (in. -lb)).  The  axial  load  is  carried 
by  ten  3  x  25  x  102  mm  (0.120  x  1  x  4  in.)  long  elements.  These  ele¬ 
ments  have  a  very  high  axial  stiffness  and  strength  but  low  torsional 
stiffness.  The  torsional  load  is  sensed  by  two  beams  which  have  a  low 
axial  stiffness.  The  axial  and  torsional  loads  are  in  effect  isolated 
so  that  torque  can  be  measured  independent  of  axial  load.  Strain  gages 
are  mounted  on  the  two  beams.  Test  results  verify  the  design  behavior. 
There  is  no  measurable  influence  of  pressure  on  torque  at  least  up  to 
300  psi  where  the  test  was  performed. 

The  torque  sensing  element  is  sensitive  to  temperature.  A  temper¬ 
ature  increase  of  the  bearing  housing  causes  a  torque  shift  of  0.109 
(Thousing  -  Treference)  in.lb/°F.  This  effect  is  repeatable  and  has 
been  corrected  in  the  data  reduction. 

Waviness  Drive 

Since  the  previous  report  was  written  several  additions  have  been 
made  to  the  system.  The  first  addition  is  the  waviness  drive.  Figure 
2-4  shows  a  motor  driven  gear  arrangement  which  connects  to  the  shaft 
of  the  worm  which  drives  the  waviness  cylinder.  The  apparatus  shown 
normally  drives  the  waviness  cylinder  at  one  rotation  per  9.47  h  for 
the  purpose  of  providing  a  slowly  moving  wave  for  testing  the  effects 
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Figure  2-4.  Waviness  Drive. 
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of  waviness.  By  sliding  gear  A,  shown  in  the  figure,  the  motor  can  be 
connected  directly  to  the  worm  so  that  the  waviness  cylinder  speed  is 
about  1/2  RPM.  This  higher  speed  arrangement  is  used  for  setup  and 
calibration  purposes. 


Leakage  Measuring  Device 

A  second  addition  to  the  test  apparatus  is  a  leakage  measuring 
device.  For  early  tests  (up  through  test  number  37),  leakage  was  mea¬ 
sured  using  a  graduated  cylinder  and  clock.  This  had  the  disadvantage 
that  a  continuous  measurement  of  leakage  rate  was  not  available.  Also 
certain  inconsistencies  occurred  because  the  leakage  could  not  be  aver 
aged  over  a  sufficiently  long  period  of  time  in  order  to  average  out 
leakage  fluctuations  (described  later). 

To  electronically  record  leakage  on  a  continuous  basis  presents 
some  problem  when  leakage  rate  may  vary  over  an  extremely  wide  range, 
say  from  a  fraction  of  a  cm3/min  to  100  cm3/min.  Since  a  computer  is 
part  of  the  system,  the  approach  used  is  to  measure  the  time  required 
for  a  given  small  quantity  of  flow  and  compute  the  average  flow  rate. 
This  has  been  accomplished  using  the  tipping  device  shown  in  Figure 
2-5.  The  device  tips  alternately  after  about  2.0  cm3  of  flow.  Each 
time  tipping  occurs,  a  photodiode  with  associated  circuitry  closes  a 
relay  which  in  turn  closes  a  second  relay.  The  computer  is  then  in¬ 
formed  of  this  event  and  the  time  is  noted.  By  noting  the  time  for 
each  tip,  the  leakage  rate  is  calculated  after  each  tip  and  recorded 
and  plotted. 

The  device  shown  in  Figure  2-5  is  fabricated  of  thin  sheet  brass 
and  the  tipping  bucket  is  supported  by  jeweled  bearings  to  minimize 
friction  and  wear.  Calibration  tests  show  the  device  is  consistent 
and  reliable.  At  high  leakage  rates  (100  ml/m)  some  error  is  intro¬ 
duced  by  the  time  response  of  the  bucket.  However,  most  of  the  leak¬ 
age  rates  measured  are  well  below  this  rate  so  that  no  significant 
error  is  introduced. 


Leakage  Inlet 


Figure  2-5.  Leakage  Measuring  Device. 


Computer  Control  and  Pat a  Acquisition 

During  August  1979  all  of  the  data  gathering,  storing,  display, 
and  control  functions  were  taken  over  by  computer.  The  components 
of  the  original  control  system  were  retained  as  a  back  up  and  so  that 
the  apparatus  can  still  be  operated  independently  of  the  computer  as 
is  needed. 

The  computer  based  system  will  now  be  described  in  some  detail. 
The  basic  system  used  is  illustrated  in  Figure  2-6.  The  computer  can 
read  voltage,  etc.  through  the  plug  in  cards  in  the  6940  or  can  output 
signals  through  the  6940  by  other  cards  such  as  the  relay  card  shown. 
In  addition,  the  computer  controls  the  plotter  and  its  own  internal 
printer  and  tape  drive.  All  logic  is  stored  in  the  computer.  The 
control -data  acquisition  system  uses  the  program  listing  which  is 
shown  in  Appendix  A.  The  functions  of  the  program  are  as  follows: 

Pressure  -  An  analog  pressure  signal  is  supplied  by  a  BLH  pres¬ 
sure  transducer  to  the  low  level  A/D  converter  card.  This  signal  is 
read  approximately  one  time  per  second.  If  pressure  is  below  the 
pressure  set  point,  a  relay  in  the  relay  output  card  is  activated  for 
50  ms  which  in  turn  activates  a  second  relay  having  an  approximate  one 
second  latch  time.  This  second  relay  controls  the  air  valve  to  the 
high  pressure  pump.  Air  is  turned  on  for  approximately  one  second 
allowing  the  pressure  pump  to  make  one  stroke.  Using  this  system, 
pressure  can  be  controlled  to  within  about  5  psi  of  a  set  point.  If 
pressure  becomes  greater  than  120  percent  of  the  set  point  or  less 
than  70  percent  of  the  set  point,  this  indicates  a  malfunction  and 
the  apparatus  is  shut  down. 

Torque  -  An  analog  signal  is  supplied  by  the  torque  transducer 
described  in  the  previous  annual  report  [6].  This  signal  is  read 
by  the  low  level  A/D  card  approximately  12  times  per  minute.  The 
signal  is  corrected  for  a  torque  shift  caused  by  bearing  housing 
temperature,  and  the  proper  calibration  constant  is  applied  to  give 
direct  torque  readings.  The  12  points  are  averaged  and  the  average 
is  plotted  and  recorded  on  tape  along  with  the  time. 
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Calibration  of  the  torque  transducer  to  obtain  a  calibration 
constant  is  accomplished  by  a  subroutine  contained  in  the  control 
program.  Torque  transducer  voltage  output  is  read  for  zero  applied 
torque  and  a  given  applied  torque,  and  a  calibration  constant  is 
established. 

Zero  level  of  the  torque  signal  is  also  established  using  the 
control  program  at  the  beginning  and  end  of  each  test.  The  zero 
level  of  torque  is  determined  by  hand  cranking  the  drive  both  clock¬ 
wise  and  counterclockwise,  taking  a  number  of  readings  in  each  direc¬ 
tion,  and  establishing  the  voltage  level  output  that  corresponds  to 
zero  torque  level. 

Face  Temperature  -  The  low  level  A/D  card  has  a  built-in  thermo¬ 
couple  reference.  The  face  temperature  thermocouple  is  read  12  times 
per  minute,  converted  to  temperature,  averaged,  plotted,  and  stored 
on  tape.  If  face  temperature  goes  above  a  set  point,  the  apparatus  is 
shut  down. 

Leakage  -  Each  time  the  leakage  measuring  device  tips,  the  time 
of  this  event  is  recorded.  At  approximately  one  minute  intervals,  the 
number  of  tips  and  the  time  between  the  first  and  last  tip  is  deter¬ 
mined.  The  leakage  rate  is  calculated,  stored  on  tape,  and  plotted. 

If  there  have  been  no  tips  within  the  one  minute  interval,  then  this 
fact  is  noted  on  tape  and  no  point  is  plotted.  When  an  interval  is 
encountered  during  which  a  tip  has  occurred,  leak  rate  is  determined 
by  the  total  elapsed  time  between  tips  even  though  these  may  be  more 
than  one  minute  apart,  and  the  leakage  rate  is  computed,  plotted,  and 
recorded  on  tape. 

Bearing  Housing  Temperature  -  This  temperature  is  read  using  the 
low  level  A/D  card  each  time  torque  is  read  for  the  purpose  of  correct¬ 
ing  torque.  If  housing  temperature  exceeds  a  certain  level,  then  the 
apparatus  is  shut  down.  The  12  readings  taken  within  one  minute  are 
averaged  and  recorded  on  tape  every  ten  minutes. 

Oil  Temperature  -  This  temperature  is  read  one  time  per  minute 
using  the  low  level  A/D  card.  Temperature  is  compared  to  a  limit. 
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The  apparatus  is  shut  down  if  temperature  goes  above  this  limit.  This 
temperature  is  recorded  on  tape  approximately  every  ten  minutes. 

Water  Temperature  -  The  procedure  for  water  temperature  is  the 
same  as  for  oil  temperature.  Water  temperature  control  is  handled  by 
the  Honeywell  controller  described  previously  [6], 

Vibration  Level  -  Vibration  level  is  monitored  continuously  by  a 
vibration  analyzer.  The  signal  is  compared  to  a  set  point  using  an 
electronic  comparator  circuit.  When  vibration  level  exceeds  a  certain 
level,  a  relay  is  closed.  On  manual  operation,  closing  of  the  relay 
shuts  down  the  test  apparatus.  Under  computer  control,  each  time  the 
relay  closes,  the  computer  is  signaled  through  the  event  sense  card. 

If  this  occurs  more  than  five  times  and  for  a  period  greater  than  ten 
seconds,  the  test  apparatus  is  shut  down. 

Plotting  -  Torque,  face  temperature,  and  leakage  rate  are  plotted 
as  functions  of  time  at  one  minute  intervals  throughout  the  test.  An 
example  of  such  a  plot  is  shown  in  Figure  2-7.  This  provides  a  visual 
indication  of  the  test  results  as  the  test  proceeds.  Data  stored  on 
the  tape  can  be  replotted  after  a  test  is  completed. 

Proximity  Probe  -  The  proximity  probe  for  measuring  the  change  in 
gap  at  the  face  of  the  seal  (described  later)  is  read  using  the  high 
speed  A/D  card.  This  provides  a  large  number  of  readings  over  a  short 
period  of  time. 

Surface  Analyzer  -  The  high  speed  A/D  card  in  conjunction  with  the 
programmable  timer  is  used  to  read  a  roughness  signal  for  the  surface 
analyzer.  The  timer  pulses  the  A/D  card  at  a  high  rate  of  speed  so 
that  5  to  10000  points  can  be  quickly  input  into  the  computer  for 
further  analysis.  The  computer  roughness  analysis  is  described  in 
detail  later. 

The  program  performs  its  various  functions  above  by  being  in  an 
infinite  loop.  There  are  several  possible  loops  and  modified  loops 
built  into  the  program  to  perform  various  tasks  at  different  times. 

To  change  from  one  loop  to  another,  certain  keys  are  defined  to  shift 
control.  To  illustrate  some  of  the  features  of  the  program,  the  key 
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Figure  2-7.  Sample  Plot. 


definitions  are  given  in  Table  2-1. 

Operation  of  the  test  apparatus  using  the  computer  control  is 
greatly  simplified  compared  to  using  chart  recorders,  etc.  The  data 
is  corrected  and  converted  automatically  to  the  proper  units.  Having 
all  output  on  the  same  graph  with  a  greatly  compressed  time  scale  is 
quite  helpful  in  interpreting  the  results.  There  is  also  considerable 
flexibility  for  changing  the  type  of  data  collected.  Control  and 
safety  features  are  easily  programmed  in. 

Problems  that  have  arisen  in  the  system  have  been  related  to  com¬ 
puter  hardware.  The  original  input/output  ROM  in  the  computer  had  a 
logic  error  which  caused  the  computer  to  hang  up  at  random.  The  time 
clock  had  a  malfunction  which  caused  it  to  reset  to  zero  at  random. 

One  loose  connector  caused  a  computer  hangup.  Since  these  problems 
have  been  corrected,  operation  has  been  reliable  and  several  tests  have 
been  run  using  the  computer  system. 

Test  Program 

Using  the  apparatus  described,  45  different  test  runs  have  been 
made  to  date.  The  time  duration  for  these  tests  ranges  from  a  few 
hours  to  several  hundred  hours.  The  types  of  tests  run  are  as  follows: 

Flat  face:  Seal  faces  lapped  flat  and  run 

Waviness:  Continuously  moving  waviness  of  differing  amounts 


Wear:  A  test  of  long  duration  to  provide  measurable 

wear 

Pressure:  Pressure  is  varied  from  standard  conditions 

Speed:  Speed  is  varied  from  standard  conditions 

Proximity  Probe:  Proximity  probe  is  used  to  measure  face  displace¬ 
ment 

Standard  or  base  case  test  conditions  have  been  established  as 
follows: 
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Table  2-1 
Key  Definitions 

Key  Function 

KO  "Shut  Down"  -  This  subprogram  trips  main  control  relay  out  and  shuts 
down  the  apparatus.  Program  shifts  to  idle  where  "CONTROL  PROGRAM 
READY"  is  displayed.  Program  can  be  shifted  back  to  the  "start" 
mode. 

K1  "Run"  -  Data  is  taken.  Time  is  set  at  zero.  "Run"  can  be  used 
any  time  during  a  test.  "Run"  will  restart  the  data  at  time  *  0 
and  previous  data  will  be  erased.  Subsequent  shifts  to  "start"  and 
"Conrun"  do  not  shift  time  axis  and  previous  data  are  retained. 

K2  "Start"  -  When  "CONTROL  PROGRAM  READY"  then  "start"  can  be  executed. 
This  closes  the  main  control  relay  so  the  machine  can  be  started  up. 
Pressure  is  controlled. 

K3  "Axis"  -  To  draw  and  set  up  axes  any  time.  "Axis"  subprogram  runs 
in  the  "start"  mode  so  that  pressure  is  controlled  so  that  test  may 
continue.  "Axis"  goes  to  "start"  after  axes  have  been  plotted.  To 
continue  test  use  "Conrun."  Data  which  were  taken  beyond  the  pre¬ 
vious  time  axis  will  be  plotted  first  before  current  data  are  plotted 
Program  will  automatically  do  this  after  an  "axis"  and  "conrun"  are 
used.  Data  can  be  stored  for  only  five  hours  beyond  the  plotter  lim¬ 
it.  Going  beyond  this  will  cause  the  program  to  exceed  dimensional 
limits. 

K4  "Rezero"  -  The  torque  can  be  rezeroed  at  any  time.  Temperature  cor¬ 
rection  is  shifted  accordingly.  Machine  is  placed  in  "CONTROL  PRO¬ 
GRAM  READY"  mode  after  rezeroing.  Test  can  be  started  over  by  using 
"start"  and  run,  or  it  can  be  continued  by  using  "start"  and  "conrun. 

K5  "Conrun"  -  Continues  running  after  a  "start"  interrupt.  Program 
takes  data  from  where  it  left  off.  Missing  data  are  represented 
by  a  gap. 

K6  "Testend"  -  Gives  value  of  current  file  and  stops  taking  data. 

Program  goes  into  "start"  mode.  Finds  torque,  temperature,  and 
leakage  averages.  Use  "shutdown"  to  stop  machine. 
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RPM  =  1800 

Albuquerque  tap  water 

38°C  (100°F)  seal  temperature  environment 

Seal  rings  were  prepared  as  follows: 

Secondary  rotating  ring:  WC,  lapped  and  polished,  0.03  to  0.08 

um  (1  to  3  uin.)  CLA,  radial  taper 
< 100. 10-6 

Primary  stationary  floating  ring:  P658RC  carbon,  as  lapped,  0.51 

um  (20  yin.)  CLA,  radial  taper  <50. 1 0“^ 

The  above  specifications  on  seal  rings  were  held  for  most  of  the 
tests.  However,  in  some  cases  where  previous  wear  was  small,  WC  rings 
were  used  for  more  than  one  test  without  relapping.  New  (as  manufac¬ 
tured)  carbon  rings  had  a  lower  surface  roughness  than  quoted  above 
for  those  prepared  on  our  own  lapping  machine.  Also,  the  radial  taper 
on  carbon  rings  prepared  in  our  own  lab  was  generally  (but  not  always) 
less  than  20. 10"6.  Radial  taper  on  some  of  the  as  manufactured  rings 
was  larger  than  this. 

Seal  roughness  (CLA)  was  measured  in  a  radial  and  tangential 
direction  both  before  and  after  each  test.  Radial  taper  was  measured 
before  and  after  each  test.  For  measuring  carbon  ring  wear,  the  nose 
height  of  the  seal  was  measured  relative  to  the  shoulder  using  the 
surface  analyzer  and  a  gage  pin.  For  later  tests,  the  carbon  was  oven 
dried  after  test  to  minimize  distortion  due  to  water  absorption.  In 
later  tests,  surface  roughness  statistics  were  determined  using  the 
surface  analyzer  in  conjunction  with  the  computer. 

The  installation  procedure  for  the  seal  faces  was  to  wipe  each 
face  with  solvent,  dry  the  faces,  and  then  wet  each  surface  with  water 
before  they  were  mated. 

Appendix  B  contains  tabulated  results  for  all  tests.  Appendix  C 
contains  performance  curves  for  those  tests  conducted  since  the  computer 
system  has  been  installed.  The  detailed  results  for  each  type  of  test 


will  now  be  presented  and  discussed. 

Flat  Face 

Figure  2-8  shows  typical  temperature  and  torque  versus  time 
performance  curves  for  a  flat  faced  seal.  There  was  no  measurable 
leakage  in  this  test  so  it  does  not  appear  on  the  graph.  The  gap  in 
the  curves  was  caused  by  a  data  collection  malfunction;  the  seal  it¬ 
self  ran  continuously.  The  torque  curve  is  characteristic  of  a  flat 
faced  test;  there  is  a  wide  unpredictable  variation  in  torque  level. 
Occasionally,  the  torque  will  increase  abruptly.  Face  temperature 
follows  torque  as  expected. 

Upon  disassembly,  surface  analyzer  traces  of  the  carbon  show  a 
divergent  coning  of  about  1 300-10“^  (see  data  in  Appendix  B).  This 
indicates  that  during  operation  thermal  coning  caused  the  seal  faces 
to  form  a  convergent  taper.  After  operation  under  this  condition, 
the  faces  wear  flat  such  that  when  cooled  they  appear  to  be  divergent. 

Figure  2-8  shows  that  initial  torque  is  relatively  low  and  grad¬ 
ually  increases  with  time.  This  behavior  agrees  with  the  theory  pre¬ 
sented  in  reference  [8]  which  suggests  that  during  initial  operation, 
thermal  coning  produces  a  convergent  film  shape  and  reduces  friction. 
As  the  face  wears  flat,  the  friction  increases. 

Table  2-2  shows  results  for  five  flat  faced  tests  (pgas  =  0). 

No  measurable  leakage  is  shown  for  any  of  the  tests.  Average  fric¬ 
tion  torque  (for  the  final  hours  of  the  test)  varied  widely  from  7.7 
to  13.7  N-m.  Although  the  gradual  increase  in  flat  faced  seal  torque 
can  be  explained  as  mentioned,  the  sudden  changes  in  the  torque  curves 
are  more  difficult  to  explain.  It  is  possible  that  sudden  increases 
are  caused  by  the  faces  drying  out,  i.e.,  the  fluid  becomes  sealed 
off  at  the  outside  diameter.  Dry  carbon  has  a  much  higher  friction 
coefficient  than  wet  carbon.  After  the  temperature  increases,  suffi¬ 
cient  convergent  coning  occurs  such  that  fluid  reenters  the  gap  and 
friction  is  reduced. 
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Table  2-2 


Flat  Face  and  Waviness  Tests 
(Test  durations  greater  than  10  hrs) 
pR  0  =  3.45  MPa  1800  RPM 


pgas 

(MPa) 

Test 

No. 

Test 

Duration 

(h) 

Average 

Face 

Temperature 

(C) 

Average 

Torque 

(N-m) 

Average 

Leakage 

(cm^/min) 

0 

5 

19 

44.4 

mam 

0.0 

0 

13 

24 

46.1 

■9 

0.0 

0 

22 

100 

48.9 

■99 

0.0 

0 

35 

100 

43.9 

8.4 

0.0 

0 

40 

111 

42.8 

9.2 

0.0 

3.45 

41 

102 

41.1 

4.6 

0.0 

4.14 

29 

25 

42.8 

5.1 

0.0 

4.76 

23 

100 

40.6 

2.0 

0.3 

4.83 

27 

21 

40.0 

2.9 

0.7 

4.83 

36 

100 

40.0 

2.4 

1.4 

5.03 

12 

24 

41.1 

2.3 

0.6 

6.90 

21 

100 

40.0 

1.5 

1.1 

6.90 

25 

43 

40.0 

2.5 

0.6 

6.90 

32 

18 

40.0 

3.1 

0.7 

6.90 

34 

100 

39.4 

1.4 

1.4 

6.90 

38 

42 

39.4 

2.0 

2.5 

6.90 

39 

100 

39.4 

1.6 

3.2 

7.07 

10 

24 

40.0 

2.9 

1.2 

8.27 

44 

73 

38.3 

2.1 

8.2 

8.52 

11 

23 

40.0 

2.3 

5.2 

Waviness 

Figure  2-9  shows  a  typical  test  result  for  a  wavy  seal.  At  the 
start  of  the  test  the  torque  is  high  and  then  levels  off  to  a  lower- 
cyclic  value.  The  initial  high  torque  is  caused  by  the  fact  that  it 
takes  several  hours  for  the  carbon  face  to  assume  its  equilibrium 
shape.  Once  this  occurs,  then  the  torque  cycle  becomes  constant  with 
increasing  time. 

Like  torque,  leakage  stabilizes  at  some  constant  cyclic  value. 

The  period  of  the  leakage  cycle  is  1.58  h.  This  corresponds  to  1/6 
of  the  time  for  one  revolution  of  the  waviness  cylinder  (9.47  h). 

The  torque  and  leakage  variations  are  being  caused  by  a  superposition 
of  the  third  harmonic  waviness  over  a  second  harmonic  waviness  which 
is  being  produced  by  the  two  drive  forces.  Thus,  the  significance  of 
the  cyclic  behavior  is  that  operation  can  be  influenced  by  unwanted 
deflections,  and  that  such  deflections  must  be  controlled.  Otherwise, 
the  extremes  of  the  cycle  could  represent  unacceptable  operating  con¬ 
ditions.  In  the  case  shown,  the  extremes  are  not  large.  Additional 
waviness  test  curves  are  shown  in  Appendix  C. 

The  greatest  difference  between  the  flat  face  and  wavy  test  re¬ 
sults  is  that  the  torque  and  temperature  of  the  wavy  seal  are  much 
lower  than  for  the  flat  face  seal.  In  Table  2-2  the  results  of  a 
number  of  waviness  tests  at  different  gas  pressures  are  shown.  The 
data  are  averages  over  the  last  9.47  h  of  the  test.  The  effect  of 
waviness  on  torque  is  large.  Torque  values  are  reduced  to  one  fourth 
of  those  values  for  the  flat  face. 

For  leakage,  the  trend  is  as  expected  with  leakage  increasing  from 
unmeasurable  quantities  in  the  case  of  the  flat  face  to  several  mJt/min 
at  the  higher  gas  pressures.  The  scatter  in  results  is  greater  than  for 
torque.  One  reason  for  this  is  that  some  of  the  earlier  leakage  measure 
ments  were  made  using  a  clock  and  graduated  cylinder  before  the  comput¬ 
erized  apparatus  was  completed.  It  is  likely  that  this  approach  intro¬ 
duced  some  error  because  there  was  no  way  to  average  over  a  complete 
leakage  cycle  as  shown  in  Figure  2-9.  In  addition,  since  leakage  is 
related  to  film  thickness  cubed,  very  small  differences  in  operating 
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equilibrium  film  thickness  lead  to  large  variations  in  leakage. 

Wear 

Table  2-3  shows  the  results  of  several  100  h  wear  tests. 

Torque,  leakage,  and  temperature  values  fit  in  with  the  trends  pre¬ 
viously  described.  Wear  results  require  some  explanation.  First, 
looking  at  the  flat  face  results,  tests  35  and  40  give  comparable 
results  whereas  test  22  shows  a  very  high  wear.  Comparing  torque 
and  temperature  of  35  and  40  to  22  shows  that  test  22  operated  at  a 
significantly  higher  temperature  and  torque.  Thus,  a  different  type 
of  equilibrium  was  established  such  that  wear  was  much  larger. 

Several  of  the  tests  show  a  negative  wear.  The  significance  of 
this  is  that  wear  was  less  than  the  distortion  in  the  3  mm  (0.125  in.) 
nose  of  the  carbon  caused  by  water  absorption  and  thermal  effects. 

This  result  occurred  in  spite  of  careful  drying  procedures.  Wear  mea¬ 
surements  on  a  given  seal  are  repeatable  so  that  this  was  not  the  source 
of  error.  It  is  known  that  strains  of  the  order  of  0.001  can  be  intro¬ 
duced  in  carbon  by  water  absorption.  For  the  3  mm  nose  height,  this 
would  correspond  to  a  dimensional  change  of  3  um  (125  yin.).  Thus, 
for  now  it  is  concluded  that  the  actual  wear  magnitude  for  the  cases 
giving  a  negative  reading  are  of  the  same  order  as  the  nose  distortion 
or  3  um  (125  yin.).  Based  on  Table  2-3,  the  wear  rate  for  pgas  =6.9 
MPa  wavy  cases  are  from  four  to  forty  times  lower  than  the  flat  face 
test. 

To  help  in  narrowing  this  range,  a  series  of  500  h  duration  tests 
are  being  run  to  attempt  to  establish  a  more  repeatable  wear  rate, 
particularly  for  the  waviness  cases. 

Radial  Taper 

In  order  to  test  the  validity  of  the  hydrostatic  theory  of  Refer¬ 
ence  [8]  which  considers  hydrostatic  pressure,  load  sharing,  and 
thermal  coning,  three  tests  (24,  42,  43)  were  run  where  a  converging 
radial  taper  was  lapped  into  the  seal  before  the  test. 

To  create  the  radial  taper,  the  device  shown  in  Figure  2-10  was 
inserted  into  the  back  of  the  seal,  pressurized,  and  the  seal  lapped 
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Figure  2-10.  Radial  Taper  Device 
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flat.  The  device  is  designed  to  create  a  uniform  moment  about  the 
ring  centroid  causing  the  faces  to  diverge.  After  lapping  and  releas¬ 
ing  the  pressure,  a  converging  taper  results.  Taper  is  proportional 
to  the  pressure  applied. 

Table  2-4  shows  the  results  from  three  radial  taper  tests  and  one 
flat  faced  test  for  comparison.  Torque  in  the  radial  taper  cases  is 
less  than  one  half  of  that  in  the  flat  faced  case.  Thus,  the  converg¬ 
ing  radial  taper  produces  a  larger  component  of  hydrostatic  lifting 
force  than  in  the  flat  faced  case.  Test  number  43  shows  that  torque 
is  insensitive  to  the  amount  of  taper.  However,  it  can  be  shown  for 
a  balanra  ratio  of  1.0;  this  result  can  be  expected  in  this  range  of 
large  radial  tapers. 

Figure  C-5  (test  42)  of  Appendix  C  shows  a  gradual  decrease  in 
torque  with  time  as  wear  progresses  across  the  face  of  the  seal. 

Also,  there  is  some  leakage  initially  which  gradually  diminishes  to 
unmeasurable  values  with  time.  The  torque  and  temperature  are  quite 
steady  compared  to  a  flat  faced  test  except  for  an  occasional  excur¬ 
sion  such  as  at  184  and  264  h  for  test  42  (Figure  C-5). 

Figure  2-11  shows  a  comparison  of  the  original  taper  to  the  worn 
taper  after  the  test.  It  is  clear  that  a  certain  fraction  of  the  seal 
face  has  worn  flat.  This  causes  a  change  in  operating  performance 
with  time.  Table  2-4  shows  that  a  greater  fraction  of  the  seal  is 
worn  flat  for  a  longer  duration  test.  This  result  suggests  that  after 
a  long  period  of  operation  the  taper  would  wear  off  completely  and  the 
seal  would  continue  operation  performing  like  a  flat  faced  seal.  A 
comparison  of  these  results  to  theory  will  be  made  in  a  later  report. 

Pressure 

Table  2-5  shows  the  effects  of  varying  sealed  pressure  for  a  wavy 
seal.  At  the  lowest  test  pressure,  leakage  is  large  in  comparison. 
This  occurs  because  the  wave  is  not  flattened  out  as  much  as  at  higher 
sealed  pressures.  At  the  highest  pressure,  torque  is  significantly 
larger  than  at  the  lower  sealed  pressures.  This  indicates  that  a  sig¬ 
nificant  fraction  of  the  face  load  is  being  carried  by  mechanical 
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Radial  Profile  (.001  filter) 
Test  #42  (Initial  Trace) 


Radial  Profile  (.001  filter) 
Test  #42  (Final  Trace) 


Figure  2-11.  Radial  Taper  Profile,  Before  and  After  Test. 
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Table  2-5 

Sealed  Pressure  Effects  -  Wavy  Seal 
pgas=6.9Mpa  1800  RPM 


Test 

Number 

— 

ph2o 

MPa 

Test 

Duration 

h 

Average 

Face 

Temperature 

°C 

Average 

Torque 

N-m 

Average 

Leakage 

cm^/tni  n 

9 

1.73 

25 

40 

2.2 

11 .0 

21 

3.45 

100 

40 

1.5 

1.1 

34 

3.45 

100 

39 

1.4 

1.4 

39 

3.45 

100 

39 

1 .6 

3.2 

17 


5.18 


32 
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contact.  However,  this  torque  is  about  one  half  of  that  of  the  flat 
faced  seal.  Results  for  test  number  18  in  Appendix  B  shows  that  by 
increasing  the  gas  pressure,  the  friction  torque  can  be  reduced  to 
one  half  of  that  for  case  17  above  while  leakage  increases  to 
6  cm3/min.  In  general,  this  behavior  is  as  expected.  Some  modifi¬ 
cations  of  the  results  in  Table  2-5  would  be  expected  if  the  test 
duration  were  made  the  same  for  all  tests. 

Speed 

Table  2-6  shows  two  sets  of  results  where  sealed  pressure  and 
waviness  pressure  are  held  constant  while  speed  is  varied,  one  set  of 
data  at  a  higher  waviness  pressure  than  the  other.  A  comparison  is 
made  somewhat  difficult  because  the  test  durations  are  different,  how¬ 
ever,  certain  observations  can  be  made.  Considering  the  results  for 
Pgas  =6.9  MPa,  torque  increases  with  increasing  speed.  Leakage 
appears  to  increase,  although  this  is  not  definite. 

The  trend  in  torque  values  is  exactly  opposite  to  what  we  would 
expect  assuming  hydrodynamic  mechanisms  are  operative.  Hydrodynamic 
theory  predicts  a  reduced  torque  at  higher  speed,  at  least  over  some 
range  of  speeds  [  1  ].  The  explanation  for  this  apparent  contradiction 
in  behavior  is  complex  and  is  taken  up  in  detail  in  Chapter  3. 

The  data  for  pgas  =  4.8  to  5.5  MPa  show  no  particular  trend  one 
way  or  the  other.  It  is  thought  that  the  low  speed  results  particularly 
may  be  influenced  measurably  compared  to  the  others  by  the  seal  not 
being  worn  in  at  its  particular  operating  conditions.  Although  large 
time  tests  are  needed  to  refine  the  data,  it  is  not  expected  that  the 
above  noted  trend  would  change. 

Proximity  Probe 

Figure  2-12  shows  the  installation  of  a  Kaman  Science  KD-2300- 
.5SU  displacement  measuring  probe.  The  probe  operates  based  on  imped¬ 
ance  variations  caused  by  eddy  currents  induced  in  a  conductive  target. 

A  change  in  displacement  of  the  carbon,  to  which  the  probe  is  bonded, 
relative  to  the  WC  face,  produces  a  change  in  DC  output  of  the  probe 
driver.  This  signal  variation  was  read  using  the  high  speed  A/D  card 
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Table  2-6 


Speed  Effects  -  Wavy  Seal 
PH  q  =  3.45  MPa 

Average 

Face 

Test  pgas  Duration  Temp. 

Number  Speed  MPa  h  °C 

26  500  6.9  4  40 

33  900  6.9  20  40 

34  1800  6.9  100  39 

39  1800  6.9  100  39 

21  1800  6.9  100  40 

30  2700  6.9  46  41 

31  3600  6.9  6  43 


28  900  4.8  5  40 

36  1800  4.8  100  40 

19  2700  5.5  24  40 
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described  previously.  Sensitivity  of  the  probe  was  found  to  be  1.4 
volts/mm  so  that  a  1  urn  displacement  change  caused  a  1.4  mv  change 
in  signal . 

Although  the  WC  used  in  these  experiments  is  slightly  magnetic, 
it  was  found  that  Kaman's  probe  KD-2300-.5SUM  for  magnetic  materials 
would  not  operate.  The  probe  referenced  above  for  non-magnetic  mate¬ 
rials  operated  successfully. 

To  operate  satisfactorily  the  probe  and  its  cable  must  be  pro¬ 
tected  from  water.  To  do  this,  the  probe  was  epoxied  into  the  carbon 
and  a  3.2  mm  OD  brass  tube  was  epoxied  at  the  rear  of  the  seal  and 
sealed  at  the  vessel  housing.  This  provided  a  completely  water  tight 
passage  for  the  probe  cable.  The  objective  of  the  experiments  was  to 
measure  the  relative  change  in  displacement  between  the  carbon  and 
the  WC  that  occurs  as  the  waviness  cylinder  is  rotated.  Such  data 
can  be  used  to  verify  deflection  theory  which  is  an  essential  part 
of  modeling  the  seal  behavior.  No  attempt  was  made  to  establish  a 
zero  displacement  reading  or  to  measure  actual  film  thickness.  Only 
relative  displacements  were  obtained. 

The  procedure  used  was  to  drive  the  waviness  cylinder  at  its 
high  speed,  about  two  minutes  per  revolution.  A  switch  was  connected 
to  the  drive  system  so  as  to  trigger  the  computer  when  exactly  one 
turn  of  the  waviness  cylinder  was  completed.  Voltage  output  at  equal 
time  steps  was  read  and  stored  in  computer  memory.  About  500  voltage 
readings  were  taken  per  turn  of  the  waviness  cylinder.  Figure  2-13 
shows  a  typical  output  signal  for  one  turn  of  the  waviness  cylinder. 
There  is  a  considerable  amount  of  noise  on  the  signal.  A  Fourier 
analysis  of  the  signal  filters  out  the  noise  and  shows  that  the  sig¬ 
nal  has  predominant  third  and  sixth  harmonic  components  as  expected. 
Table  2-7  shows  the  Fourier  analysis  results  of  all  of  the  various 
proximity  probe  tests. 

Tests  were  run  with  the  seal  under  spring  load  only  and  under 
water  pressure  and  spring  load  combined  and  at  various  gas  pressures. 
Results  show  that  the  waviness  is  repeatable  and  that  water  pressure 
.causes  considerable  flattening  of  the  wave.  The  tests  were  performed 
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Figure  2-13.  Output  From  Proximity  Probe. 


on  an  unworn  carbon  ring.  A  comparison  of  these  results  to  theory  is 
made  in  Chapter  3. 

The  original  plan  was  to  measure  waviness  and  particularly  the 
change  in  waviness  during  an  operating  waviness  test.  This  was  not 
done  however  because  of  a  bearing  seizure  of  the  waviness  cylinder 
caused  by  the  higher  speed  operation  of  the  waviness  cylinder.  The 
repair  of  this  failure  required  that  the  proximity  probe  be  destroyed. 
These  tests  may  be  performed  at  a  later  time  after  the  proximity  probe 
is  replaced. 

Wear  Profile 

After  each  wear  test  a  wear  profile  of  the  carbon  surface  was  made 
using  the  surface  analyzer.  Three  such  profiles  are  shown  in  Figure 
2-14.  These  profiles  are  for  p^O  =  3.45  MPa,  Pgas  *  6.9  MPa,  1800 
RPM,  100  h  tests.  The  significant  feature  of  all  three  wear  profiles 
is  that  both  the  inside  and  outside  corners  of  the  seal  face  are  worn 
tapered.  The  approximate  amount  of  taper  is  shown  in  the  figure.  The 
tapering  on  the  third  case  is  not  nearly  as  pronounced  as  in  the  first 
two  cases.  There  is  no  apparent  explanation  for  this.  It  should  be 
noted  that  wear  and  taper  of  the  WC  mating  face  was  insignificant  rela¬ 
tive  to  the  taper  shown. 

The  significance  of  the  tapering  on  both  corners  is  that  the  seal 
face  must  rock  or  tilt  from  a  converging  configuration  to  a  diverging 
configuration  as  the  waviness  cylinder  passes  by.  Thus,  from  these 
results  and  the  proximity  probe  results,  it  is  concluded  that  the 
waviness  producing  device  creates  waviness  as  well  as  an  alternating 
tilt  around  the  seal.  The  faces  do  not  remain  radially  parallel  as 
was  assumed  in  the  hydrodynamic  theory  contained  in  Reference  [  1  ]. 

A  full  explanation  of  this  type  of  deformation  is  contained  in 
Chapter  3. 

Surface  Roughness 

Surface  roughness  was  measured  before  and  after  each  test.  These 
data  are  shown  in  Appendix  B.  After  a  test  the  roughness  of  the  carbon 
ranges  from  0.2  to  0.5  um  CLA.  The  WC  is  generally  around  0.05  urn  CLA. 
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After  studying  the  roughness  values,  no  particular  pattern  or  relation¬ 
ship  with  other  variables  seem  apparent.  However,  after  more  data  has 
been  accumulated,  a  closer  examination  is  needed. 

Since  the  time  the  computer  was  installed  in  the  system,  it  has 
been  possible  to  perform  a  statistical  analysis  on  data  obtained  in  the 
surface  roughness  analysis.  To  do  this  the  surface  analyzer  output  is 
connected  directly  to  the  high  speed  A/D  card  described  previously.  As 
the  surface  analyzer  makes  its  trace  across  the  surface,  several  thou¬ 
sand  ordinate  values  are  read  and  stored  into  the  computer.  A  computer 
program  has  been  written  to  take  these  points  and  calculate  CLA,  stan¬ 
dard  deviation,  skewness,  Kurtosis,  a  histogram,  and  an  autocorrelation 
function.  These  data  provide  much  more  information  about  the  surface 
than  is  provided  by  the  CLA  reading  given  by  the  surface  analyzer. 

Examples  of  the  output  data  are  given  in  Figures  2-15  and  2-16 
for  the  carbon  ring  of  test  number  41.  The  histogram  is  characteris¬ 
tic  of  a  surface  having  a  series  of  valleys  with  peaks  being  chopped 
off,  typical  of  carbon.  The  negative  skewness  characterizes  the  shape 
of  the  histogram  shown.  The  autocorrelation  function  shown  in  Figure 
2-15  decays  toward  zero  showing  that  the  roughness  has  no  periodicity, 
which  is  expected.  Characteristic  asperity  size  is  taken  at  R(x)  =  0.5 
which  corresponds  to  9.5  urn. 

Table  2-8  shows  a  sunmary  of  results  from  several  statistical 
analyses  of  the  typa  described.  Test  40  was  a  flat  face,  whereas 
test  41  was  a  Pgas  *  3.45  MPa  waviness  test.  Comparing  test  40  at  0°  to 
test  41  at  0°  shows  no  significant  difference  in  the  surface  statistics. 
This  suggests  that  the  wear  processes  may  be  similar,  even  though  seal 
performance  was  greatly  different. 

The  carbon  ring  from  test  number  41  was  analyzed  at  three  differ¬ 
ent  angular  locations  to  see  if  there  would  be  any  difference  in  surface 
characteristics  dependent  upon  where  within  one  wave  the  sample  was 
taken.  Again,  there  appears  to  be  no  significant  difference.  Thus, 
once  the  carbon  is  worn  in,  surface  statistics  do  not  greatly  change 
as  a  function  of  whether  or  not  contact  has  been  recent  or  not. 
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Figure  2-16.  Autocorrelation  Function  -  Test  41. 
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Table  2-8 

Carbon  Surface  Statistics 


Test  No. 

c 

(ym) 

Skewness 

*0.5 

(ym) 

Location 

Radial 

Tangential 

Radial 

Tangential 

Radial 

Tangential 

O 

o  o 

0.56 

0.47 

-2.95 

-2.48 

13.4 

9.9 

41 

0° 

0.59 

0.45 

-2.40 

-2.49 

9.5 

10.2 

41 

40° 

0.60 

0.43 

-1.94 

-2.68 

9.1 

7.4 

41 

80° 

0.76 

0.59 

-2.56 

-2.86 

12.0 

10.2 

Avg 

0.63 

0.49 
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One  significant  feature  can  be  noted  from  Table  2-8.  In  all  cases, 
the  standard  deviation  is  lower  when  measured  in  the  tangential  direc¬ 
tion  than  when  measured  in  the  radial  direction.  This  may  be  caused 
by  characteristic  tangential  grooving  of  the  carbon  surface.  However, 
the  roughness  in  both  directions  is  similar  so  that  the  seal  surface 
can  be  considered  as  isotropic  for  analysis  purposes. 

At  the  present  time,  only  a  limited  number  of  seal  statistical 
analyses  such  as  above  have  been  made.  As  more  data  becomes  available, 
it  is  expected  that  these  data  will  be  useful  in  comparing  surfaces 
generated  by  wear  in  different  types  of  sealing  environments. 
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Chapter  3 

Thermoelastic  Instability  in  Face  Seals 


As  a  part  of  an  overall  interest  in  understanding  mechanical  seal 
fundamentals,  some  effort  has  been  spent  in  attempting  to  better  under¬ 
stand  the  role  of  thermoelastic  instability  in  face  seals.  Interest  in 
this  phenomenon  is  heightened  given  the  opportunity  to  do  testing  using 
the  experimental  seal  test  apparatus  described  previously. 

Thermoelastic  instability  in  a  face  seal  occurs  when  unbounded  cir¬ 
cumferential  pressure  and  temperature  variations  occur  in  an  otherwise 
parallel  face  contact.  Previous  investigators  have  shown  [9,10,11]  that 
for  a  given  material  pair  and  geometry,  there  is  a  certain  critical 
speed  of  sliding  above  which  such  pressure  and  temperature  variations 
grow.  This  is  called  the  threshold  of  instability  or  critical  speed. 
Below  this  critical  speed,  such  variations  will  diminish  with  time. 

Given  this  basic  knowledge,  the  theory  of  thermoelastic  instabil¬ 
ity  has  been  applied  to  various  mechanical  configurations  in  an  attempt 
to  accurately  predict  critical  speed,  and  predictions  have  been  compared 
to  experimental  results.  Mechanical  face  seals  represent  one  applica¬ 
tion  where  the  prediction  of  thermoelastic  instability  is  of  great  in¬ 
terest  to  designers. 

Netzel  [12]  has  reported  evidence  of  patch  contact  on  the  hard 
faces  of  two  high  speed  anhydrous  anmonia  seals  and  one  water  seal. 
Existing  theory  [9]  based  on  an  insulator  sliding  on  a  conducting  body 
shows  that  the  critical  speed  in  these  applications  is  much  lower  than 
reality.  Theory  [9]  based  on  two  conducting  bodies  gives  critical 
sliding  speeds  higher  than  those  observed. 

In  order  to  more  accurately  predict  critical  speeds  in  mechanical 
seals,  a  modified  theory  of  thermoelastic  instability  for  mechanical 
face  seals  is  developed  herein.  The  effects  of  ring  deflection  and 
heat  transfer  between  the  rings,  heretofore  neglected,  are  included. 

It  is  shown  that  at  lower  numbers  of  waves,  ring  deflection  predominates 
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and  completely  changes  the  character  of  the  instability  relationships. 
Thermal  contact  resistance  between  the  faces  acts  to  decrease  the  crit¬ 
ical  sliding  speed.  A  similar  approach  was  taken  in  reference  [13]. 
However,  certain  assumptions  made  concerning  heat  transfer  led  to  re¬ 
sults  which  are  inconsistent  with  those  of  Burton  [9].  In  this  paper, 
these  assumptions  are  avoided,  and  the  results  are  consistent  with  those 
of  Burton. 

Theory 

Figure  3-1  shows  the  basic  mechanical  face  seal  configuration. 

As  one  ring  rotates  relative  to  the  other,  sliding  occurs  at  the  seal¬ 
ing  interface.  The  faces  are  loaded  against  each  other  by  spring  and 
fluid  pressures.  Both  seal  rings  are  assumed  to  be  of  thickness  h  and 
mean  radius  R. 

Figure  3-2  shows  the  faces  of  the  seals  developed.  Ring  H  is 
rotating  at  angular  speed  w  and  ring  S  is  fixed.  It  is  assumed  that 
there  exists  a  sinusoidal  temperature  disturbance  on  face  S  moving  rel¬ 
ative  to  face  S  at  angular  speed  $.  It  is  also  assumed  that  a  differ¬ 
ent  sinusoidal  temperature  disturbance  occurs  on  face  H  also  moving  at 
speed  $.  There  are  n  periods  or  waves  around  the  seal,  and  the  condi¬ 
tion  n  >  2  must  be  satisfied  for  the  ring  geometry. 

Relative  to  the  fixed  ring,  the  temperature  variations  are  assumed 
as  follows 

T^(e,t)  =  T^  sin  n(e  -  <j>t)  ,  (3-1) 

TH(e,t)  =  Th  sin  n(e  -  <|>t  +  n)  ,  (3-2) 

where  the  angle  n  allows  for  a  phase  shift  between  the  two  temperature 
waves.  Relative  to  the  rotating  ring,  TH(e,t)  becomes 

TH(e,t)  =  Th  sin(e  +  8t  +  n)  .  (3-3) 


where 


$  +  S  s  ui  . 


(3-4) 


57 


Figure  3-1.  Mechanical  Face  Seal 


H 


(tf 


Figure  3-2.  Seal  Rings. 


Given  the  sinusoidal  temperature  distribution  at  the  face,  the 
temperature  distribution  within  the  ring  can  be  found.  With  reference 
to  Figure  3-3,  it  is  assumed  that  the  ring  is  thin  such  that  radial 
temperature  variations  can  be  neglected.  It  is  also  assumed  that 
there  is  heat  transfer  to  surrounding  fluid  media  as  shown.  Thus, 
the  governing  equation  for  heat  transfer  becomes 


32T  32T  1 

3y2  302  ? 


(T 


11  1 
3t  k  ' 


(3-5) 


Given  the  form  of  the  face  temperature,  equation  (3-1)  or  (3-3),  the 
exact  solution  for  the  above  equation  for  a  finite  length  seal  ring  is 
difficult  to  find,  so  the  solution  for  a  semiinfinite  body  is  used 
where 


gy  0  as  y  •*  °°  .  (3-6) 

Similar  to  that  given  by  Burton  [1],  the  solution  for  body  S  becomes 

"h<-yc 

Ts(e,t,y)  =  Tse  sin  n(e  -  $t  +  a$ys)  ,  (3-7) 

where 
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It  has  been  found  that  e  <<  1  for  most  values  of  encoun¬ 
tered  in  this  investigation.  That  is,  the  temperature  disturbance  is 
very  shallow.  However,  for  n  *  2  or  3,  in  some  cases  b^  is  as  small 
as  unity,  so  some  error  is  introduced.  Nevertheless,  the  temperature 
distribution  given  by  equation  (3-7)  for  a  semiinfinite  thin  plate  will 
be  considered  as  a  reasonable  first  approximation  for  a  finite  length 
thin  plate  or  cylinder. 

The  signs  on  the  square  roots  in  equations  (3-8)  and  (3-9)  were 
selected  to  give  real  roots  with  the  proper  phase  shift  for  the  assumed 
direction  of  <t>. 

For  body  H,  the  solution  for  equation  (3-5)  gives 


-buy  u 

THU,t,y)  =  THe  sin  n(8  -  <j>t  +  n  +  a^) 
relative  to  the  fixed  reference. 


Heat 

q  = 


transfer  per  unit  area  into  the  faces  is  given  by 


v  aT 
‘K  ay 


y=0 


(3-11) 


(3-12) 


(3-13) 


(3-14) 


(3-15) 
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or 


q$  =  -KjT^C-  sin  n(e  -  4>t)  +  na^  cos  n(e  -  t ) ]  ,  (3-16) 

qH  =  -K^C-  bH  sin  n(e  -  $t  +  n)  +  na^  cos  n(e  -  $t  +  n)]  (3-17) 

both  relative  to  the  fixed  reference. 

As  shown  in  Figure  3-2,  it  is  assumed  that  heat  is  transferred 
between  the  seal  faces  by  conduction.  This  conduction  represents 
thermal  contact  resistance.  It  is  also  assumed  that  all  of  the  heat 
generated  by  friction  occurs  midway  between  the  faces  and  is  conducted 
to  the  two  faces  through  the  thermal  contact  resistance  Uc. 

Calling  the  temperature  at  the  midway  point  T^,  then 

(TM  '  TS)Uc  =  qS  ’  (3-18) 

(Tm  -  Th)Uc  =  qH  .  (3-19) 


Eliminating  T^, 

(TH  -  TS>Uc  ■  "S  -  qH  •  (3-20> 

Upon  substituting  for  Ts,  T^,  q^,  and  q^  from  equations  (3-1),  (3-2), 
(3-16),  and  (3-17)  and  separating  the  sine  and  cosine  terms,  the  fol¬ 
lowing  expressions  are  obtained 

(K^  +  Uc)Th  cos  nn  +  KHTHnaH  S1n  nn  "  TS^KSbS  +  uc)  =  0  »  (3_21) 
-  KHTHnaH  cos  nn  +  (K^  +  uc)th  sin  nn  +  KSTSnaS  =  0  ’  (3-22) 


From  these  expressions. 


(Uc  +  Ksbs)KHnaH  -  (U£  +  KHbH)Ksnas 

tan  nn  =  - * -  , 

(Uc  +  Ksbs)(Uc  +  KHbH)  +  KhICsR^s 


(3-23) 
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(3-24) 


w 


11  i  11  -i 


][h  =  Uc  +  KSbS _ r__ 

"  (Uc  +  KHbH)  cos  nn  +  K^na^  sin  nri  ’ 


Thus,  given  the  material  properties  and  speeds  <p  and  e  and  using  equa¬ 
tions  (3-8),  (3-9),  (3-12),  and  (3-13)  to  solve  for  a^,  b^,  aH  and  bH> 
then  the  phase  shift  n  and  the  temperature  ratio  can  be  determined. 

The  total  heat  generated  at  the  interface  must  be  transferred  into 
the  faces.  Friction  heat  at  the  interface  is  proportional  to  the  inter 
face  pressure  p  and  the  sliding  speed.  Thus  an  energy  balance  gives 

pyR(<j>  +  g)  -  qH  +  qs  .  (3-25) 


Equation  (3-25)  provides  an  expression  for  p  in  terms  of  T$  and  T^. 

Considering  now  the  thermal  deflection  that  occurs  due  to  the 
temperature  distribution  in  the  rings,  it  will  be  assumed  that  for  low 
harmonic  numbers  n  ring  deflection  theory  is  applicable  whereas  for 
high  harmonic  numbers  n,  semiinfinite  plate  theory  is  applicable.  It 
should  be  noted  that  at  all  harmonic  numbers  the  semi  infinite  plate 
theory  gives  larger  thermal  deflections  than  the  ring  theory. 

Using  ring  theory  presented  in  reference  [14],  and  with  reference 
to  Figure  3-4  for  ring  S,  it  can  be  shown  that 


M 


x 


ESJxS 


M 


0 


GSJeS  /  , 
R 


v^ 

R 


) 


+  M;  =  o  , 

-  mx  =  °  , 


(3-26) 

(3-27) 

(3-28) 

(3-29) 
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Figure  3-3.  Seal  Cross  Section. 


Figure  3-4.  Geometry  for  Ring  Deflec 
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where  the  primes  denote  derivatives  with  respect  to  e,  and  v  and  y 
represent  the  face  displacement  and  rotation,  respectively.  Using 
equation  (3-7)  for  the  temperature  distribution  in  the  ring,  equa¬ 
tions  (3-26),  (3-27),  (3-28)  and  (3-29)  can  be  solved  simultaneously 
to  give  the  displacement  v  which  represents  the  thermal  deflection 


STS  =  RasTs[fls(n)  sin  n(e  -  $t)  +  f2$(n)  cos  n(e  -  4>t) ]  ,  (3-30 


where 


flS(n) 


(!i .  y  )e-Vs 

— — ^ — S  ^  2 -  (-  b$  cos  na$ys  +  nas  sin  nasy$) 


JxS(n  -  !)  bS+naS 


"Vs 


(4 *  "24)‘ 


f7b2  -  n2a| )  cos  nasys  -  2nasb$  sin  na$ysJ 


(3-31 


(*  -  *) 


f2S(n)  =  1  ,  2  n  )  s*  ,  ('  ^  $in  n3SyS  "  naS  C0$  naSyS} 

JxS(n  •  1}  (  bS  +  n  aS 


(bs  +  "2a!)‘ 


|^(b2  -  n2aj:)  sin  na$ys  +  2na$bs  cos  nasys] 


(3-32 


For  body  H, 


=  RaHTH[f1H(n)  sin  n(e  -  <*>t  +  n)  +  f2H(n)  cos  n^0  *  +  ’ 


Hi 


where  the  expressions  for  and  f£H  are  the  same  as  equations  (3-31) 
and  (3-32)  except  for  the  subscripts. 

For  the  semi  infinite  plate  assumption,  based  on  the  work  of 
Burton  [9]  it  can  be  shown  that 


flS(n)  R(t>  ’ 

(3-34) 

2k, 

f2S^n^  =  nR$  ^bS  "  * 

(3-35) 

f  /.%  .  2kHaH 
^ih  n  Re  » 

(3-36) 

f2H^  =  '  nRe  ^bH  "  ^  • 

(3-37) 

For  elastic  deflection  it  is  again  assumed  for  small  n  that  ring 
theory  applies  and  for  large  n  semiinfinite  plate  theory  applies.  Based 
upon  previous  work  [13]  and  assuming  that  face  pressure  is  sinusoidal, 
elastic  deflection  is  given  by 

6PS  =  E^  f3S^  ’  (3-38) 


6PH  =  f3H^  ’ 


where  for  ring  theory 


n2  4-  ESJxS 
f  (n  1  -  R3h  GSJ9S 
3S(  }  "  JxS  n2(n2-l)2 


and  for  the  semi  infinite  plate 


f3$(n) 


2 

n  ’ 


(3-39) 


(3-40) 


(3-41) 
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f3H  is  obtained  by  changing  the  subscripts  in  equations  (3-40)  and  (3-41). 

In  order  for  the  rings  to  maintain  contact,  the  sum  of  the  thermal 
and  elastic  deflections  must  be  zero  (assuming  no  wear). 


0  =  6TH  +  6ts  -  SpH  -  <sps  . 


(3-42) 


Equations  (3-30),  (3-33),  (3-38)  and  (3-39)  are  substituted  into  equa¬ 
tion  (3-42).  Pressure  p  is  eliminated  using  equations  (3-25),  (3-16) 
and  (3-17).  Then  for  the  resulting  equations  to  be  valid  in  general, 
the  coefficients  of  the  sin  n(e  -  <j>t)  and  cos  n(e  -  $t)  terms  must 
vanish  independently.  This  gives 


uU  +  S)R 


th  T 

aSflS^  +  —  aHflH(n)  C0S  nn  '  r  aHf2H(n)  Sin  nn 

-  s  s 


TH  Th 

Ksbs  +  —  K^b^  cos  nn  +  |-  KHnaH  sin  nn 
*  S  S 


f3H^  ,  f3S^ 


E  1  ’ 

H 


(3-43) 


u  ( <J>  +  6)R 


H  H 

aSf2S^n^  +  ~aHflH^  sin  nn  +  r  aHf2H(n)  C0S  nn 
.5  S 


TH  T 

-Ksnas  +  y-  Kj^bj_j  sin  nn  -  y2-  KHnaH  cos  nn 

S  S  J  L  H 


rf3H(n)  +  f3S^  "I 


E  1  ’ 
tS 

(3-44) 


*  U 

where  ^  and  n  are  given  by  equations  (3-23)  and  (3-24). 
*S 


Equations  (3-43)  and  (3-44)  represent  two  equations  in  two  un¬ 
knowns  <f>  and  6.  Since  only  the  ratio  T^/T^  appears,  T^  and  T<.  can 
take  any  arbitrary  value.  Thus,  solutions  to  these  equations  repre¬ 
sent  critical  speeds  where  the  system  is  on  the  threshold  of  insta¬ 
bility.  Numerical  methods  are  required  for  a  solution. 

To  provide  a  check  on  the  above  equations,  the  semi  infinite  body 
assumption  was  used  and  Uc  -*•  «  and  =  0.  These  conditions  make  the 
above  equations  conform  precisely  to  those  of  Burton  [9],  Using  the 
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above  equations,  solutions  were  carried  out  for  graphite  running  against 
cast  iron  as  in  8urton’s  [9]  Table  III.  The  results  from  the  above 
equations  agreed  with  those  of  Burton. 

Results 

Netzel  [12j  found  evidence  of  hot  spots  in  three  different  high 
speed  seals  using  a  combination  of  carbon  versus  tungsten  carbide.  The 
data  and  observations  are  given  in  Table  3-1. 


Table  3-1. 

Hot  Spot 

Observations 

Fluid 

Pressure 

(psi) 

Speed 

(RPM) 

Shaft  Size 
(in.) 

No.  of 
Spots 

Anhydrous 

Ammonia 

600 

7200 

4-1/2 

3 

Anhydrous 

Ammonia 

850 

3600 

4 

2 

Water 

600 

7200 

4-1/2 

5-6 

In  view  of  these  observations  it  is  useful  to  perform  stability 
evaluations  on  a  seal  having  similar  materials  and  a  comparable  size. 
With  S  representing  the  carbon  and  H  the  tungsten  carbide,  the  base 
case  data  chosen  foTow: 


k$  =  3.5  x  10"6  m2/s 
K$  =  5.0  W/mK 
Es  =  2.0  x  1010  N/m2 
G$  =  0.83  x  1010  N/m2 
as  =  5.0  x  10-6/K 

i j  =  0.02  m 

JxS  =  3.3  x  10-9  m4 
;  .  >  6.9x  10-10  m4 

*  j 


k^  =  2.8  x  IQ’5  m2/s 
<H  =  70.0  W/mK 
Eh  =  5.0  x  1011  N/m2 
Gh  =  2.1  x  10H  N/m2 
aH  =  5.0  x  10_6/K 
=  0.02  m 

JxH  =  3.3  x  TO'9  m4 
J9$  =  6.9  x  10'10  m4 


(3-45) 
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If 


Data  conmon  to  both  seal  rings  are 

R  =  0.05  m 

u  =  0.2 

h  =  0.005  m  (3-46) 

H  =  0 

oo 

Uc  =  1.0  x  104  W/m2K 

The  value  of  contact  resistance  represents  a  case  where  a  0.38  um 
(15  yin.)  RMS  roughness  carbon  is  sliding  on  perfectly  smooth  tungsten 
carbide  with  the  contact  voids  filled  with  ammonia  gas  as  calculated 
by  formulas  developed  by  Fenech  and  Rohsenow  [15]. 

To  provide  some  frame  of  reference  for  the  model  proposed,  it  will 
be  compared  to  the  other  models  mentioned  previously.  Burton's  model 
based  on  a  conductor  sliding  on  an  insulator  [9]  will  be  termed  model 
I.  Burton's  model  which  considers  two  conducting  semiinfinite  plates 
with  Uc  *  ■»  [9]  will  be  termed  model  II.  The  present  model  will  be 
termed  model  IIIR  and  model  1 1 IS  based  on  ring  theory  and  semi  infinite 
plate  theory,  respectively. 

The  predictions  for  each  of  the  models  using  the  base  case  data 
are  shown  plotted  in  Figure  3-5.  The  results  from  model  I  would  sug¬ 
gest  that  nearly  all  seals  would  be  unstable,  a  result  which  is  con¬ 
trary  to  fact.  This  model  cannot  be  applied  to  this  particular  prob¬ 
lem.  Model  II  considers  the  effect  of  a  second  conducting  body.  This 
result  shows  that  the  second  body,  even  through  it  may  be  a  poor  con¬ 
ductor  (which  is  the  case  here  with  carbon),  greatly  increases  the 
critical  speed.  This  model  also  shows  that  the  lowest  harmonics  could 
be  expected  to  occur  first  as  speed  is  increased.  For  a  ring  the 
lowest  number  that  can  occur  is  two,  so  model  II  predicts  a  critical 
sliding  speed  of  500/s. 

Model  IIIR  shows  a  quite  different  behavior.  It  shows  that  the 
critical  speed  increases  rapidly  as  n  decreases.  This  is  caused  by 
ring  deflection.  As  the  number  n  decreases,  the  relative  s.tiffness 
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3-5.  Comparison  of  Models. 


of  the  ring  decreases  dramatically  compared  to  the  edge  stiffness  of  a 
semiinfinite  plate.  On  the  other  hand,  model  IIIR  shows  that  the  crit¬ 
ical  speed  becomes  quite  low  as  n  increases.  This  occurs  because  the 
ring  stiffness  becomes  extremely  high  at  large  n. 

Model  IIIS  is  basically  the  same  as  model  II  except  that  thermal 
contact  resistance  between  the  faces  is  considered.  The  result  shows 
that  contact  resistance  lowers  the  critical  speed.  The  shape  of  the 
curve  for  model  IIIS  is  otherwise  the  same  as  model  II. 

Now,  following  the  original  assumptions  concerning  the  validity 
of  the  ring  and  semiinfinite  plate  deflection  models,  it  can  be  argued 
that  there  should  be  a  transition  between  curves  IIIR  and  IIIS  as 
suggested  by  the  dashed  curve  in  Figure  3-5.  The  precise  nature  of 
this  transition  curve  is  not  known,  however,  it  can  be  argued  that 
if  the  thermal  and  elastic  deflection  of  both  the  ring  and  plate  mod¬ 
els  are  added  together,  a  first  approximation  for  the  entire  curve  will 
be  obtained.  This  can  be  rationalized  by  noting  that  at  low  n  the  ring 
deflection  will  dominate  and  the  effect  of  the  semiinfinite  plate  deflec¬ 
tion  will  be  small  in  comparison  and  vice  versa  at  high  n.  The  result 
is  shown  in  Figure  3-6.  Model  IIISR  represents  the  combination  of  the 
two  models. 

The  result  of  this  combination  is  a  curve  having  a  nose  which 
determines  the  lowest  critical  speed.  For  the  base  case  results,  this 
speed  is  about  60/s  at  n  =  6-7.  The  shape  of  this  curve  suggests  that 
as  speed  increases  there  will  suddently  occur  an  instability  with  n  as 
prescribed  by  the  curve.  A  further  increase  in  speed  would  cause  n  to 
either  increase  or  decrease  depending  on  which  way  the  curve  is  followed. 

For  the  experimental  data,  the  curve  suggests  that  either  n  equal 
2  or  3  or  that  n  be  very  large  (beyond  the  range  of  the  curve). 

The  effect  of  thermal  contact  resistance  is  also  shown  on  Figure 
3-6.  Contact  resistance  was  reduced  a  factor  of  100  to  simulate  the 
effect  of  having  a  liquid  between  the  faces,  and  the  curves  were  re¬ 
calculated.  For  this  case  the  first  critical  speed  moves  up  to  500/s. 

At  the  operating  speed  of  the  water  seal  described  previously,  7200  RPM, 
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Model  III 


r 


the  curve  predicts  values  of  n  of  5  or  6  or  again  very  high  (16  to  18). 

The  model  provides  for  convection  losses  through  H^.  The  effect 
of  convection  losses  is  to  reduce  the  amplitude  of  the  thermal  deflec¬ 
tion.  The  results  presented  so  far  were  for  =  0.  Introducing 
-  1000  W/m^K  into  the  model  causes  #  -*■».  This  result  suggests 
that  convection  cooling  could  be  very  important  in  reducing  or  elimin¬ 
ating  thermoelastic  instability  in  seal  rings,  however,  further  studies 
are  needed  to  understand  this  effect. 


I 


Chapter  4 

Wavy  Seal  Analysis 


Wavy  Seal  Theory 

As  will  be  shown  later  in  this  chapter,  the  modeling  of  the  exper¬ 
imental  seal  described  herein  requires  a  more  comprehensive  model  than 
that  developed  for  hydrodynamic  lubrication  in  previous  reports  and 
papers  [5,1,6,16,17].  The  primary  difference  between  the  behavior  of 
the  actual  seal  and  the  ideal  seal  assumed  in  the  previous  analysis  is 
that  in  the  actual  seal  there  is  a  tilt  superimposed  over  the  waviness. 
In  the  previous  analysis  it  was  assumed  that  the  seal  faces  alway  re¬ 
main  radially  parallel  in  spite  of  waviness  (except  for  some  of  the 
results  contained  in  Chapter  3  of  reference  [1]).  As  was  shown  by 
these  results,  tilt  of  the  faces  greatly  changes  the  hydrodynamic 
effect  in  the  seal.  Hydrostatic  effects  may  become  much  more  signi¬ 
ficant  than  hydrodynamic  effects.  Thus,  in  this  more  comprehensive 
model,  certain  features  have  been  added  to  account  for  a  completely 
arbitrary  film  shape  and  wear. 

The  development  of  this  more  comprehensive  model  to  the  level  of 
sophistication  of  the  previous  model  [1]  is  quite  difficult,  so  certain 
approximations  are  made.  First,  a  simplified  approach  to  evaluating 
the  effects  of  surface  roughness  is  taken.  Second,  because  the  cavity 
shape  becomes  more  complex,  approximate  cavity  boundary  conditions  are 
used  to  simplify  the  solution.  The  model  will  now  be  described  in 
detail . 

Reynolds  Equation  for  Rough  Surfaces 

The  Reynolds  equation  predicts  the  pressure  distribution  for  an 
arbitrary  film  shape.  Both  hydrostatic  and  hydrodynamic  effects  are 
accounted  for.  In  previous  reports  [5,1,6]  considerable  attention 
was  given  to  how  the  Reynolds  equations  should  be  modified  to  account 
for  surface  roughness.  A  contribution  to  the  literature  on  the  subject 
was  made  [7].  Certain  discrepancies  between  the  Patir-Cheng  theory  [18] 
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and  the  Teale  and  Lebeck  interpretation  of  their  theory  were  found  [7]. 
By  including  the  effects  of  side  leakage,  it  was  found  that  flow  resis¬ 
tance  decreased  compared  to  the  Patir-Cheng  results.  Recently,  Tonder 
[19]  also  computed  distinctly  different  results  than  found  by  Patir  and 
Cheng.  Even  though  they  disagree,  both  theories  show  that  the  surface 
roughness  correction  used  in  previous  wavy  seal  model  [1],  as  origi¬ 
nally  introduced  by  Christensen  [20]  and  Tonder  [21],  provides  too  littl 
flow  resistance.  Thus,  before  making  calculations  based  on  a  new  seal 
model,  some  resolution  of  these  issues  is  needed. 

To  provide  some  insight  into  this  issue,  the  flow  factors  as  com¬ 
puted  by  Patir  and  Cheng  [18]  and  Tonder  [19]  are  compared.  For  one¬ 
sided  roughness  with  the  smooth  surface  moving  (which  simulates  the 
carbon  versus  hard  face  roughness  combination  quite  well)  and  an 
h/a  =  2  or  F  =  h/c  *  0.67  the  shear  flow  factors  calculated  are  0.7  and 
0.8  (compared  to  1.0  for  smooth  surfaces)  from  the  Patir-Cheng  [18] 
and  Tonder  [19]  results,  respectively.  Pressure  flow  factors  are  0.75 
and  1.2  (compared  to  1.0  for  smooth  surfaces),  respectively.  While  the 
shear  flow  factors  agree  to  some  extent,  the  pressure  flow  factors  show 
opposite  trends. 

Given  this  irreconci liable  difference  in  theory  in  conjunction 
with  the  necessity  to  produce  a  useful  though  approximate  seal  model, 
it  was  decided  that  it  would  be  best  to  avoid  choosing  either  of  these 
models  until  the  basic  issues  are  resolved.  It  was  decided  to  com¬ 
pletely  neglect  the  effects  of  surface  roughness  on  the  Reynolds  equa¬ 
tion.  This  choice  has  the  advantage  that  when  the  above  issue  is 
resolved,  then  the  results  can  be  recomputed  based  upon  correction  due 
to  roughness  as  opposed  to  being  recomputed  because  the  roughness  model 
chosen  was  incorrect.  The  former  is  much  preferred. 

One  can  in  fact  neglect  surface  roughness  effects  on  seal  lubrica¬ 
tion  while  introducing  only  a  small  error.  Experience  shows  that  in  a 
wavy  or  tapered  seal,  the  film  thickness  is  as  low  as  F  =  0.67  for  only 
a  small  fraction  of  the  face  area.  For  most  of  the  seal,  F>  1,  where 
no  significant  surface  roughness  effects  occur.  Thus,  any  error  intro¬ 
duced  occurs  in  only  a  small  fraction  of  the  seal  face.  The  error 


introduced  is  not  large  to  begin  with  since  at  h  =  0.67,  typically  a 
lower  limit  for  a  seal,  =  0.75  or  1.25,  a  25  percent  error. 

Seal  Model 

Even  though  the  effect  of  roughness  on  the  Reynolds  equation  is 
to  be  neglected,  computation  of  mechanical  load  support  and  friction 
require  the  use  of  surface  roughness  concept.  Thus,  surface  roughness 
is  shown  in  Figure  4-1,  along  with  seal  geometry,  h  is  the  nominal 
film  thickness  at  a  point  r,e  and  in  general  in  a  function  of  r  and  e. 
h  is  defined  to  the  mean  level  of  the  rough  surface.  It  is  assumed 
that  the  roughness  is  distributed  over  a  range  -c  to  +c  and  that  the 
distribution  is  polynomial. 

f(h)  3  (c2  -  h2)3  .  (4-1) 

32c' 


Tctal  film  thickness  is  given  by 

H  3  h  +  6  ,  (4-2) 

where  5  represents  the  random  part  of  the  film  thickness. 

Based  upon  the  assumption  that  roughness  effects  are  neglected, 
the  Reynolds  equation  is 


Flow  is  given  by 


(4-4) 

(4-5) 


where  h  represents  the  nominal  film  thickness  as  defined  above..  Sym¬ 
bols  are  nondimensional ized  as  shown  in  the  list  of  symbols.  Because 
of  the  definition  of  h,  it  is  possible  for  h  to  become  negative.  When 
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Section  AA 


Figure  4-1.  Seal  face  Geometry  and  Surface 
Roughness. 
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this  occurs,  h  is  set  equal  to  zero. 
Boundary  conditions  are 


p(ri )  =  Pi  > 

(4-6) 

P(ro}  =  po  • 

(4-7) 

In  the  previous  model  [1],  cavitation  was  accounted  for  by  an  exact 
procedure  which  provided  flow  continuity  across  the  cavity.  A  unique 
numerical  technique  was  developed  for  the  solution.  For  the  present  case 
the  cavity  boundary  condition  is  treated  by  first  solving  for  the  com¬ 
plete  pressure  distribution  disregarding  cavitation  and  negative  pres¬ 
sures.  Then,  all  pressures  less  than  the  cavity  pressure  are  set 
equal  to  the  cavity  pressure.  This  procedure,  equivalent  to  the  half 
Sommerfeld  condition  in  journal  bearings,  introduces  some  error  into 
the  load  support  and  leakage  calculations.  However,  because  of  the 
complications  of  multiple  and  arbitrarily  shaped  cavities  introduced 
in  this  more  general  problem  and  because  of  the  greater  numerical 
effort  required  to  obtain  deflection  solutions,  the  simplification 
was  found  necessary  at  least  for  the  computation  of  the  preliminary 
results  presented  here. 

Equation  (4-3)  is  solved  numerically  by  over  relaxation  methods 
as  before  [1],  Care  must  be  taken  in  setting  up  the  equation  so  that 
film  thickness  slope  discontinuities  do  not  introduce  large  errors. 

Load  support  is  calculated  as  before  [1]  and  consists  of  the  fluid 
pressure  component  plus  the  mechanical  component 


+  pbf]  r  dr  de  , 


(4-8) 


where 


bm  =  P(H  >  h)  , 


(4-9) 
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the  probability  that  contact  occurs  or  the  fraction  of  asperity  area 
in  contact.  For  the  polynomial  distribution 


-  [16  +  35fi  -  35h3  +  21  h5 


(4-10) 


b 


f 


(4-11) 


which  is  the  fraction  of  area  subjected  to  fluid  pressure.  For  practi¬ 
cal  purposes, 

bf  =  1.0  .  (4-12) 

Pm  is  the  compressive  strength  of  the  asperities  as  before  [1]. 

The  applied  load  is  given  by 

W*  =  ir(l  -  ??)  [p0B  +  p.(l  -  B)  +  psp]  (4-13) 

and  for  equilibrium,  the  film  thickness  must  be  adjusted  so  that 

W*  =  W  .  (4-14) 

Leakage  is  given  by 


q  .  f2*.S LIE 
Qf  Jo  12  3r 


r.  de  , 


(4-15) 


or 


y»  2tt  r3  - 

0  12  »r 


r  de  . 
o 


(4-16) 


When  cavitation  occurs  and  negative  pressures  are  set  equal  to  the 
cavitation  pressure  as  described,  some  discrepany  between  Q..  and  0Q 
is  expected. 


Friction  is  calculated  as  before  [1].  Shear  stress  for  the  fluid 


is  given  by 


t  =  nru)E(i-)  +  JF  |f  E(h)  • 


Assuming  a  friction  radius  r^  where 


r?> 


(4-17) 


(4-18) 


Then  the  total  fluid  friction  force  due  to  the  fluid  is 


dr  de  . 


(4-19) 


E(jj-)  was  given  previously  and  represents  the  expected  value  of  one  over 
the  fluid  film  thickness.  E(H),  the  expected  value  of  the  film  thick¬ 
ness,  is  taken  as  hi n  accordance  with  the  present  approximation.  In 
the  cavity  region  =  0,  and  the  first  term  in  the  above  integral  is 

0  U 

multiplied  by  the  ratio  of  the  film  thickness  at  the  beginning  of  the 
cavity  to  the  current  film  thickness  to  account  for  the  fraction  of 
the  full  film. 


The  factor 


(4-20) 


is  used  as  before  [1]  to  limit  the  viscous  shear  stress  to  the  mechan¬ 
ical  shear  strength  and  is  used  in  the  integration  for  E[— ). 

'fi; 

Friction  due  to  mechanical  contact  is  given  by 


T 


P  b 
ys  m 


(4-21) 


p.b  r 
s  m 


dr  de  , 


(4-22) 
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where  p$  is  the  mechanical  shear  strength  and  bm  again  is  the  fraction 
of  area  in  mechanical  contact. 

Total  friction  is  given  by  the  sum  of  the  above  two  terms,  and 
coefficient  of  friction  is  thus. 


F,  +  F  F,  +  F 
f _ m  _  f _ m 


(4-23) 


As  shown  in  previous  work  [1]  the  average  wear  rate  of  the  soft 
face  is  proportional  to 

wear  as  (1  -  F)  ,  (4-24) 

where  F  represents  the  fraction  of  the  load  carried  by  fluid  pressure. 

F  is  computed  using  equation  (4-8)  with  the  mechanical  load  support 
term  omitted. 


Film  Thickness  Shape 

As  will  be  shown  later,  the  film  thickness  function  consists  of 
some  portion  due  to  wear  and  some  portion  due  to  net  deflection. 
Deflection  creates  changes  in  film  thickness  by  displacement  of  the 
centroid  of  the  ring  plus  rotation  of  the  ring  about  its  centroidal 
axis.  It  is  assumed  that  deflection  occurs  only  in  the  soft  (low 
modulus)  face.  Figure  4-2  shows  these  two  components  of  deflection. 
Coordinate  system  x-y  is  fixed  in  space  a  constant  distance  C  from 
the  nondeflecting  ring.  The  ring  centroid  deflects  in  the  y  direc¬ 
tion  and  in  the  <j>  direction  about  the  centroidal  axis.  For  small 
angles  $,  one  can  show  that 

h  =  const  -  v  +  (r  -  r  )  $  .  (4-25) 


In  general,  the  ring  deflections  v  and  $  can  be  represented  as  a 
Fourier  series.  The  v  term  can  be  written  in  terms  of  an  equivalent 
h  term.  Thus, 
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h  =  h  +  w(r)  +  l  (h  cos  ne  +  h.  sin  ne) 

0  n=l  an  Dn 

oo 

+  (r  -  r  )  l  (<|>  cos  ne  +  <t>.  sin  ne)  .  (4-26) 

c  n=1  an  Dn 

The  term  hQ  is  the  minimum  value  of  h.*  In  general  load  support  W  in¬ 
creases  with  decreasing  h  .  Therefore,  hQ  is  determined  by  the  con¬ 
dition 

W(hQ)  *  W*  ,  (4-27) 

where  W(hQ)  and  W*  are  given  by  equations  (4-8)  and  (4-13),  respectively. 

w(r)  represents  the  wear  term.  For  a  seal  with  a  moving  waviness 
or  an  axisynmetric  seal,  wear  becomes  a  function  of  r  only.  For  the 
moving  waviness  seal  the  wear  simply  becomes  uniformly  distributed  in 
the  0  direction. 

Assuming  h  ,  h.  ,  $  and  <j>.  are  known,  w(r)  can  be  easily  found 
an  °n  an  °n 
for  the  fully  worn  in  condition  by 


w(r)  =  -min 


l  (h  cos  ne  +  h.  sin  ne) 
n=l  an  Dn 


+  (r  -  r  )  l  U  cos  ne  +  <f>b  sin  ne) 
n=l  n  n 


(4-28) 


This  condition  requires  that  at  all  values  of  r,  the  seal  will  be  just 
touching  at  some  angle  e.  This  condition  is  necessary  for  uniform  wear 
of  the  seal  ring  which  is  the  condition  reached  after  the  wave  has 
moved  around  the  seal  a  number  of  times.  Partial  wear  conditions 
(where  only  a  fraction  of  the  width  of  the  face  has  worn)  can  also  be 
determined  as  shown  later. 

The  deflection  terms  h  ,  h.  ,  a  and  arise  from  all  loads 

an  Dn  an  Dn 

imposed  on  the  ring  r.  These  loads  include  deliberately  applied  loads 
such  as  waviness  pressure  and  face  loads  generated  by  the  face  pressure 

*w(r)  is  adjusted  to  make  this  true. 
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distribution.  Before  considering  complications  of  deflection,  the 

solution  to  the  previous  seal  model  equations  will  be  given  for  an 

arbitrary  set  of  h  ,  h.  ,  ,  <f>.  to  see  how  seal  lubrication  and 

“n  Dn  an  Dn 

leakage  behave  for  this  general  type  of  deflection.  For  this  purpose, 

it  will  be  assumed  that  all  significant  disturbances  h  ,  h,  ,  $  ,  <j>. 

an  n  an  Dn 

are  of  the  third  harmonic.  The  three  will  be  dropped  to  shorten  the 
subscript.  Thus, 

h  =  hQ  +  w(r)  +  (ha  +  (r  -  rc)  4>a)  cos  39 

+  (hb  +  (r  -  rc)  4>b)  sin  3e  .  (4-29) 


Solution 

For  a  given  hfl,  hb,  <j>a,  <j>b ,  material  properties,  roughness,  sealed 
pressure,  balance  ratio,  etc.  the  general  method  of  solution  is  as 
follows: 

1)  Calculate  w(r)  using  equation  (4-28) 

2)  Assume  a  value  for  hQ 

3)  Solve  Reynolds  equation  (4-3)  fo.'  the  given  boundary  condi¬ 
tions  using  over  relaxation 

4)  Set  all  negative  pressures  equal  to  zero 

5)  Calculate  load  support  using  equation  (4-8) 

6)  Compare  W  to  W*  and  adjust  hQ  using  a  secant  root  finding 
technique.  Return  to  step  3)  and  repeat  until  W  =  W* 

7)  Calculate  friction  and  leakage 

To  demonstrate  the  solution,  a  base  case  common  to  the  test  con¬ 
ditions  of  Chapter  2  has  been  established  as  follows: 

pi  =  0 

pQ  =  3.45  MPa  (500  psi) 
pm  =  262  MPa  (38000  psi) 

Ps  *  0.1  pm 

Psp  *  0.207  MPa  (30  psi)  (continued  on  following  page) 
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c  *  0.51  pm  (0.14  pm  CLA)  (20  pin.)  (5.4  pin.  CLA) 

B  *  1.0 

n  =  6.83-10"4  Pa-s,H20  at  38°C  (100°F)  (0.99-10"7  Ib-s/in.2) 

RPM  =  1800  (188.5/s) 
r..  s  48.26  mm  (1.900  in.) 
rQ  =  53.04  m  (2.088  in.) 
r  =  52.83  nm  (2.080  in.)  (4-30 

v 


Corresponding  nondimensional  constants  are: 


p0  =  2.46*10 

pm  =  1.87-10"1 

ps  =  1.87-10'2 

Pcn  =  1.48-10'4 

rQ/c  -  1 . 044 • 1  O' 

r.  =  0.9100 

r.  ■  0.9962 
c 


(4-31 


The  waviness  selected  for  this  example  is 


R  =  5.0 

a 


where 


Rb  -  0.0 


(4-32 


♦  =  42.43  ib  *  42.43 


<P  s  <t> 


(4-33 


Using  the  solution  method  outlined  with  a  20  x  20  grid,  the  pres¬ 
sure  distribution  is  shown  in  Figure  4-3.  Note  the  presence  of  the 
odd-shaped  cavity  (denoted  by  zero  pressure)  and  the  nonsymmetric  pres 
sure  distribution.  This  cavity  shape  is  much  different  than  that 
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Figure  4-3.  Pressure  Distribution  for  Base  Case 


obtained  for  radially  parallel  faces  [1].  The  arbitrary  nature  of  the 
cavity  shapes  produced  by  combined  waviness  and  tilt  necessitates  the 
use  of  the  approximate  cavity  boundary  conditions  discussed  previously. 

Based  on  equation  (4-28)  the  wear  was  calculated.  The  film  thick¬ 
ness  is  shown  plotted  in  Figure  4-4  at  four  different  angular  positions 
At  e  =  0,  the  film  shape  is  convergent  but  maintains  a  gap  at  r^ .  At 
30°  the  film  shape  is  convergent  but  touches  at  r^ .  At  60°  the  film 
shape  is  divergent  and  touches  at  rQ;  whereas  at  90°  the  film  shape  is 
divergent  and  does  not  touch  at  r  .  The  face  is  worn  into  a  curved 
shape.  From  the  results  shown  in  Figures  4-3  and  4-4,  it  is  apparent 
that  combined  waviness  and  tilt  may  in  general  lead  to  nonsymmetric 
pressure  distributions,  curved  wear  surfaces,  and  film  shapes  which 
vary  dramatically  with  e. 

The  computed  results  for  the  base  case  are 

n  =  0.026 

Q  =  37.0  ml /min 

Tq  *  7.3  N-m  (4-34) 

%  =  77.9 

hQ  =  0.577 

While  torque  is  about  one-half  of  that  of  a  comparable  flat  faced  seal, 
leakage  is  quite  large.  Leakage  is  large  because  of  the  gaps  shown  in 
Figure  4-4.  For  a  flat  faced  seal  %  =  49.  Thus,  average  wear  rate 
for  the  above  seal  would  be  expected  to  be  (1  -  0 . 78 )/ ( 1  -  0.49)  »  0.43 
or  43%  of  that  of  the  flat  faced  seal  according  to  equation  (4-24). 
Hydrodynamic  action  is  clearly  present  as  is  indicated  by  the  cavita¬ 
tion  region  shown  in  Figure  4-3.  Hydrodynamic  action  in  the  base  case 
is  also  clearly  shown  by  comparing  to  a  solution  for  the  above  base 
case  when  the  speed  was  set  to  0.18  RPM.  For  0.18  RPM,  the  results 


u  =  0.047 
Q  =  31.2  ml /min 

(4-35) 

Tq  =  13.3  N  m 
%  =  52.8 

The  friction  torque  is  much  larger  and  the  %  is  much  lower  indicating 
that  the  fraction  of  load  supported  by  fluid  pressure  in  this  static 
case  is  much  smaller  than  in  the  dynamic  case. 

General  Studies 

To  evaluate  the  potential  wear  and  friction  reduction  using  com¬ 
bined  parallel  face  waviness  and  tilt,  several  parameter  studies  were 
made.  First,  Figure  4-5  shows  leakage  and  torque  as  functions  of  J 

a 

for  a  fixed  ha.  Note  that  both  Fi^  and  are  zero  so  that  the  wave 
and  the  tilt  are  in  phase.  Considering  first  the  0  RPM  or  static 
case,  it  is  clear  that  a  certain  tilt  combined  with  the  given  waviness 
produces  a  minimum  value  of  torque.  Leakage  values  are  quite  large 
until  this  critical  value  is  reached,  after  which  they  remain  quite 
small.  At  1800  RPM  the  general  behavior  is  similar  except  that  torque 
values  are  much  lower  than  for  the  static  case,  and  leakage  values  are 
slightly  higher.  A  comparison  of  the  0  and  1800  RPM  friction  curves 
clearly  indicates  that  hydrodynamic  effects  are  occurring  and  increas¬ 
ing  the  fraction  of  fluid  pressure  load  support. 

The  curve  also  shows  for  the  1800  RPM  case  that  operation  at  what 
will  be  called  the  optimum  tilt  (the  minimum  point  on  the  torque  curve) 
provides  much  better  seal  performance  than  for  parallel  face  waviness 
where  =  0.  Leakage  is  much  lower  and  torque  is  somewhat  lower. 

Q 

This  result  points  out  a  method  whereby  hydrodynamic  effects  may  be 
used  to  great  advantage  while  leakage  remains  controlled  at  low  levels. 

Considering  the  optimum  point,  a  detailed  study  shows  that  this 
occurs  when 
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1 


*,  =  - - —  (4-36) 

rc  -  ri 

which  is  the  case  where  double  touching  just  begins.  An  illustration 
of  wear  profiles  for  some  of  the  cases  shown  in  Figure  4-5  provides  an 
explanation  for  the  optimum  point  as  well  as  other  behavior.  For 
$a  =  -60,  Figure  4-6  shows  that  the  faces  are  either  parallel  or  diver¬ 
gent.  This  explains  the  high  torque  and  leakage  for  both  the  static 
and  dynamic  cases  at  negative  When  5  =  60  (very  close  to  the 

a  a 

optimum  point),  the  faces  are  either  parallel  or  convergent.  This  con¬ 
vergence  explains  why  friction  is  so  low  in  both  the  static  and  dynamic 
cases.  The  seal  faces  touch  at  r.  all  around  the  seal.  This  explains 
the  low  leakage,  i.e.,  there  is  no  signficant  gap  at  any  point  around 
the  seal.  As  becomes  larger,  leakage  remains  low  because  touching 

a 

occurs  at  some  point  everywhere  around  the  seal.  However,  the  gap 
becomes  divergent  for  a  fraction  of  the  seal  face  and  considerable  load 
support  is  lost  compared  to  the  <t>.  =  60  case. 

CL 

From  this  illustration  the  principles  for  operation  of  a  static 
or  dynamic  seal  utilizing  combined  waviness  and  tilt  where  friction 
is  to  be  at  a  minimum  and  leakage  must  be  small  becomes  clear.  For 
minimum  leakage,  touching  must  occur  all  around  the  seal.  This  does 
not  require,  as  shown  by  the  middle  case  of  Figure  4-6,  that  a  large 
radial  fraction  must  contact;  a  small  fraction  as  shown  will  suffice. 

It  is  only  necessary  that  a  touching  sealing  dam  be  established  every¬ 
where.  Concerning  friction,  where  the  face  departs  from  the  parallel 
condition  it  must  become  convergent  rather  than  divergent.  The  con¬ 
vergent  face  causes  the  fluid  pressure  to  decrease  at  a  lower  rate 
radially  across  the  face  than  for  a  divergent  face,  thus  giving  much 
greater  load  support.  In  fact,  load  support  is  better  than  for  par¬ 
allel  faces  even  in  the  hydrodynamic  case. 

It  also  becomes  clear  from  this  discussion  that  there  is  much 
merit  in  simply  designing  a  seal  face  with  a  converging  radial  taper 
uniform  around  the  seal,  i.e.,  no  waviness.  Indeed,  as  shown  by  Snapp 
and  Sasdelli  [22]  this  approach  does  lead  to  low  friction.  However, 
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there  is  no  mechanism  by  which  to  maintain  the  taper  on  an  indefinite 
basis  as  in  the  case  for  a  wavy  seal.  Therefore,  to  minimize  wear 
and  degradation  of  operation,  the  convergent  taper  seal  must  be  de¬ 
signed  having  a  large  gap  and  consequently  large  leakage.  The  wavy 
approach  eliminates  this  disadvantage  because  the  gap  can  be  essen¬ 
tially  closed  and  leakage  values  made  quite  low. 

In  Figure  4-5,  the  waviness  and  tilt  were  in  phase.  Now  the 
effect  of  a  phase  shift  is  considered.  For  the  results  shown  in 
Figure  4-7  the  following  conditions  were  maintained. 

h  =  5.0 

a 

hb  =  o 

*  =  60  (4-37) 

$  COS  a 

$  sin  a 

This  has  the  effect  of  moving  the  tilt  an  angle  a  relative  to  the  waviness 
h  .  The  results  in  Figure  4-7  clearly  show  that  leakage  is  at  a  minimum 
only  when  the  waviness  and  tilt  are  in  phase  for  both  the  static  and 
dynamic  cases.  Torque  for  the  static  case  also  reaches  a  minimum  under 
these  same  conditions.  Torque  for  the  dynamic  case  does  reach  a  rela¬ 
tive  minimum  at  a  =  ir,  however,  the  case  where  a  =  0  clearly  gives  the 
lowest  torque.  In  conclusion,  based  on  this  study,  waviness  and  tilt 
should  be  in  phase  to  obtain  the  lowest  leakage  and  torque. 

Optimum  Conditions 

In  Figure  4-8,  the  results  of  a  study  of  performance  for  the  base 
case  seal  under  the  optimum  condition  of  equation  (4-36)  is  shown, 
h  *  0  represents  a  parallel  face  case.  As  waviness  (and  therefore 
tilt)  are  increased,  torque  decreases  and  leakage  increases.  In  the 
dynamic  case,  low  values  of  torque  are  reached  at  relatively  low 
leakage  values.  The  selection  of  an  operating  point  is  clearly  a 
trade-off  question.  If  relatively  large  amounts  of  leakage  can  be 
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Figure  4-7.  Torque  and  Leakage  as  Functions  of  Phase  Angle. 
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Figure  4-8.  Leakage  and  Torque  at  Optimum  Conditions. 
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tolerated,  then  R  should  be  large.  If  only  a  small  leakage  can  be 

a 

tolerated,  then  some  compromise  on  torque  must  be  made. 

Figure  4-8  shows  an  additional  advantage  to  this  design  concept 
not  pointed  out  previously.  The  static  torque  is  significantly  reduced 
by  the  same  changes  that  reduce  the  dynamic  torque.  Thus,  even  at 
lower  operating  speeds,  the  seal  will  have  lower  friction  than  a  parallel 
seal.  This  feature  would  minimize  wear  at  low  speeds  and  during  start 
up  conditions.  Such  a  seal  design  is  more  versatile  than  a  seal  that 
derives  its  lift  purely  hydrodynamically. 

Ring  Deflection 

The  combined  waviness  and  tilt  model  just  described  is  applicable 
to  the  experimental  seal.  However,  before  the  model  can  be  used,  the 
waviness  in  the  experimental  seal  must  first  be  determined.  The 
waviness  in  the  seal  occurs  due  to  various  deflections.  The  forces 
producing  these  deflections  are  the  waviness  producing  mechanism, 
the  system  of  pressures  acting  on  the  face,  and  the  drive  forces.  The 
resulting  deflections  are  complex,  and  a  thorough  analysis  must  be 
made. 

Loads 

As  a  general  approach  to  calculating  deflection,  it  is  first 
assumed  that  all  loads  are  distributed  such  that  the  load  can  be  re¬ 
presented  by  a  Fourier  series.  This  has  the  advantage  of  greatly  sim¬ 
plifying  the  calculations  for  the  simultaneous  solutions  for  deflection 
needed  as  the  final  result,  because  certain  harmonics  or  terms  in  the 
Fourier  series  become  negligible. 

A  segment  of  the  carbon  ring  is  shown  in  Figure  4-9.  It  is  assumed 
that  all  significant  deflections  occur  in  the  carbon  ring  compared  to 
the  W-C  because  of  its  low  Youngs  modulus.  The  loads  shown  are  assumed 
to  be  distributed  and  functions  of  0.  The  distributed  force  loads  are 
assumed  to  be  acting  through  the  centroid.  Given  the  waviness  cylinder 
arrangement  shown  in  Chapter  2,  which  applies  pressure  to  the  inside 
surface  of  the  seal  ring,  and  the  pressure  and  friction  on  the  seal  face. 
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the  sources  of  the  various  loads  shown  are  as  follows. 


Px  -  caused  by  waviness  pressure 
Py  -  caused  by  face  pressure 

py  -  caused  by  friction  between  the  0-rings  of  the  waviness 

cylinder  and  the  inside  surface  of  the  ring  and  by  friction 
at  the  contact  face 

mx  -  caused  by  friction  of  the  0-rings  and  friction  at  the  con¬ 
tact  face 

m  -  same  as  m 

y  * 

m.  -  caused  by  the  wav i ness  pressure 

It  can  be  shown  that  deflection  due  to  friction  forces  is  small 

relative  to  deflection  due  to  the  other  loads.  This  is  true  because 

friction  forces  on  the  face  and  between  the  ring  and  waviness  cylinder 

are  relatively  small.  Therefore,  one  need  only  consider  p  ,  p  and 

x  y 

v 

Now,  considering  the  waviness  pressure.  Figure  4-10  shows  px  and 
m.  caused  by  the  six  gas  pressure  regions  acting  on  the  carbon  ring. 

It  is  assumed  that  the  force  distribution  can  be  approximated  as  a 
uniform  force  acting  over  a  distance  d,  the  diameter  of  the  hole, 
and  having  a  total  force  equal  to  that  of  the  gas  pressure  times  the 
area  of  the  hole.  The  factor  r.j/rc  puts  the  force  and  moment  in  terms 
of  per  unit  length  of  centroidal  circumference.  Moment  arm  e  is  the 
distance  in  the  y  direction  between  the  center  of  the  hole  and  the 
centroid. 

A  Fourier  analysis  for  the  pressure  and  moment  shown  in  Figure 
3-10  gives 


Pv  =  Pn  \  sin  cos  6e  +  j  sin  ^  cos  12e  +  •••! 
x  9  rc  [2  r.  4  r.  j 


=  pva  cos  6e  +  pv,  cos  12e  +  •••  , 
xa6  xal 2 


(4-38) 
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mun  -  cos  3e  +  m  cos  6e  +  ••• 

ep  °Pa3  °Pa6 


+  m9pb3  sin  30  +  mepb6  cos  69  + 


(4-43) 


Only  harmonic  multiples  of  three  appear  because  there  are  three  waves. 
Both  sin  and  cos  terms  appear  because  the  wave  is  in  general  nonsymmetric 
about  the  origin. 

Loading  due  to  both  the  waviness  pressure  and  face  load  appear  in 
series  form.  Deflection  caused  by  these  loadings  is  now  considered. 
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Deflection 

Based  upon  equation  (2-30)  of  reference  [23],  one  may  solve  for 

the  deflection  of  a  ring  caused  by  general  harmonic  loads  p  ,  p  and 

y  * 

mQ.  For  the  nth  harmonic  of  deflection  the  result  is 
0 


P  r 
Kxa„  c 
n 

EJ.. 


1)‘ 


"ea^c 

EJ„ 


(n2  -  l7 


p  r 
Pyan  c 
EJ.. 


_A+_n__ 
n2(n2  -  1)‘ 


(4-44) 


xan  c  Jxv  1 

ej  J  71  77 

x  y  (n  -  1 ) 
3 

p  r 

yan  c  A  +  1 


Vc  An2,! 


EJ. 


(n^ 


EJx  (7o7‘ 


(4-45) 


<j>.  and  v.  are  given  by  the  same  expression  as  above  but  substituting 
bn  bn 

b's  for  a's.  The  symbol  A  above  is  given  by 


A  = 


GJ9  * 


(4-46) 


Jx  Jy  Jxy  Je  are  standard  section  properties  for  a  ring  cross  section 
and  are  given  in  the  list  of  symbols.  The  net  deflection  then  becomes 


v  *  l  v  cos  ne  +  l  v.  sin  ne  , 
n  an  n  n 

p  *  l  p  cos  ne  +  l  pb  sin  ne  . 
n  an  n  Dn 


(4-47) 

(4-48) 


Based  on  the  actual  seal  ring  of  interest,  the  ratio  Jxy/Jy  is 
quite  small.  Therefore,  the  effect  of  px  is  insignificant  and  this 
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term  will  be  dropped  for  the  balance  of  this  work.  Also  note  that  all 
deflection  terms  contain  a  factor  (n  -  1)  in  the  denominator.  This 
means  that  deflection  of  the  6th  harmonic  will  be  5  percent  of  the 
magnitude  of  deflection  of  the  3rd  harmonic  assuming  comparable  loading 
terms.  Therefore,  terms  of  the  6th  harmonic  and  higher  will  be  neglected 
for  the  remainder  of  this  work.  Since  all  remaining  terms  will  be  for 
n  =  3,  the  n  will  be  dropped  from  the  symbols. 

Now  assuming  that  the  reference  point  (e  =  0)  for  the  system  is  at 
the  center  of  one  of  the  waviness  gas  pressure  circles,  and  with  the 
help  of  equations  (4-25)  and  (4-26)  one  can  combine  equations  (4-44) 
through  (4-48)  to  obtain  the  following. 


h  =  hQ  +  w(r)  +  (ha  +  (r  -  rQ)  $a)  cos  3e 
+  (hfa  +  (r  -  rc)  $b)  sin  3e  , 


where 


ha  “  hma  +  hga  +  hpa  * 
hb  '  "mb  +  hgb  +  hpb  • 


where 


m 


ma  epa  EJ 


1  +  A 
x  (n2  -  1) 


h  -  m  c  1  +  A 
hga  -  %a  EJ,  ^ ~ ’ 


rc  A  +  n2 


hpa  *  -pya  EJV  „2,„2 


\  n  (n  -  1)' 


(4-49) 

(4-50) 

(4-51) 

(4-52) 

(4-53) 

(4-54) 


(4-57) 


Pb 

and  where 


hnw  =  -P 


rc  A  +  n2 


*b  EJx  n2(n2-l)2 


b  =*  +  A  +  A  , 

ra  *ma  *ga  vpa  ’ 


^*b  ”  ^mb  +  ^gb  +  ^pb  * 


where 


.  _  rc  An2  +  1 

nB  9pa  EJx  („2  -  I)2 

r  2  2 

_ _  c  An^  +  1 

*ga  '  mega 

_  3 


’pa  ‘  ‘pya  EJX  ^7  .  ,,a 

.  .  _  rc  An2  +  1 

♦mb  -  »epb  EJX  ^777  • 


A  +  1 


r  2 
c  An4  +  1 


♦gb  *  magb  E0X  ^777  ■ 
_  3 


pb  =  ’Pyb  EJ  77 


c  A  +  1 


x  (nc  -  I)1 


(4-58) 

(4-59) 

(4-60) 

(4-61 ) 

(4-62) 

(4-63) 

(4-64) 

(4-65) 


The  term  m.„  is  defined  by  equation  (4-39);  p„  .  p  .  ,  and  m„  , 
oya  ya  yo  opd 

m0pb  are  derived  by  making  a  Fourier  analysis  of  the  face  pressure  as  de¬ 
fined  by  equations  (4-40)  and  (4-41)  and  (4-42)  and  (4-43).  m0gb  is  zero 
if  the  9  coordinate  begins  at  the  center  of  one  of  the  gas  pressure  holes. 
The  above  terms  may  be  thought  of  as  follows. 
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hma  waviness  due  to  a  moment  of  the  face  pressure 

hga  waviness  due  to  a  moment  of  the  gas  pressure 

hpa  waviness  due  to  the  force  of  the  face  pressure 
4>ma  tilt  due  to  a  moment  of  the  face  pressure 

if>ga  tilt  due  to  a  moment  of  the  gas  pressure 

<j>pa  tilt  due  to  the  force  of  the  gas  pressure 

The  b  subscript  terms  simply  represent  the  sin  component  of  the  wave. 
Rewriting  equation  (4-49)  in  the  above  terms, 

h  ■  ho  ♦  «(<•)  +  +  "ga  *  "pa  ♦  (<~  -  rc)(*m  *  *ga  +  ♦„,)]  cos  38 

+  [hmb  *  hgb  +  "pb  *  (r  -  rc)(*n.b  +  *gb  +  *pb)]  sin  36  •  (4'66) 
Nonparallel  Faces 

From  equations  (4-53)  and  (4-54)  and  (4-61)  and  (4-62),  deflection 
ratios  are 


A  +  n2 
n2(l  +  A) 


(4-67) 


(0.12)  . 


(4-70) 


Now  for  a  given  ratio  of  the  face  pressure  harmonic  component  p^a  to 
the  gas  pressure  moment  harmonic  component  m0ga>  the  face  pressure 
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caused  waviness  hpa  is  larger  than  the  face  pressure  caused  tilt.  This 

means  that  when  the  net  waviness  component  (h  +  h  )  is  flattened  to 

'pa  ga 

nearly  zero  by  the  face  pressure,  the  net  tilt  component  will  not  be. 

In  the  original  concept  of  the  experimental  apparatus,  it  was 
assumed  that  both  tilt  and  waviness  would  be  reduced  a  proportional 
amount  by  the  face  pressure,  thus  resulting  in  relatively  small  net 
tilt  as  well  as  relatively  small  net  waviness.  The  above  relationship 
shows  that  this  is  not  true.  Even  though  waviness  flattens  to  zero, 
the  tilt  does  not  go  to  zero.  Thus  in  the  experimental  seal,  the  film 
thickness  function  is  quite  different  than  first  anticipated  [1  ]. 
Analysis  of  the  experimental  face  seal  based  upon  wavy  but  radially 
parallel  faces  led  to  predictions  that  totally  disagreed  with  the 
experiment.  After  all  above  effects  were  included,  agreement  became 
much  better  as  will  be  shown. 

After  a  thorough  analysis  of  the  deflection  behavior  it  becomes 
clear  that  it  is  very  difficult  to  deflect  a  seal  face  so  that  the 
waviness  remains  radially  parallel.  Tilt  is  almost  inevitably  intro¬ 
duced.  However,  as  the  preceeding  results  show,  a  combination  of 
waviness  and  tilt  may  actually  lead  to  a  more  improved  seal  concent 
than  parallel  face  waviness  alone. 


Solution  for  Test  Seal 


To  solve  the  combined  waviness  and  tilted  seal  problem  for  a  real 

case  becomes  much  more  difficult  than  for  the  hypothetical  cases  studied 

previously.  The  difficulty  that  arises  is  that  h,,  h  h  ,  h  .  <*  a 

ma  pa  mb  pb  ma  pa 

^mb  ^pb  are  un*(nown  and  depend  on  the  solution  itself.  Thus,  one 
must  assume  a  set  of  values,  evaluate  the  pressure,  and  then  check  the 
resulting  deflections  against  the  assumed  deflections. 

To  expedite  the  solution  the  following  approach  was  taken.  First, 
since  the  position  of  the  gas  pressure  wave  can  be  considered  to  be 
arbitrary. 


mepb  =  ’megb  ’ 


(4-71) 
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(4-72) 


m. 


epa 


+  m 


ega 


ta 


Now  1S  generally  much  smaller  than  m„  ,  so  if  one  assumes  an 

a pa  yga 

m9ga  ^ased  on  9as  pressure,  only  a  small  change  in  effective  pressure 
results.  Conditions  (4-71)  and  (4-72)  change  equation  (4-66)  to  the 
following. 


h  *  h0  ♦  »(<■)  ♦  +  hpa  ♦  (r  -  rc)(m9aK2  +  *pJ)]  cos  3e 

+  [hpb  +  (r  -  rc)  4>pb]  sin  3e  ,  (4-73) 

where  K-j  and  K«  are  given  in  equations  (4-52)  and  (4-60).  Now  since 
<t>pa  is  proportional  to  hpa  and  4>pb  is  proportional  to  hpb>  there  are 
only  three  arbitrary  constants  remaining  in  equation  (4-73).  These 
are  maa,  hpa  and  hpb.  Given  values  for  these  parameters,  the  solution 
to  the  problem  then  becomes  the  same  as  before.  Thus,  a  general  solu¬ 
tion  for  an  actual  seal  can  be  accomplished  as  follows. 

1)  Assume  a  value  for  mQa  based  on  the  gas  pressure  at  which  the 
problem  is  to  be  solved. 

2)  Guess  values  for  h„,  and  h„K  calculate  <f>  and  <p  .  from  these. 

pa  pD  Tpa  Tpb 

3)  Solve  for  the  pressure  distribution  and  equilibrium  hQ  based 
on  load  as  before. 

4)  Evaluate  integrals  (4-40)  and  (4-41)  as  functions  of  e  and 
make  a  Fourier  analysis  to  obtain  pya  p^b  mQpa  m0pb. 

5)  Find  hpa  and  hpfa  using  equations  (4-54)  and  (4-57). 

6)  Compare  the  calculated  values  of  hpa  and  hpb  to  the  assumed 

values.  If  they  agree,  then  calculated  deflection  is  equal 

to  assumed  deflection  and  the  solution  is  correct.  If 

they  do  not  agree,  then  new  values  of  h  and  h  .  must  be 

pa  pb 

assumed  and  steps  2)  through  6)  are  repeated. 

In  essence,  finding  the  correct  hpa  and  hpb  represents  solving 
two  nonlinear  equations  in  two  unknowns,  where  the  function  to  be 
evaluated  each  time  represents  the  entire  solution  to  the  seal 
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1J* 


3E 


lubrication  problem.  This  procedure  was  facilitated  by  using  a  non¬ 
linear  root  finding  technique  based  upon  fitting  approximate  planes  to 
the  surfaces  represented  by  the  two  conditions,  and  solving  for  the 
condition  where  the  intersection  of  the  two  planes  satisfies  both  con¬ 
ditions.  This  is  equivalent  to  a  two-dimensional  version  of  the  se¬ 
cant  method  for  finding  a  root.  The  method  works  quite  well  if  one 
can  guess  the  correct  neighborhood  for  the  solution. 

Once  the  solution  above  is  found,  then  the  actual  pressure  to 
which  the  solution  corresponds  is  determined  as  follows. 


megb  =  "mepb  ’ 


mega  ~  mea  m9pa  ’ 


(4-74) 

(4-75) 


where  m0  and  m^  are  computed  from  the  pressure  as  described.  The 
amplitude  of  meg  is 


eg 


* 


T 

ega 


Since 


dr .e 

meg  =  p9  —  s1n  ST 


(4-77) 


from  equation  (4-39),  the  corresponding  gas  pressure  can  be  calculated. 

It  should  be  pointed  out  for  a  static  solution  all  of  the  subscript 
b  terms  vanish  because  of  symmetry,  and  the  solution  involves  only  one 
unknown. 


Partial  Wear  Conditions 

In  some  of  the  results  to  follow,  a  partial  wear  condition  was 
used.  For  a  completely  worn  in  seal,  the  condition  given  by  equation 
(4-28)  for  100  percent  wear  was  used.  For  an  unworn  seal,  the  term 
w(r)  of  equation  (4-66)  was  set  equal  to  a  constant.  Conditions  in 
between  these  extremes  were  also  evaluated  for  static  cases  where  the 
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tilt  and  waviness  were  in  phase.  This  was  done  using  the  approach 
shown  in  Figure  4-11  where  the  fraction  of  interference  between  the 
faces  was  calculated  for  the  two  extreme  deflection  positions.  When 
the  average  height  hQ  has  been  adjusted  to  provide  the  desired  frac¬ 
tion  of  wear,  then  the  wear  function  w(r)  is  defined  as  the  shaded  por¬ 
tion  in  the  figure  that  must  be  removed  to  eliminate  interference. 

Parameter  Values  for  Test  Seal 

Equations  (4-30)  and  (4-31)  show  the  values  of  most  of  the  param¬ 
eters  required  for  the  computer  solutions  for  the  test  seal.  Addi¬ 
tional  parameters  and  information  are  given  in  Figure  4-12.  Several 
of  the  parameter  value  selections  require  some  explanation. 

The  surface  roughness  used  (equation  (4-30))  for  the  computed 
results  is  C  =  0.51  urn  (0.14  pm  CLA).  The  surface  finish  of  the  car¬ 
bon  after  test  was  found  to  vary  from  0.2  to  0.5  pm  CLA,  somewhat 
larger  than  the  value  used  in  the  computed  results.  It  is  thought 
that  the  0.14  pm  CLA  value  used  is  justified  as  follows. 

1)  After  a  running  in  period  no  leakage  of  flat  faced  seals  has 
ever  been  measured.  Leakage  rates  lower  than  0.1  ml/pin. 
have  been  measured  so  that  the  measurement  of  low  level  leak¬ 
age  rates  is  known  to  be  possible.  Based  upon  a  theoretical 
analysis  using  0.14  pm  CLA.the  flat  face  leakage  should  be 
0.1  ml/min.  In  a  parallel  face  seal  the  leakage  rate  is 
determined  primarily  by  the  surface  roughness  (a  cubic  rela¬ 
tionship).  Since  no  leakage  has  been  measured,  it  is  thought 
that  the  effective  roughness  is  even  lower  than  the  assumed 
value. 

2)  Close  examination  of  a  tangential  surface  profile  shows  that 
the  surface  is  characterized  by  peaks  which  are  quite  flat 

on  top  and  valleys  which  are  quite  deep  relative  to  roughness 
which  appears  on  the  flat  peaks.  The  valleys  are  spaced  far 
enough  apart  such  that  most  likely  they  do  not  connect  to 
each  other.  Given  this  situation,  flow  would  be  controlled 
primarily  by  the  roughness  magnitude  associated  with  the 
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31.75 


Jx*  =  1.840-104  mil4  (0.0442  in4) 
Jy*  *  2.46 -TO3  mn4  (0.00591  in4) 

*  -1.11 -103  mn4  (-0.00267  in4) 
09*^  =  5.95-103  mm4  (0.0143  in4) 
e  *  4.85  mn  (0.191  in) 

E  *  20.7  GPa  (3.0-106  lb/in2) 

G  =  9.6  GPa  (1  .25-106  lb/in2) 

EJx 

A  =  7.44 

GJ9 


♦Calculated  assuming  x/r  <<  1  (creates  1%  error). 

♦♦Calculated  numerically  using  torsion  theory  and  finite  differences. 


Figure  4-12.  Test  Seal  Cross  Section  and  Properties. 


flat  peaks.  This  value  amounts  to  approximately  0.1  ym  CLA, 
much  less  than  the  roughness  measured  taking  the  entire  sur¬ 
face  into  account.  This  value  is  approximately  the  roughness 
value  used  for  the  computed  results. 

The  roughness  question  is  very  important.  Table  4-1  shows  the 
equilibrium  operating  condition  for  the  test  seal  at  base  operating 
conditions  for  three  roughness  values.  Torque  varies  a  factor  of  two, 
while  leakage  varies  a  factor  of  ten.  The  effect  of  roughness  on  hy¬ 
drodynamic  load  support  was  thoroughly  analyzed  in  reference  [1]  and 
shown  to  be  quite  significant.  Even  though  a  rational  basis  was  used 
for  the  effective  roughness  used  in  the  computed  results  herein,  it 
is  clear  that  a  more  thorough  analysis  of  the  actual  roughness  is 
needed.  To  this  end,  investigation  of  the  actual  surface  using  the 
previously  developed  flow  model  [7]  is  now  underway. 

Mechanical  friction  coefficient  is  taken  at  0.1.  Based  on  measure¬ 
ments,  this  value  may  be  somewhat  high  as  discussed  later.  However, 
friction  coefficient  varies  over  a  wide  range  depending  upon  conditions, 
so  a  representative  value  is  appropriate.  In  Figure  4-12  the  value  of  E 
is  taken  as  the  reported  value  for  P658-RC  carbon.  Deflection  tests  have 
been  conducted  confirming  this  value  to  within  10  percent.  The  value  of 
6  is  based  upon  assuming  a  Poisson's  ratio  of  0.2  for  the  carbon. 

Comparison  of  Theory  and  Experiment 
Deflection 

Using  the  method  outlined,  solutions  for  deflection  and  pressure 
distribution  were  obtained  for  the  test  conditions  used  for  the  proxim¬ 
ity  probe  tests  outlined  in  Chapter  2.  A  zero  fraction  of  wear  was 
assumed.  For  the  zero  sealed  pressure  case,  only  the  spring  pressure 
acted.  The  tests  were  conducted  under  static  conditions.  Table  4-2 
shows  the  calculated  results  for  the  zero  pressure  case.*  Table  4-3, 
zero  wear,  shows  the  computed  results  for  the  pressurized  case.  Deflec¬ 
tion  at  the  midpoint  of  the  face  (the  location  of  the  proximity  probes) 
was  calculated  using  equation  (4-66).  Figure  4-13  shows  a  comparison  of 

♦Zero  pressure  is  approximated  by  0.007  MPa  (1  psi ) . 
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Figure  4-13.  Comparison  of  Measured  Face  Deflection 
Theory. 
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the  theoretical  deflection  to  the  measured  deflection  given  in  Table 
2-7.  In  the  lower  curve,  where  most  of  the  gas  pressure  caused  waviness 
was  flattened  out,  agreement  is  quite  good.  In  the  upper  curve  where 
very  little  of  the  wave  was  flattened,  agreement  is  good  at  3.5  MPa  but 
begins  to  deviate  at  7  MPa.  Actual  error  at  this  level  is  about  10  per¬ 
cent.  Considering  the  numerous  approximations  and  uncertainties  in  both 
theory  and  experiment,  this  error  appears  to  be  quite  acceptable.  The 
results  shown  in  Figure  4-13  show  that  the  deflection  theory  used  is  in 
all  probability  valid  for  the  purpose  of  this  analysis. 

Zero  Speed  Results 

Because  of  the  relatively  large  gaps  that  develop  in  the  seal 
under  the  combined  waviness  and  tilt  conditions,  and  because  no  reduc¬ 
tion  in  friction  with  increasing  speed  was  observed  in  the  tests,  it 
was  originally  thought  that  hydrodynamic  effects  were  playing  a  minimal 
role  in  seal  performance  and  that  a  static  analysis  might  provide  an 
adequate  prediction  of  performance.  Thus,  the  results  in  Table  4-3 
were  calculated.  Oifferent  wear  fractions  were  evaluated  to  attempt 
to  simulate  the  type  of  wear  observed  and  shown  in  Figure  2-14.  The 
most  significant  feature  of  these  results  indicating  disagreement  with 
experiment  is  that  the  torque  changes  very  little  between  3.5  and  7.0 
MPa  gas  pressure.  This  is  contrary  to  experimental  results  which  show 
that  torque  decreases  as  gas  pressure  increases.  Also,  the  calculated 
leakage  is  too  high  at  least  for  the  100  percent  wear  case.  These 
results  do  have  some  use  however  in  that  they  show  how  the  seal  might 
be  expected  to  behave  at  low  speeds. 

Dynamic  Results  -  Gas  Pressure  Effect 

Table  4-4  shows  the  results  for  different  gas  pressures  calculated 
using  the  model  described  and  assuming  100  percent  wear.  Note  that  all 
of  the  various  deflection  terms  discussed  show  up  in  the  table  and  that 
the  gas  pressure  is  corrected  according  to  the  procedures  described. 

The  theoretical  results  of  Tables  4-3  and  4-4  are  shown  compared  to 
experimental  results  of  Table  2-2  in  Figures  4-14  and  4-15.  Consid¬ 
ering  first  Figure  4-14  for  leakage,  the  1800  RPM  theoretical  results 
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Figure  4-15.  Wavy  Seal  Torque. 


clearly  follow  the  experimental  trend,  and  considering  the  scatter  in 
the  leakage  results,  can  be  considered  to  give  adequate  agreement.  The 
zero  speed  theoretical  result  clearly  overpredicts  leakage.  This  points 
out  that  hydrodynamic  effects  are  significant,  even  though  the  seal  does 
not  operate  in  a  classically  hydrodynamic  manner.  The  added  deflection 
components  and  load  support  provided  by  the  effect  of  hydrodynamic  pres¬ 
sure  buildup  and  cavitation  clearly  alter  the  seal's  behavior  relative 
to  the  static  case.  Leakage  is  smaller  in  the  dynamic  case  because  the 
net  waviness  is  in  fact  smaller  than  in  the  hydrostatic  case. 

Now  considering  Figure  4-15,  the  1800  RPM  theoretical  results 
clearly  predict  the  trend  in  torque,  that  of  a  gradually  decreasing 
torque  with  increasing  pressure.  The  static  results  do  not  predict 
this  trend  as  mentioned  previously.  Since  the  theoretical  prediction 
of  torque  requires  that  the  coefficient  of  friction  between  carbon  and 
W-C  be  assumed,  and  the  value  assumed  in  the  calculation  was  not  actu¬ 
ally  measured  (since  it  varies  widely  depending  on  conditions),  the 
torque  prediction  based  on  a  lower  coefficient  of  friction  (0.06)  is 
also  shown.  This  theoretical  result  intersects  the  experimental  results. 

Thus,  so  far,  prediction  of  both  leakage  and  torque  by  the  theore¬ 
tical  model  appears  to  be  reasonable  and  suggests  that  the  model  may  be 
satisfactory.  However,  further  verification  is  desirable  and  the  wear 
profiles  are  considered  next. 

Wear  Profiles 

Using  the  test  results  for  100  h  wear  tests  at  3.45  MPa  water  pres¬ 
sure  and  6.9  MPa  gas  pressure  tabulated  in  Tables  2-2  and  2-3,  average 
wear  slopes  were  obtained  from  the  surface  traces  shown  in  Figure  2-14. 
These  are  tabulated  in  Table  4-5.  Based  on  the  theoretical  wear  profile 
for  100  percent  wear,  the  theoretical  slopes  were  computed  and  are  also 
shown  in  Table  4-5.  The  comparison  for  the  outside  slope  is  quite  good. 
The  inside  slope  does  not  compare  as  well,  however,  these  slopes  are 
rough  estimates  and  good  agreement  would  not  be  expected. 

Perhaps  more  important  than  the  values  of  the  slopes  is  the  fact 
that  the  shape  of  the  theoretical  wear  surface  shown  in  Figure  4-16  does 


Table  4-5 


Wear  Slope  Comparison 

p _a _  =  6.9  MPa  Pu  n  =  3.45  MPa  1800  RPM 
100  h 


L.H 

Outside 

L.H 

Inside 

R.H 

Inside 

R.H 

Outside 

21 

710*x  10"6 

390  x  10‘6 

400  x  10"6 

875  x  10"6 

34 

670  x  10"6 

670  x  10"6 

860  x  10"6 

500  x  10"6 

39 

1200  x  10"6 

1300  x  10“6 

— 

600  x  10"6 

Slope 

Experimental 

Theory 

(1800  RPM  cas^ 

Average  of  Outside  Slopes 

760  x  10‘6 

680  x  10'6 

Average  of  Inside  Slopes 

730  x  10"6 

430  x  10'6 

★ 

Data  taken  from  wear  profiles  in  figure  2-14. 
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match  the  general  shape  of  the  wear  profiles  shown  in  Figure  2-14.  It 
is  clear  that  the  double-sided  touching  observable  in  Figure  2-14  is 
predicted  by  the  theory  in  Figure  4-16.  Given  this  theoretical  wear 
shape,  it  is  also  probable  in  Figure  2-14  that  in  tests  34  and  39,  the 
100  percent  wear  condition  was  not  achieved  (see  also  wear  test  results 
and  compare  between  tests  21,  34,  and  39). 

Figure  4-16  illustrates  one  final  limitation  of  the  static  anal¬ 
ysis.  At  100  percent  wear,  the  wear  profile  is  quite  different  from 
what  is  observed.  At  50  percent  wear  however,  the  profile  is  similar. 
However,  it  is  clear  that  the  1800  RPM  case  predicts  the  shape  better 
than  the  static  result. 

Speed  Effects 

Table  4-6  shows  the  theoretical  prediction  for  100  percent  wear  at 
various  speeds.  The  unexpected  result  shown  here  is  that  torque  in¬ 
creases  with  speed  even  though  it  is  known  that  hydrodynamic  effects 
are  present.  The  reason  for  this  is  that  the  lubrication  behavior  is 
interactive  with  deflection.  As  speed  changes  the  pressure  distribution 
changes  and  changes  the  deflection.  Given  the  sensitivity  of  perform¬ 
ance  to  the  deflection  or  net  waviness  and  tilt  described  previously, 
it  is  quite  reasonable  to  assume  that  as  speed  increases  for  this  parti¬ 
cular  seal,  the  net  waviness  and  tilt  are  altered  to  a  less  favorable 
configuration  for  load  support  than  at  the  lower  speeds.  Therefore, 
what  initially  appears  to  be  a  surprising  and  contradictory  trend  is 
quite  explainable  in  the  context  of  the  total  interactive  mechanism  of 
the  seal . 

Theoretical  seal  leakage  and  torque  are  compared  to  experimental 
results  in  Table  4-7.  The  experimental  torque  shows  the  same  trend 
as  the  theoretical  torque,  an  increase  with  speed  beyond  a  certain 
speed;  agreement  for  leakage  is  not  good.  It  is  assumed  that  this 
poor  agreement  stems  from  the  fact  that  the  speed  tests  were  short 
term  tests  and  the  seal  faces  did  not  approach  100  percent  wear. 
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6.90  M 


Figure  4-16.  Theoretical  Carbon  Face  Wear  Profiles. 


Table  4-7 
Speed  Effect 

Comparison  of  Experimental  and  Theoretical  Results 
pu  n  =  3.45  MPa,  p„ae  =  6.90  MPa,  100%  Wear 

H0U  yaS 


Speed 

Experimental 

Theoretical 

Tq* 

N*m 

Q* 

ml  /min 

Tq 

N-m 

Q 

ml /min 

900 

1.8 

0.2 

3.1 

10.4 

1800 

1.5 

1.9 

4.5 

4.4 

2700 

3.2 

4.5 

6.1 

2.6 

3600 

5.0 

1.8 

6.6 

2.8 

♦Average  of  values  in  Table  2-6 
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Chapter  5 

Summary,  Conclusions,  and  Recommendations 

In  this  chapter  both  specific  and  general  conclusions  are  drawn. 
Experimental  Results 

1)  Test  Apparatus  -  The  test  apparatus  described  with  the  wavi¬ 
ness  mechanism  does  provide  the  capability  of  measuring  the  effects  of 
waviness  on  seal  performance.  No  major  operational  limitations  have 
been  found.  While  the  waviness  device  works  satisfactorily  on  this 
experimental  apparatus, there  is  wear  on  the  waviness  0-rings  which 
would  limit  life  on  commercial  applications.  Other  means  of  applying 
waviness  must  be  developed  in  order  to  utilize  the  concept  commercially 
One  aspect  that  does  need  improvment  is  finding  the  zero  torque  level. 
An  improvement  would  incorporate  a  means  of  pulling  the  seal  faces 
apart  after  the  vessel  was  sealed. 

2)  Leakage  -  The  use  of  the  leakage  measuring  device  described 
greatly  improves  the  quality  of  the  leakage  data.  It  does  provide 
leakage  measurements  over  a  very  wide  range  of  leakage. 

3)  Computer  -  The  computer  data  acquisition  system  greatly  im¬ 
proves  the  quality  of  data  which  can  be  collected  and  stored  for  later 
analysis.  The  computer  also  provides  the  ability  to  both  control  the 
test  apparatus  and  monitor  performance  from  a  safety  standpoint. 

4)  Flat  Face  Results  -  Torque  is  high  compared  to  the  wavy  seal 
and  fluctuates  with  time.  Leakage  is  zero.  Torque  does  show  an  ini¬ 
tial  increase  with  time  indicating  that  thermal  coning  causes  a  con¬ 
verging  radial  taper  which  wears  off. 

5)  Wavy  Seal  Results  -  Torque  level  is  about  one-fifth  of  that 
for  a  flat  faced  seal.  Leakage  increases  with  increasing  waviness. 
Torque  level  is  quite  steady  except  for  a  pulsing  corresponding  to  the 
superposition  of  the  third  and  second  harmonic  waves. 
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6)  Wear  -  Wear  rates  are  from  four  to  forty  times  lower  for  the 
wavy  seal  compared  to  the  flat  faced  seal.  A  better  method  of  measuring 
wear  is  needed  to  obtain  more  consistent  results. 

7)  Radial  Taper  -  A  convergent  radial  taper  reduces  friction 
torque  as  expected.  Results  show  that  the  radial  taper  wears  away 
with  time  as  would  be  expected  [8]. 

8)  Sealed  Pressure  -  Tests  show  that  the  waviness  imposed  de¬ 
creases  with  increasing  sealed  pressure.  That  is,  the  seal  waviness 
flattens  out  with  increased  load. 

9)  Speed  -  For  the  wavy  seal  torque  increases  with  increasing 
speed.  The  reason  for  this  behavior  appears  to  be  that  the  seal  deflec¬ 
tion  changes  to  a  less  favorable  configuration  as  the  hydrodynamic 
effects  become  more  significant.  These  tests  need  to  be  repeated  to 
make  certain  the  seal  is  worn  in  at  each  operating  condition. 

10)  Proximity  Probe  -  Tests  show  that  it  is  possible  to  measure 
changes  in  carbon  face  displacement  relative  to  the  tungsten  carbide. 

11)  Wear  Profile  -  The  wear  profile  on  the  wavy  seal  shows  that 
the  seal  faces  tilt  as  well  as  wave.  This  requires  a  more  comprehen¬ 
sive  theory  to  properly  analyze  behavior  than  presented  previously. 

The  wear  profiles  for  flat  face  tests  clearly  show  evidence  of  thermal 
coning  as  expected. 

12)  Surface  Roughness  -  So  far,  no  relationships  between  seal  per¬ 
formance  and  surface  roughness  has  been  found.  However,  tests  where  sur¬ 
face  roughness  has  been  deliberately  altered  have  not  yet  been  run. 

Thermoelastic  Instability 

1)  The  model  developed  shows  that  ring  deflections  have  a  pre¬ 
dominant  effect  over  semi  infinite  plate  deflection  at  low  n  and  must 
be  considered  as  a  part  of  a  model  to  describe  thermoelastic  insta¬ 
bility  in  seal  rings. 


2)  The  effect  of  thermal  contact  resistance  between  the  faces  on 
critical  speed  is  large.  However,  the  model  for  heat  transfer  between 
the  faces  used  herein  requires  further  investigation.  For  example,  it 


was  assumed  that  all  heat  is  generated  at  the  middle  of  the  interface. 
Also,  the  effect  of  tangential  mass  transfer  of  the  sealed  fluid  might 
be  significant. 

3)  The  effect  of  convection  losses  on  $  is  also  large.  This 

effect  requires  further  evaluation.  Preliminary  results  suggest  that 

a  -*■  »  even  with  H  small, 
c  00 

4)  Although  some  of  the  critical  speeds  and  number  of  waves  pre¬ 
dicted  agree  with  the  experimental  observations  reported,  in  view  of 

2)  and  3)  above  this  agreement  should  be  regarded  as  coincidence.  How¬ 
ever,  the  agreement  does  encourage  the  refinement  of  the  model. 

5)  The  combined  stability  curves  of  Figure  3-6  are  only  to  be 
considered  as  first  approximations.  The  curves  are  based  on  the  sum 
of  two  limiting  cases.  Further  refinement  will  require  that  elastic 
and  thermal  deflection  terms  be  computed  using  finite  element  numeri¬ 
cal  techniques.  It  is  thought  that  such  a  refinement  would  result  in 
a  more  accurate  stability  curve,  particularly  for  the  nose  portion. 

6)  At  this  point  in  the  development  of  thermoelastic  theory  for 
rings  there  is  a  definite  need  for  further  experimental  results  to 
determine  if  the  basic  shape  of  the  stability  curve  of  Figure  3-6  is 
correct. 

Wavy  Seal  Analysis 

1)  Reynolds  Equation  for  Rough  Surfaces  -  Since  there  is  still 
controvery  over  how  to  model  the  effect  of  roughness  on  the  Reynolds 
equation,  certain  simplifying  assumptions  were  made  herein  to  avoid 
the  problems  associated  with  selecting  one  model  over  another  model. 

2)  Waviness  Model  -  A  new  seal  waviness  lubrication  model  was 
developed.  The  fact  that  there  is  a  tilt  in  addition  to  the  wave 
necessitated  this  model.  This  model  is  more  general  but  not  as  accu¬ 
rate  because  the  precise  cavity  boundary  conditions  could  not  be  used 
in  conjunction  with  this  more  general  case.  The  new  model  accounts 
for  alteration  of  the  face  shape  brought  about  by  tilting  of  the  faces 
and  wear. 
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3)  Film  Thickness  Shape  -  For  the  more  general  case  of  waviness 
combined  with  tilt,  it  is  shown  that  the  film  thickness  shape  becomes 
complex.  Around  the  circumference  of  the  seal  there  are  in  general 
regions  where  the  faces  are  radially  parallel,  converging,  and  diverg¬ 
ing. 

4)  General  Solution  -  Using  the  model,  a  family  of  solutions  was 
made  for  various  combinations  of  waviness  and  tilt.  It  was  found  in 
general  that  torque  and  leakage  are  lowest  when  waviness  and  tilt  are 
in  phase.  The  combination  of  waviness  and  tilt  provides  lower  torque 
and  leakage  than  if  the  faces  are  simply  wavy.  This  occurs  because  a 
favorable  hydrostatic  pressure  distribution  results. 

5)  Optimum  Conditions  -  When  tilt  is  such  that  seal  face  contact 
goes  between  the  extremes  of  contact  at  the  inside  radius  in  conjunc¬ 
tion  with  a  converging  taper  to  flat  faced  contact,  an  optimum  condi¬ 
tion  giving  minimum  torque  results.  Both  hydrodynamic  and  hydrostatic 
pressure  distributions  become  favorable  to  provide  maximum  load  support. 
Leakage  is  quite  small  because  there  exists  a  continuous  ring  of  contact 
around  the  seal . 

6)  Ring  Deflections  -  It  is  shown  that  the  predominant  deflection 
in  the  seal  ring  is  caused  by  waviness  gas  pressure  and  pressure  on  the 
face.  Deflection  equations  are  developed.  There  are  twelve  signifi¬ 
cant  components  of  deflection,  six  sine  and  six  cosine  components.  Of 
the  six,  there  is  waviness  due  to  gas  pressure,  moment  of  the  face 
pressure,  and  the  face  pressure  itself.  Each  of  these  forces  also 
produces  corresponding  tilt  terms. 

7)  Nonparallel  Faces  -  The  tilt  term  and  waviness  term  are  not 
in  general  proportional .  This  means  that  even  though  waviness  may  be 
completely  flattened,  tilt  will  not  be  flattened.  The  implications  of 
this  is  that  in  general  the  wavy  seal  will  operate  with  a  large  com¬ 
ponent  of  tilt  in  conjunction  with  waviness. 

8)  Test  Seal  Solution  -  The  solution  for  the  test  seal  requires 
the  simultaneous  solution  for  two  deflection  terms.  This  presents 
some  numerical  difficulty  but  a  method  for  doing  this  was  found  so 
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that  the  performance  of  the  actual  seal  could  be  predicted. 

9)  Surface  Roughness  -  Theoretical  results  show  that  surface 
roughness  has  a  large  effect  on  seal  torque  and  leakage.  It  is  shown 
how  the  effective  surface  roughness  is  much  lower  than  the  measured 
surface  roughness  because  of  the  nature  of  the  surface.  It  is  pro¬ 
posed  that  further  theoretical  work  be  done  to  model  the  surface 
roughness  effects  for  the  actual  surface  so  that  the  assumptions  made 
can  be  supported  or  rejected. 

10)  Comparison  of  Theory  and  Experiment  -  Wavy  Seals 

a)  Deflection  -  Calculated  deflection  was  compared  to  deflection 
measured  using  a  proximity  probe.  Agreement  was  quite  good  thus  val¬ 
idating  the  deflection  theory. 

b)  Static  Results  -  Theoretical  results  based  on  a  static  anal¬ 
ysis  (0  RPM)  do  not  compare  well  to  experiment.  This  indicates  that  in 
spite  of  the  large  taper  in  the  seal,  hydrodynamic  effects  must  still  be 
significant.  While  the  static  analysis  is  much  simpler  to  use,  predic¬ 
tions  are  inadequate. 

c)  Dynamic  Results  -  Using  the  model  developed,  a  comparison  of 
theoretical  and  experimental  results  show  good  agreement  for  leakage 
and  a  good  prediction  for  the  trend  for  torque. 

d)  Wear  Profiles  -  The  model  predicts  the  shape  of  the  wear  pro¬ 
file  observed. 

e)  Speed  Effects  -  The  model  predicts  the  observed  increase  in 
torque  with  speed. 

General 

1)  The  comparison  of  theoretical  to  experimental  results  is  gen¬ 
erally  quite  good.  While  tnere  are  known  areas  where  the  theoretical 
model  requires  improvement,  the  agreement  observed  suggests  the  basis 
for  the  model  is  valid. 

2)  Hydrodynamic  effects  do  occur  in  the  wavy  face  seal.  This 
conclusion  is  supported  by  the  fact  that  only  after  hydrodynamic  fea¬ 
tures  were  included  in  the  model  was  agreement  between  theory  and 
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experiment  satisfactory. 


4)  Waviness  reduces  friction  and  wear  in  a  face  seal.  This  is 
now  known  experimentally  as  well  as  theoretically.  Given  the  results 
herein,  it  is  quite  likely  that  water  seals  can  be  designed  which  have 
on  the  order  of  ten  times  the  life  of  a  comparable  flat  faced  seal. 

5)  The  moving  waviness  concept  does  distribute  the  wear  and  does 
preserve  the  wavy  shape  as  predicted.  The  moving  wave  must  be  used  to 
benefit  from  wav i ness  in  the  long  term. 

6)  The  discovery  of  the  optimum  combinations  of  waviness  and  tilt 
points  the  way  toward  designing  a  moving  waviness  seal  having  even  better 
performance  (i.e.,  lower  torque  and  leakage)  than  those  tested  herein. 
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Computer  Control  and  Data  Acquisition  Program 


134 


S  t  CONTROL  PROGRAM 
10  1 

IS  l 

20  STANDARD 
25  OPTION  BASE  1 

20  COM  Zeco,Inlbav, Relay, FlagS,?h2e>,Rpm,Tqmax,Toaax,TfaaxfTbmax,Tvmax 

35  COM  T.Tineref, SHORT  Thout,Toout,Twout,R£AL  Lkold,Tlkold,Tllt(2Q0>  ,Thce£,Nl)c 

40  COM  160, I60on,U2,Ngo, SHORT  TiBe,T£out,Tqout,Lkout,Zil,Z12,Z13,Z14,Vibset,N 

7ib,Tvib,0ays(12) 

45  SBORT  Tiaeovec,Tqover,Tmpover,LJ(over ,Filedu» 

50  DIM  Tiaeover ( 300 ) ,Tqover(300) .Tapover (300) .LJcover  (300) ,Oa(10) 

55  DIM  A$< 40) 

60  READ  AS (*) 

65  OATA  ?1,F2,F3,F4,F5,F6,F7,F8,P9,P10 

70  DATA  F11,F12,F13,F14,?15,F16,F17,F13,F19,F20 

75  DATA  F21,F22,F23,F24,F25,F26,F27,F28,F29,?30 

80  OATA  F31,F32,F33,F34,F35,F36,F37,F38,F39,F40 

85  REAO  D ay s ( * ) 

90  OATA  0,31,60,91,121,152,182,213,244,274.305,335 
95  ON  ERROR  GOSUB  240 

100  t  Be lay "DECIMAL ( 20000 ) 

105  OUTPUT  5  WHS  USING  *  4, W , *-4000, Relay 
110  CONTROL  MASK  5; 128 
115  PLOTTER  IS  7,5,*9372A" 

120  LIMIT  0,270,0,190 
125  ON  KEY  40,14  GOSUB  10 
130  ON  KEY  11,10  GOSUB  11 
135  ON  KEY  42,10  GOSUB  12 
140  ON  KEY  43,10  GOSUB  13 
145  ON  KEY  44,10  GOSUB  14 
150  ON  KEY  45,10  GOSUB  IS 
155  ON  KEY  46,10  GOSUB  16 
160  Ngo-7 
165  GOTO  250 
170  I Is  Ngo-1  !  RUN 
175  RETURN 

180  12:  Ngo-2  !  START 

135  RETURN 

190  13:  Hgo«3  !  AXES 

195  RETURN 

200  14:  Ngo-4  1  REZERO 

205  RETURN 

210  15:  Ngo*S  ICONRUN 

215  RETURN 

220  16:  Ngo-6  !  TEST  END 

225  RETURN 

230  10:  Mgo-7  S  IDLE 

235  RETURN 
240  PRINT  ERRMS 
245  RETURN 

250  INPUT  "Teat  no?",Testno 

255  INPUT  "Typetest?"  .TypetestS 

260  INPUT  "Date?*, oateS 

265  INPUT  *Ph2o?",Ph2o 

270  INPUT  "Pgas?" , Pgaa 

275  INPUT  "Rpa?" ,Rpn 

230  INPUT  "Twmax?" ,Twmax 

235  INPUT  "Tqraax?" ,Tqmax 

290  INPUT  "Toilmax?" ,Tomax 

295  INPUT  "Ttacenax?" ,T£»ax 
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300  INPUT  "Thousi.ngmax7"  ,T6max 

30 S  INPUT  ‘First  fils  for  esst  or  surrtnt  fils  for  continuation?" ,Nfila 

310  ASSIGN  AS(Nfils)  TO  U 

313  GOTO  340 

320  Fils:  Ufils-Nfile-rl 

323  ASSIGN  AS(Nfils)  TO  II 

330  PRINT  ‘FILE  l-"jNfils 

333  RETURN 

340  INPUT  ‘IS  THIS  A  CONTINUATION?" ,Y$ 

345  IP  YS-*Y"  THEN  360 
330  St-0 
355  GOTO  475 

360  INPUT  "REFERENCE  TIME-?" , St 
363  PRINT  "•  CONTINUATION  •“ 

370  PRINT  "REFERENCE  TIME-"jSt 
375  INPUT  "Tlkold?"  ,Tllcold 
380  INPUT  "Lkold?" ,Lkold 
385  Typs-TYP(-l) 

390  IF  Typs-2  THEN  465 

395  IF  (Typs-L)  OR  (Typs-6)  THEN  405 

400  IF  (Typs-3)  OR  CTyps-4)  OR  (Typs-7)  THEN  435 

4QS  READ  lljValus 

410  FOR  I-L  TO  8 

415  Da(  I) -Os ( 1*1) 

420  NEXT  1 
425  0a(9)-Valus 
430  GOTO  385 
435  Na-6 

440  IF  0a(5)—  2000  THEN  Na-1 
445  Tioeold-Oa(Na) 

450  Tapold-Oa(Na+l) 

455  Tqold-Oa(Na+21 
460  GOTO  480 
465  READ  lljValusS 
470  GOTO  335 

475  PRINT  41jTestno,OatsS,TypstsstS,Ph2o,Pqas,Rp*,Twaa*,Tq«a*tToaax.Tfaax,Tli*ax 
, Nf ila 

480  PRINT  "TEST  NO  -".-Tsstno 

485  PRINT  "DATE  -".-OatsS 

490  PRINT  "TYPETEST  -"sTypstsstS 

495  PRINT  "Ph2o  -*;PH2o 

500  PRINT  "Pgas  -";Pgas 

505  PRINT  "Rps  -";Ppm 

510  PRINT  "Twnax  -'jTwcax 

515  FRINT  "Tqniax  -"sTqraax 

520  PRINT  "Toiaax  -";Tonax 

525  PRINT  "Tfnax  «";Tf»ax 

530  PRINT  "Thnax  «" ; Tftna x 

535  PRINT  "FIRST  PILE  «";NfUs 

540  FIXED  2 

545  ON  END  II  GOSUB  Fils 

550  INPUT  "00  YOU  KANT  TO  CALIBRATE  TORQUE?", YS 
535  IF  YS— " Y"  THEN  Zc 

56u  INPUT  "WHAT  IS  CALIBRATION  CONSTANT  IN  IN-LB/KV" , Inlbnv 

505  GOTO  Za 

570  Zc:  CALL  Calibrate 

575  Zd:  PRINT  "IN-LB/WV-"  ;  InlbtT.v 

560  INPUT  "DC  YOU  KANT  TO  FI.NO  ZERO?" , YS 

565  1!  YS-"N  *  THEN  Zb 

590  CALL  Zsrotub (Thcsf ) 

595  GOTO  Za 

600  Zb i  INPUT  "KHAT  IS  ZERO  TOFQUE  IN  IN-LB?", Zsrc 

C 05  isrpur  "'..HAT  IS  f.SF  TEMP?" , Th ref 

61ci  Za:  s'FINT  "Zaco  »":Zaro 

•515  PRINT  "T.’.csf  -"jtr.rsf 

5 * 5  INP'il  "CN5CK  ZCIND  AGAIN  ,YZ 
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czi  :?  ys-*y*  tuc:;  590 

63G  3  TV. .  3 3  3 

335  INPUT  "00  YOU  WANT  TO  DRAW  COORDS?", 25 
$40  IP  2$-*N*  THEN  Ca 
645  FlagS-*Idla* 

650  CADE  Coecd(Tiaamin,Tieanax,TEqmin,Tcqmax,LkBax,Tnpmax,T!apmin) 
655  GOTO  Cb 

660  Ca:  INPUT  "Tiaanin?* ,Ti*aain 

665  INPUT  *Tia«max?* .Tlacoax 

670  INPUT  *Trqain?“,Trqoin 

675  INPUT  ’Trqmax?* ,Trqaax 

680  INPUT  *Lksax?" ,Lkaax 

685  INPUT  *Tmp«ax?* .Tapoax 

690  INPUT  ’Tapain?* ,Tapain 

695  Cb*  A-10+120/26*2 

700  LOCATE  A, 130,10, 97 

70S  ON  INT  *5,14  CALL  16940 

710  l 

715  l 

720  Idlat  t 

725  D-PNRoff(l) 

730  D-FNRo££ ( 2) 

735  0-FNRof £ ( 4) 

740  Flags-* Idla* 

745  l 

750  Idlaloop:  DIS?  ‘CONTROL  PROGRAM  READY* 

755  ON  Ngo  GOTO  Run, Scarc,Axes , Razaro ,Conrun,Tastend, Idlaloop 
760  l 
765  ! 

770  Run:  1 
775  Vibaae-0 
780  I£ic»t-0 
785  T-0 

790  1600-0 

795  160-0 

800  160on-0 

805  112-0 

810  Nlk-0 

815  IF  StOO  THEM  875 
820  CALL  Tiae(Thrs, TinaS) 

825  Tiaacef-Thcs 

83Q  PRINT  "TEST  START  TIME-* ;Tiaa$ 

835  PRINT  ’REFERENCE  TIME-* » Tiaatbf 

840  Tlknld-Tlanain 

845  Tlkplot-Tiaaain 

850  Tianold-Timnain 

855  Lkold-0 

860  Tqold-0 

865  Tapold-Taomin 

870  GOTO  890 

875  Tiaaraf-St 

880  CALL  Tiaa  (Thrj,TiaaS) 

385  Tlkold-Cbes-^Tiaacaf 
390  Ndata-0 
895  Ioploc-0 
900  Ilkcaa-0 
905  0-PNRon(4) 

910  Flags-* Run* 

915  Iloop-1 
920  Iovac-0 
925  GOTO  Runloop 
930  1 

935  1 

940  Stact:  1 
945  Vlbaat-O 
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¥30  0-f NROt £ ( 4 ) 

955  D-FNRon(l) 

960  ! 

965  1 

970  Startloop:  PlaqS-'Star t" 

975  CALL  Tint 

960  015?  "Startioop* 

965  ON  Nqo  GOTO  Run,Staetloop,Axa»,Razaro,Conrun,Tastand, Idla.Raplot 
990  i 
99S  I 

1000  RunlOOp !  1  INFINITE  LOOP 

1005  OISP  'Runloop* 

1010  CALL  Tint 

1015  IP  16001  THEN  Rloopand 

1020  IT  I60on<>l  THEN  Rloopand 

1025  ir  Ifirst-0  THEN  1390 

1030  ON  Hoop  GOTO  Ll,L2,L3,L4,LS,L6,L7,L6 

1035  1 

1040  LI:  PRINT  *l;Tiaa,TJout,Tqout,LkOut 
104 5  Ndata-Ndata+1 
1050  Iloop«Iloop+l 
1055  GOTO  Rloop 
1060  1 

1065  L2:  SCALE  Tinaain.Tiaaaax ,Trqain, Ttqaax 

1070  MOVE  Tiaaold ,Tqold 

1075  Iloop-Iloop+l 

1080  GOTO  Rloop 

1085  t 

1090  L3 :  IF  Tina>Tinamax  THEN  Ml 
1095  IF  Tqout>Tcqaax  THEN  Tqout-Trqnax 
1100  Hp«Tqout"2*PI» ( 1/12) •Rpo/33000 
1105  IF  Hp<5.5  THEN  1125 

1110  PRINT  "SHUTDOWN  -  HORSEPOWER  IS*;Hp  ■ 

1115  PRINT  "TIME-'sTiaaS, “Lkold-*  sLkold,  "Tlkold-'jTlkold 
1120  Nqo>7 

1125  IF  Tqout<Trqain  THEN  Tqout-Trqain 

1130  PLOT  Tina,Tqout«-l 

1135  PENUP 

1140  Iloop-Iloop+1 

1145  GOTO  Rloop 

1150  l 

1155  L4 :  SCALE  Tia«nin,Tia«aax,TapBin,Tapnax 

1160  MOVE  Timaoid.Tapold 

1165  Iloop-Ilaop+l 

1170  GOTO  Rloop 

1175  ! 

1160  L5:  IF  Tfou&>Tapmax  THEN  TJout-Tnpnax 
1185  I?  Tfouc<Tapmin  THEN  Tfout-Tapnin 
1190  PLOT  Tiaa.Tfout.-l 
1195  PENUP 

1200  move  Tiaaold, Two Id 
1205  PLOT  Tina,Twouc,-l 
1210  PENUP  / 

1215  Iloop-lloop+l  . 

1220  GOTO  Rloop 
1225  1 

1230  Ml:  IF  IOvac>3CG  THEN  GOTO  1260 
1235  Iov«r*lovac*’l 
1240  Timeovac ( lovar) -Time 
1245  Tqovar (lovar) -Tqout 
1250  Tapover (lovar ) -Tfcut 
1253  Lkovat ( lovar ) -Lkout 
1260  180-0 
1265  1600-160041 
1270  Iloop-l 
1275  GOTO  Rlooo 


1280  t 

1285  L6:  Iloop-Iloop+1 
1290  Noplot-0 

1295  SCALE  TineminrTi:n«max(0,Lknax 
1300  IF  Tlkplot<7iaeain  THEN  Nop  lot-1 
1305  HOVE  Tlkplot,Lkold 
1310  GOTO  Rloop 
1315  I 

1320  L7 :  Iloop-Iloop+1 

1325  IF  Lkout>Lkaax  THEM  Lkout-Lkaax 

1330  IF  Lkout— 1000  THE2I  M2 

1335  IF  Noplot-1  THEM  GOTO  1350 

134G  PLOT  Tlaa, Lkout, -1 

1345  PENUP 

1350  Tlkplot-Tlas 

1355  Lkold-Lkout 

1360  : 

1365  M2:MOVE  Tiaaoax , Lkmax 

1370  Timsold-Tiae 

1375  Tqold-Tqoue 

1380  Tapold-Tfout 

1385  I600-I600+1 

1390  160-0 

1395  IF  If icst-0  THEN  Ifirst-1 
1400  Iloop-1 
1405  GOTO  Rloop 
1410  ! 

1415  L8s  1 

1420  Piledum— 2000 

1425  PRINT  #1; F iled  urn,  Tima  , Twout , Toout ,Thout 
1430  160-0 
1435  Iloop-1 
1440  GOTO  Rloop  ‘ 

14  45  ! 

1450  Rloop i  I60on-0 
1455  IF  1600-10  THEH  160-1 
1460  IF  1600-10  THEH  IlOOp-6 
1465  IF  1600-10  THEN  1600-0 
1470  ! 

1475  Rloopcnd:  1 
1480  IF  160-1  THEN  1495 
1485  IF  I60on«0  THEN  1495 
1490  IF  IoplOt-1  THEN  Oplot 

1495  ON  Hgo  GOTO  Runloop, Start, Axa3,Raxeco,Canrun,Ta3t«ndf Idla 
1500  l 
1505  1 

1510  Axss:  FlagS-’Star t* 

1515  O-FHRof f ( 4) 

1520  INPOT  "ARE  TOO  READY  TO  DRAW  C0CR0.?*,Y 

1525  CALL  Coocd(Tineain,Tia«aax,Trqain,Trqaax,Lknax,Tnpaax,Tftpmin) 

1530  Nqo-2 

1535  GOTO  Stack 

1540  t 

1545  t 

1550  Concun:T-0 
1555  1600-0 
1560  160-0 
1565  I60on-0 
1570  112-0 
1575  Nlk-0 
1580  If lest— 0 
1385  FlaqS-*Rua" 

1590  Iloop-1 

1595  CALL  Tir.a  (Thcx,Tia«S) 

1600  Tlkold-Thrs-Tlxacaf 
1505  0-FHRon ( 4) 
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1610  IF  lomOO  THEN  Pi 
1615  Ngo»l 
1620  Ioplot-0 
162S  GOTO  Runloop 

1630  PI:  t  ASSUMES  TIMEMMIN  HAS  SEEN  SHIFTED  SO  DATA  PITS 

163S  I  a  ax- lover 

1640  Iovae-0 

1645  I op lot- 1 

1650  Ngo-1 

1655  Iplot-1 

1660  Igo-1 

1665  M-l 

1670  Tiaaold-Tiaaovac(Iflax) 

1675  Tqold-Tqovar ( I max ) 

1680  Tapold-TBpovat (Iaax) 

1685  FOR  I-Iaax  TO  1  STEP  -1 
1690  J-I 

1695  IP  Lkovar(J)  0-1000  THEN  ML5 
1700  NEXT  I 
1705  GOTG  M16 

1710  M15:  Tlkplot-Tiaeovar (J) 

1715  Lkold*Lkover  (0) 

1720  M16:  i 

1725  GOTO  Runloop 

1730  Op lot:  ON  Igo  GOTO  S1.S2.S3 

1735  SI:  SCALE  Tiaesin.Tiaanax.Trqmin.Trqmax 

1740  IF  Tlaaovat (M) <Tla«aln  THEN  GOTO  1960 

1745  I 

1750  Tqdua-Tqovac ( M) 

1755  IP  Tqdua>Tqaax  THEN  Tqdua-Tqaax 
1760  MOVE  Tiaaovac (M) .Tqdua 
1765  Tqdun-Tqovac (M+l) 

1770  IP  Tqdum>Tqaax  THEN  Tqdua«Tqmax 
1775  PLOT  Tiaaovae (M+l) , Tqdua, -1 
1780  PENU7 
1785  Igo-Igo+1 
1790  GOTO  1495 

1795  S2:  SCALE  Timaain.Tiaaaax.TBpain.TBpmax 
1800  Tapdua-Tapcvac (K) 

1305  IF  TapduaJTapaax  THEN  Tmpdum-Tmpmax 

1810  IP  Tapdua<Trapmin  THEN  Tmpdua-Tnprain 

1815  MOVE  Timaovae (M) .Tmpdum 
1820  TBpdua*Tnpovac (M+l) 

1825  IF  Tapdua>Tapaax  THEN  Tapdua-Tapmax 

1330  IF  Tnpdun<Tr.pain  THEN  Tapdua-Tapain 

1835  PLOT  Timeevae (K+l) ,TE?dua,-l 
1840  PENUP 
1345  Igo*Igo+l 
1850  GOTC  1495 

1855  S3:  IF  Lkovar (M) >-1000  THEN  M6 
1360  IF  Iplot-0  TEEN  R5 
1865  Atovac-Tiaeovcc  (JI) 

1S70  Alkovar+Lkover (H) 

1375  Iplot-0 
1880  1 

1885  }15:Lkdua-Alkovat 

1890  IF  Lkdua>Lkttax  THEN  Lkdua-Lkmax 

1895  SOME  Tir.:rlr.,Tl-i.-t;v, 

1?C0  HEVC  Atover.Lkdua 
1905  M6:  ! 

1910  IP  Lkovar  (M)— 1000  THEN  M7 
1915  Lkdua-Lkovar <M) 

1920  IP  Lkdua>Lkmax  THEN  Lkdua-Lkaax 
1925  PLOT  Tlmaovar(M) ,Lkdua,-l 
1930  PENUP 

1935  Atovac-Tiaaovar (M) 


1940  AlKovec-Lkover  (M) 

1945  M7s  HOVE  Tiaemax  .L'tmax 

1950  IF  M+1-I»ax  THEN  M9 

1955  IF  M-Imax  THEN  NS 

1960  Iqo-l 

1965  M-M+l 

1970  GOTO  1495 

1975  M9 :  lqo-3 

1980  N-M+l 

1985  GOTO  1495 

1990  MS:  loolot-0 

1995  GOTO  1495 

2000  2 

2005  1 

2010  Testend:  3-FNSoff(4) 

2015  Iflaq-0 
2020  J£lag-0 
2025  Kflag-0 
2030  Lflag-0 
2035  Nn-0 
2040  Delt-0 

2045  PRIHT  "LAST  FILE  FOR  TEST  «";U£ile 
2050  PRIST  "MO. OF  DATA  POINTS  OH  TAPS  «*;Ndata 
2055  PRINT  "Tllcold-";  TlKcld , "Lkold-"  jLltold 
2060  PRIHT  "RELATIVE  TEST  END  TIME  ■" jTima 
2065  CALL  Tina  (Ttir  a , TimeS) 

2070  PRIHT  "TEST  END  TIME-" ,Tia*S 

2075  INPUT  "BEGIN  SEARCH  WITH  FILE*-?" ,N£ile 

2080  ASSIGN  AS(N£ila)  TC  *1 

2085  ON  END  41  GOTO  Filechan 

2090  Idal-0 

2095  IP  Tima<l  THEN  INPUT  "TEST  END  TIME-" ,Ti»a 
2100  IHPOT  "DELETE  ANT  OATA  7", 75 
2105  IF  «-"T"  THEN  Sat 

2110  GOTO  2130 

2115  Sat:  INPUT  "START  TIME  TO  0ELETE  0ATA-",Sd 

2120  INPUT  "END  TIME  TO  OELETE  OATA-" , Ed 
2125  Idal-1 

2130  INPUT  “PERIOD  OF  OATA  AVERAGING-" , Period 
2135  PRINT  "DATA  AVG-" ; Pec iod; "R" 

2140  INPUT  "READ  TESTS ,SATE , Set.  ?",X$ 

21-tS  I?  7S-"7"  THEN  Read 

2150  GOTO  I £5 

2155  Read:  READ  II;  A » ES»CS« C *  E, ? » vj »H » J f K »L 

2160  —  0 TC  If  5 

2165  III :  Iflag-i 

2170  Lflan-i 

2175  READ  #1 ;A 

2180  READ  #1;E 

2165  RE AO  *1;C 

2190  READ  *1:0 

2195  If 5:  Iflag-S 

2200  READ  »1;E 

2205  IF  E— 2000  THEN  Ifl 

2210  READ  SI; F 

2215  READ  S1;G 

2220  READ  *1;H 

2225  IF  Time-S>?eciod  THEN  t£5 
2230  IF  Idel-l  THEN  Delete 
2235  GOTO  2245 

2240  Oeleta:  IF  (S>Sd)  AND  (E<Ed)  THEN  l£S 
2245  IF  Jflaq-0  THEN  EO-E 
2250  3£lag*l 

2255  IF  Lflaq -1  THEN  Dait-0elt+(F-3 ) 

226Q  IF  Lflag-l  THEN  Sn-Nn+1 

2265  IF  Lflaq-1  THEN  Lflad-0 
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2270  fi-fi*f 

2275  G1-G1-K2 

2280  IF  HO-1000  THE:!  2290 
2285  GOTO  2335 
2290  IF  K£lag  *0  THEN  2320 
2295  ai-tH+dC)/2*(£-£l) 

2300  Hsum-Bsum+Bl 
2305  HO-B 
2310  81-B 
2315  GOTO  2335 
2320  El-E 
2325  HO-B 
2330  K£lag-1 
2335  N-N+l 

2340  IF  Tiae-E<.l  THEN  2365 
2345  GOTO  125 

2350  Fllechan:  N£ile-N£il#+1 
2355  ASSIGN  AS(Nfile)  TO  *1 
2360  GOTO  125 
2365  FIXED  5 

2370  PRINT  *T£avg-*;Fl/H 
2375  PRINT  *Tqavg-* ; Gl/N 
2380  Havg-Esmn/(E-SQ) 

2385  PRINT  * Lkavg -* ; Bavg 
2390  PRINT  "Oalt  T«*;Delt/Nn 
2395  FIXED  2 

2400  IF  Idal-l  THEN  PRINT  "DATA  DELETED  FROM* ;Sd; *TO*;£d; *H" 

2405  INPUT  'INITIAL  2ERO  VALUE-*, Iz 

2410  INPUT  ’INITIAL  REF.  TEMP.  VALUE-*, ire 

24 IS  INPUT  'FINAL  ZERO  VALUE-*, Fa 

2420  INPUT  'FINAL  REF.  TEMP.  VALUE*, Frt 

2425  PRINT 

2430  PRINT  "INITIAL  ZERO-*;Iz 

2435  PRINT 

2440  PRINT  'INITIAL  REF.  TEMP-*; It t 

2445  PRINT 

2450  PRINT  'FINAL  ZERO-"jFZ 

2455  PRINT 

2460  PRINT  "FINAL  REF.  TEMP.-*;Ftt 

2465  C t-Gl/N- F  z +1 Z  - . 109* (Fct-Irt) 

2470  PRINT 

2475  PRINT  "CORRECTED  TORQUE-* ;Ct 

2480  STANDARD 

2485  Hgo-2 

2490  GOTO  seart 

2495  ! 

2500  1 

2505  P.azaro:  D-FNRo££(l) 

2510  D-FNRo££ (2) 

2515  D-FSRo££ (4) 

2520  INPUT  *AFX  YOU  PEAOY  TO  RE-ZERO?*, IS 
2525  CALL  Zaco«ud(TNta£) 

2530  PRINT  *Zato  -  * ; Zato 
2535  PRINT  *Thca£  -  *;Thce£ 

2540  Xgo-7 

2545  GOTO  Idle 

2550  1 

2555  i 

2560  ENO 

2565  t 

2570  l 

2575  1 

2560  l 

2585  ! 

2590  I 

2595  SUB  Calibrata 
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2600  ON  ERROR  GOSUB  2610 
2605  GOTO  2620 
2610  PRINT  ERF.M? 

2615  RETURN 

26  20  COM  Zeso, Inliav, "slay , Flag?, Ph2o, Kps, Tqsax,TOEax,TJaax,ThBax,Tvmax 

2625  FOR  J-l  TO  10 

2630  -AIT  50 

2635  BEEP 

2640  NEXT  J 

2645  OIS?  "SET  TORQUE  AT  ZERO’ 

2630  WAIT  5000 

2655  INPUT  ’REA0Y7’,ReadyS 

2660  Avgl-0 

2665  FOR  1-1  TO  50 

2670  Avgl*Avg 1+FNV ( 1 ) 

2675  NEXT  I 
2630  FOR  J-l  TO  10 
2625  (AIT  50 
20SC  SEEP 
2695  NEXT  J 

2700  2 ISP  ’APPLY  50  i:>L2’ 

2705  WAIT  5000 

2710  INPUT  ’REA0Y?' .Ready? 

2715  Avg2«C 

2720  FCF.  1-1  TC  *0 

2725  Avg 2-A vg 2  +FNV { 1 } 

2730  NEXT  I 

2735  Inlbcv*50*5C/(P vg2-Avgl) 

2740  SUBEND 
2745  I 
2730  ! 

2755  DEF  FNTq 

2760  2«  ERROR  GCSUB  2770 

2765  GOTO  2760 

2770  PRINT  ERRH? 

2775  RETURN 

2720  CON  ZaEO,Inlbir.v,Reiay,FlagS,Ph2a,Rptn,Tqnax,ToBax,Tfxax,Thmax,Tvnnax 

2725  V-FNV(i) *Inlb=v-Zeco 

2730  RETURN  V 

2795  FNEHO 

2300  1 

2805  1 

2010  SU3  Zesasub ("beat ) 

2215  CN  ERROR  CC2U3  2825 

22  20  vjCTC 

23  25  PRINT  CRRIiS 

^  -» i 

•  J  vW  UtWlUi 

2325  SIR  T’u(2wC)  ,T1(200) 

2340  COM  Zero,  Inlbav,  Relay  .Flag  S,Ph2o,P.pa,Tgmax,To»ax,TSinax,Ttuiiax,Twinax 
2345  INPUT  ’ARE  YOU  REAOY7’ , Ready? 

2850  OIS?  'CRANK  C.W.  ON  SEEP* 

2855  WAIT  5000 
2360  FOR  J-L  TO  10 
2865  WAIT  50 
2370  BEEP 
2875  NEXT  J 
2880  WAIT  2000 
2385  FOR  M— 1  TO  200 
2890  Th(N)-FNV(lJ 
2395  Thigd-Thighvrh(N) 

2900  NEXT  U 

2905  Tiliqn-Thign/200 

2910  FOR  J-l  TO  10 

2915  SEEP 

2920  WAIT  50 

7.4P4  NRXT  .1 
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2930  OIS?  'STOP* 

293 S  WAIT  5000 

2940  OISP  "CRANK  C.C.W.  ON  BEEP" 

2943  R-FHTref 

2930  Thr  e£-FNTeinp { 4 , R) 

2933  WAIT  5000 
2960  FOR  J-l  TO  10 
296S  WAIT  50 
2970  BEEP 
2973  NEXT  J 
2980  WAIT  2000 
2983  FOR  N»1  TO  200 
2990  TllN)-FNV(l) 

2993  Tlow-Tlow+Tl (N) 

3000  NEXT  N 
3003  Tlow«Tlow/200 
30X0  FOR  J-l  TO  10 
30X3  BEEP 
3020  WAIT  50 
3025  NEXT  J 
3030  015?  "STOP" 

3035  WAIT  5000 
3040  FOR  N-X  TO  200 
3045  S1«S1+<T1(N)-T1ow)“2 
3050  Sh-Sb-MTSiigb-ThlN) )  “2 
3035  NEXT  N 

3060  Zaco- (Thigh+Tlow) /2"lnlbmv 
3063  OISP  "ZERO  CALOJLATED" 

3070  S1-(S1/200)*.3 
3075  Sh-(Sh/200) *.5 

3080  PRINT  “SOL*";  SI  "Inlbav , "SDH-"  ;Sh-Inlbmv,  "Avg. High- ".-Thigh" Inlbnrv  ,"Avg.  Low* 
;Tlo«*tnXbiBv,  “N-“ ;N-X 
3083  SUBEHO 
3090  ! 

3093  ! 

3X00  OSF  FNV(Mch) 

3X05  ON  ERROR  GOSUB  3XX3 
3X10  GOTO  3125 
3XX5  PRINT  ERRMS 
3X20  RETURN 

3X23  1  A-OECIHALUOOIO+Nch) 

3130  OFF  INT  *5 
3X33  waiTS  10  5>  5;32 

3X40  OUTPUT  5  WHS  USING  "» ,W" ; -3984 ,A ,-3920 

3X43  WRITE  10  5,4; 4096 

3X50  REAO  10  5,4; V 

3X35  WRITE  10  5, 5; 128 

3X60  ON  INT  *5,14  CALL  16940 

3X65  OUTPUT  5  WHS  USING  " I ,W" ;-4064 ,32767 

3170  OUTPUT  5  WHS  USING  "* ,W" ;-3936 ,23672 ,-3792 

3X75  CARO  ENABLE  5 

3180  RETURN  V-3C-3 

3X85  FNEND 

3X90  t 

3X95  1 

3200  OEF  FNRoff(N) 

3203  ON  ERROR  GOSUB  3220 

32X0  CCH  Zaco, I nlbmv.Ralay, Flags, Ph2o,Rpn,Tqnax,Toaax,T£fflax,T&max,Twfflax 

3215  GOTO  3230 

3220  PRINT  ERRMS 

3223  RETURN 

3230  A-2“(N-L) 

3235  Jl*l-«-l 

3240  Q-BIT (Relay ,Nml) 

3245  IF  3-0  THEN  GOTO  3300 
3250  l  A— OCTAL  (A)  "OCTAL  (Relay) 


3255  !  Relay-OECIMAL(A) 

3260  OPP  INT  *5 
3265  WRITE  10  3, 5; 32 

3270  OUTPUT  5  WHS  USING  *  * , W* ; -4000 .Relay 

3275  WRITE  10  5,5; 123 

3280  OH  INT  »5,14  CALL  16940 

3235  OUTPUT  5  WHS  USING  *♦ ,W* ; -4064 , 32767 

3290  OUTPUT  5  WHS  USING  ** ,W';-3936, 26572, -3792 

3295  CARO  £«A3LS  5 

3300  RETURN  A 

3305  PNEHO 

3310  l 

3315  t 

3320  QEr  PNRon (N| 

3325  ON  ERROR  GOSUB  3340 

3330  COM  Zeeo.Inlbov, Relay, Plag$,Ph2o,RpB,Tqnax,ToiBax, Tines, Th*ax,Twmax 

3335  GOTO  3350 

3340  PRINT  ERPttS 

3345  RETURN 

3350  A-2*(M-1) 

3355  Nsl-N-1 

3360  B-BIT ( Relay ,Nal) 

3365  IP  8-1  THEN  GOTO  3420 
3370  1  A-OCTAL  ( A)  OCTAL  ( Re  lay ) 

3375  l  Relay “0ECIMAL (A) 

3380  OPP  INT  15 
3385  WRITE  10  5,5; 32 

3390  OUTPUT  5  WHS  USING  * t ,W“ 1-4000 , Relay 

3395  WRITE  10  5,5; 128 

3400  ON  INT  f 5 ,14  CALL  16940 

3405  OUTPUT  5  WHS  USING  "# ,W* ; -4064 , 32767 

3410  OUTPUT  5  WHS  USING  * »,W* ;-3936 ,28672,-3792 

3415  CARO  ENABLE  5 

3420  RETURN  A 

3425  PNENO 

3430  1 

3435  I 

3440  • 

3445  1 

3450  SUB  Tint 

3455  ON  ERROR  GOSUB  3485 

3460  OPTION  BASE  1 

3465  COM  Zero, Inlbov, Relay ,Plag$,Ph2o,Rna,Tqnax,Tomax,T£aax .Thoax ,Tvraax 
3470  COM  T.Tiaeeoi, SHORT  Thout.Toout.Twout .REAL  Lkold.Tlkold , Tile ( 200) , Three , Nik 
3475  COM  160, I60on, 112, Nqo, SHORT  Tiae,TCout,Tqout,Lkout,Tf ,Tq,Th,Z14 ,Vibset,Nvib 
,  Tv  lb ,  Days  ( 12) 

3480  GOTO  3510 

348S  PRINT  ERRMS 

3490  IP  ERRNO160  THEN  3505 

3495  Sr»Ec*l 

3500  IP  Ec "2  THEN  4010 

3505  RETURN 

3510  IP  Plaq$-*Idle*  THEN  GOTO  Exit 
3515  P-PNV<2) '66.66666667 
3520  IP  ?<Ph2o  THEN  Pulse 
3525  GOTO  Check 
3530  t 

3535  Pulse:  D-PNRon(7) 

3540  WAIT  50 
3545  Q-PNRO«f(7) 

3550  Check:  IP  FlaqS-'Stac t*  THEN  Checkh 
3555  Checkl:  IP  P<.7*?h2o  THEN  Down 
3560  Cheekh:  I?  ?>1.2'Ph2o  THEN  Down 
3565  GOTO  0k 
3570  ! 

3575  Down:  PRINT  "SHUT  DOWN  PRESSURE  IS  ";P 


3500  CALI.  Tin*  (The*. Tima S) 

3585  PRIST  *TIME-*;TiaeS,*Lkold-*  ;LkOld,*Tlttold-*;Tlkold 
3590  Nq 0*7 
3595  SOTO  Exit 

3600  Ok:  If  FlagS-’Starc*  THEM  Exit 

3505  T-T+l 

3810  I60on«l 

3615  IP  T-l  THEN  Tae 

3620  IF  T-2  TEES  Rpmo 

3625  I?  T-3  THEt:  Lie 

3630  IF  T-4  THEM  Oilwater 

3635  IP  T-5  TEEM  Ptgtp 

3640  IF  T-6  THEM  160-1 

3645  IP  T-10  THEM  Ptgtp 

3650  IP  T-IS  THEM  Ptgtp 

3655  IF  T-20  THEN  Ptgtp 

3660  IF  T-25  THEN  Ptgtp 

3665  IF  T-30  THEM  Ptgt? 

3670  IP  T-35  the::  Ptgtp 
3675  I?  T-40  THEN  Ptgtp 
3680  IF  T-4 5  THEM  Ptgtp 
3635  IF  T-50  THEN  Ptgtp 
3690  IF  T-35  THEN  Ptgt? 

3655  IF  T-60  THEN  T-0 
3700  IF  T-0  THEN  Ptgt? 

3705  GOTO  Exit 
3710  1 

3715  ptgtp:  I60on-0 

3720  I?  T-5  THEN  Zero 

3725  IP  112-0  THEN  Zero 

3730  GOTO  Aa 

3735  Zero:  Ti-0 

3740  Tti-0 

3745  Tg-0 

3750  112-0 

3755  Aa:  UZ-I12+1 

3760  R-PNTrei 

3765  Tf T-PNTemp (3 . E) 

3770  Tf-re+TCf 
3773  ThL-?NTaap<4,R) 

3780  TM-Th+Thh 

3785  Tqq-?NTq*. 1090* (Thh-TNre£) 

3790  Tq-Tqvrqq 

3795  !  ASSUMES  ZERO  AT  Thref 
3800  IF  Ttf>Tfaax  THEN  Ab 
3803  GOTO  AC 

3810  AH:  PRINT  *SHCT90WN  -  FACE  TEMP  IS  * ,T«f 
3615  CALL  Tiae (Thrs, TineS) 

3820  PRINT  "TINE-*: TimeS,  "Lkold-*  ;Lkold,  "Tlkold-* ;  T1  kald 
3825  Sgo-7 
3830  GOTO  Exit 

3835  Ac:  IF  Thb>Thaax  T3EN  Ad 
3340  GOTO  A« 

3845  Ad:  PRINT  "SHUTDOWN  -  HOUSING  TEMP  IS  *;Thh 
3650  CALI  Tia«(Thcs,Ti:aeS) 

3855  PRINT  *Ti::E-*  ;Tise$ ,  "Lkold-*  ,-Lkold ,  ’Tlkold-* ;  T1  keid 
3C-5 v.  ..go-7 
3:;65  GOTO  Exit 

3870  A«:  IF  Tqq>Tqaax  THEN  A£ 

3875  GOTO  Ag 

3880  At :  PRINT  'SHUTDOWN  -  TORQUE  IS  *?Tqq 
3885  CALL  Tiae (TMcs, Tiae S) 

3890  PRINT  "TIME-" ;TiaeS» "Lkold-* ; Lkold, "Tlkold-" ; Tlkold 
3893  3go-7 
390Q  GOTO  Exit 
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3905  Aq:  IP  T-C  TEES  Out 
3910  IP  I12>12  THEN  Out 
3915  GOTO  Exit 
3920  Out:  Tfout-Tf/II2 
3925  Thout-TSi/1 12 
3930  Tqout-Tq/I12 

3935  Hp-Tqout*2’PI* (1/12) *Rpa/33000 
3940  IF  Hp<5.5  THES  3970 
3945  PRIST  'SHUTDOWN  -  HORSEPOWER  IS’jHp 
3930  CALL  Tine <Th es, TineS) 

3955  PRINT  'TIME-*  ;Tiae$ ,  ’Lkold-’ ;  Lkold,  Tlkold-* r  T1  Isold 

3960  Sqo-7 

3965  GOTO  Exit 

3970  112-0 

3975  GOTO  Exit 

3980  Rpat:  I60on-0  ! AVAILABLE  LATER 

3985  GOTO  Exit  * 

3990  Oilvatec*.  IS0on«0 

3995  R-Ft!Tref 

4000  Toout-FNT«np( S,R) 

4005  Twout-FNTaap(6,R) 

4010  IF  Toout>Tomax  THEN  3a 
4015  GOTO  3b 

4020  3a:  PRINT  ’SHUTDOWN  -  OIL  TEMP  13  *;Toout 
4025  CALL  Tiae  <Thrs,Ti.?.e$ ) 

4030  PRINT  ’TIME-’; TiceS, ’Lkold-* ; Lkold, ’Tlkold-’; Tlkold 
4035  Kqo-7 
4040  GOTO  Exit 

4045  3b:  IF  “wout>Twraax  THEN  Be 
4050  GOTO  Exit 

4055  Be:  PRINT  ’SHUTDOWN  -  V.'ATFP  TEMP  IS  ’;TvOUt 
4060  CALL  T iae (Thrs ,TiraeS ) 

4065  PRINT  ’TIME-* ; Times, "Lkold-’ ; Lkold, ’Tlkold-* ; T1 kold 
4070  Nqo-7 
4075  GOTO  Exit 

4080  Tne:  CALL  Time (Thrs .Tine S) 

4085  Vlbset-0 

4090  Tine-Thrs-Tiaeref 

4095  I60on-0 

4100  GOTO  Exit 

4105  Lk:  ,OFP  1ST  »5 

4110  WRITE  IO  5,5:32 

4115  I60on-0 

4120  IF  Slk-0  THEN  Ea 

4125  Lkout-Nlk’2/(Tlk (Nik) -Tlkold) /60 

4130  Tlkold-Tlk(Nlk) 

4135  Nlk-0 

4140  GOTO  Eb 

4145  Ea:  Lkout— 1000 

4150  Eb:  WRITE  IO  5,5:128 

4155  ON  INT  45,14  CALL  X6940 

4160  OUTPUT  5  WHS  USING  *♦ ,W’ ;-4064, 32767 

4163  OUTPUT  5  WHS  USING  ’♦, W :-3936 , 28672 ,-3792 

♦170  CARO  ENABLE  5 

4175  GOTO  Exit 

4180  Exit:  1 

4185  SUSEXIT 

4190  ! 

4195  1 

4200  OSr  PNTref 

4203  ON  ERROR  GO SUB  4240 

4210  OPTION  BASE  1 

4215  COM  A,a,C.0$,£,F,G.a,l,J,R,L,M 

4220  COM  A1,A2, SHORT  A3,A4,A5 .REAL  A6  ,A7 ,A8 (  200) ,A9 , A10 
4225  COM  31,32,33,34 .SHORT  35 , 36 , 37 , 38 ,39 ,310 , 31 1, RX 
4230  1  SLOT  401,  J  THERMO  0-400C 
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4US  SOTO  4250 
4240  PRINT  ERRMS 
4245  RETURN 
4250  OFF  1ST  1 5 
4255  WRITE  10  5, 5; 32 

4260  OUTPUT  5  WHS  USING  ** ,W* j -3984 ,4103 , -3920 

4265  WRITE  10  5,4; 4096 

4270  READ.  10  5 ,4  ;  R1 

4275  WRITE  10  5,5;128 

4280  ON  I NT  45,14  CALL  16940 

4285  OUTPUT  5  WHS  USING  *» ,W* ; -4064 , 32767 

4290  OUTPUT  5  WHS  USING  •«,W* j-3936 ,28672,-3792 

4295  CARO  ENABLE  5 

4300  IP  (RK3900)  AND  (Rl>3600)  THEN  R3-R1 
430 5  R-RS 

4310  R-(R-4096> *10/4352 

4315  R-3923.7225/(13.1S8-fLOG(-R))-273.2 

4320  R«R* (50.373743+R*( 3. 016701lE-2-R*7 .4293513E-5) ) *lE-6 

4325  RETURN  R 

4330  FHENO 

4335  1 

4340  1 

4345  OEP  FNTetnp(N,R) 

4350  ON  ERROR  GOSUB  4365 
4355  !  SLOT  401  JTHEMO  0-400C 
4360  GOTO  4375 
4365  PRINT  ERRMS 
4370  RETURN 

4375  l  A •DECIMAL (10000+N) 

4380  OFF  INT  »5 
4385  WRITE  10  5,5; 32 

4390  OUTPUT  5  WHS  USING  *»  ,'.J*  ;  -3984  ,A  ,-3920 

4395  WRITE  10  5,4; 4096 

4400  READ  10  5.4;T 

4405  WRITE  10  5, 5; 128 

4410  ON  INT  45.14  CALL  16940 

4415  OUTPUT  5  WHS  USING  * #,W* j -4064 , 32767 

4420  OUTPUT  5  WHS  USING  *  * ,W* ;-3936 , 28672, -3792 

4425  CARO  ENA3LE  5 

4430  T» (T-<T>2047)  *  4096} *  IE-5 

4435  T-(R+T)*IE6 

4440  T*T* ( 1, 9750953E-2+T* (-1. 8542600E-7+T* ( 8 . 368 39 5E-1 2 -T*l. 3280568 E-16)  ) ) 

4445  T-1.3*T+32 

4450  RETURN  T 

4455  PNENO 

4460  1 

4465  1 

4470  l 

4475  l 

4480  SUC  Cooed (Tinenin.Tiasnax ,Trqmin,Trqsiax,L)cnax,Tnpmax,Tspain) 

4485  ON  ERROR  GOSUB  4 SCO 
4490  UNIT  0,270,0,190 
4495  GOTO  4510 
4500  PRINT  EARNS 
4505  RETURN 
4510  l 

4315  FOR  I«1  TO  5 
4520  CALL  Tint 
4525  NEXT  I 

4530  LOCATE  10,130,10,97 

4535  MOVE  70,95 

4340  FOP.  I«1  TO  5 

4545  CALL  Tint 

4530  NEXT  I 

4355  LCRG  3 

4560  CSI2E  3, ,5 


4565  INPUT  *TE3T!r-?*,C 
4570  LABEL  *TEST»*;C 
4575  FOR  1-1  TC  5 
4560  CALX.  Tint 
45<iS  NEXT  I 
4590  IH9UT  *3.V2C?*,.\S 
4««*  faABEL  AS 
-.600  FOR  >1  TC  5 
♦605  CALL  Tint 
4610  ME XT  I 

4615  INPUT  ’MAX  TORQUE-7* ,Ymax 

4620  INPUT  *MIN.  TORQUE-7*  , lain 

4625  Trqmin-Yain 

4630  Trqnax-Ymax 

4635  INPUT  *)1AX  TIME-?*  ,Xaax 

4640  Tiaaaax-Xaax 

4645  INPUT  *KIN  TIME  -7“ .Tiaanin 

4650  Total-Tiaemax-Tinamin 

4655  Q— (1/12)  Total+Tiaemin 

4660  SCALE  Q , Xaax , Ymin , Ymax 

4665  INPUT  "TORQUE  CIV.-?*,Yt 

4670  INPUT  "TIME  DIV?*,Xt 

4675  AXES  Xtt,Yt,Tia«ain,Trqain,  1, 1 

4630  FOR  1-1  TO  5 

4685  CALL  Tint 

4650  NEXT  I 

4695  CSIZE  3, .3 

4700  LORO  6 

4705  :4iddle-(Tia«nax+Tlaeain) /2 

4710  MOVE  Middle, Trqmin-2. 5/50 * (Ymax-Ynir.) 

4715  LABEL  USING  "X" ; "TIME (HRS .) " 

4720  FOR  1-1  TO  5 

4725  CALL  Tint 

4730  NEXT  I 

4735  CSIZE  3,. 5 

4740  LORG  7 

4745  Xnin-Tiaerain 

4750  INPUT  "TIME  STEP-?* ,T 

4755  FOR  Xl-Xain  TO  Xaax  STEP  T 

4760  MOVE  X1+. 25/24‘Tatal+Tiaemin/Tiaenax* (. 05/24*Total) ,Trqmin-l. 5/50* (Tmax-Tmi 
n) 

4765  IE  T<1  THEN  GOTO  4730 
4770  LABEL  USING  “MOOC'jXl 
4775  GOTO  4790 
4780  ! 

4785  LABEL  USING  "MOOU.DD*  ;X-1 

4790  NEXT  XI 

4795  FOR  1-1  TO  5 

4800  CALL  Tint 

4805  NEXT  I 

4810  CSIZE  3,. 5 

4815  LORG  7 

4620  INPUT  'TORQUE  5TEP-?*,C 
4325  FOR  Yl-Yoin  TO  Yaax  STEP  C 
4830  Z-.1/24 

4835  Tocque-Tiaanin-Z'Total-Tiaanin/Tiaeffiax* (. l*Total/24) 

4840  MOVE  Torqua,Yl+.3/50* (Yaax-Yain) 

4845  LABEL  USING  *MCD0"; YX 
4850  NEXT  Y1 
4855  FOR  1-1  TO  ! 

4360  CALL  Tint 
4865  NEXT  I 
4870  CSIZE  3, .5 
4375  LOIR  PI/2 
4660  LORG  S 
4885  Z-l.3/24 
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4390  Torque-Tiaemin-Z^otal-Tiaenin/Tiaenax"  (.  l"Total/24) 

4895  MOVE  Torque , (Yaax+Ymin) /2 

4900  LABEL  USING  * K*  ;* TORQUE  (IN.  -LBS. )  * 

4905  FOR  1-1  TO  5 
4910  CALI.  Tint 
4915  NEXT  I 

4920  INPUT  "MAX  TEMP. -?" ,Yaax 

4925  Tapaax-Yaax 

4930  INPUT  "MIN  TEMP-?" ,Yain 

4935  Tmpmin-Yain 

4940  SCALE  Q,XBax,Yain,Yaax 

4945  INPUT  "TEMP  DIV.«?",Yt 

4950  AXES  Xt.Yt.l'inenin-l/U'Total, 0,1,1 

4SS5  POE  1-1  TO  5 

4960  CALL  Tint 

4965  NEXT  I 

4970  CSIZE  3,. 5 

4975  LOIR  0 

4980  LORG  7 

4985  INPUT  "TEL?  STSP-?",0 
4990  FOR  Yl«Yain  TO  Yaax  STEP  0 

4995  Tem»Tiaeain-l/12*Total-. 1/24 "Total-. l/24"Total*TineBin/Tinienax 
5000  MOVE  Tea, Y1+. 3/50* (Yaax-Yain) 

5005  LABEL  USING  "MDOD";Yl 

5010  NEXT  Y1 

5015  FOR  I»1  TO  5 

SO 20  CALL  Tint 

5025  NEXT  1 

50  30  Yaax-Yaax-t-Yain 

S03S  CSI2E  3, .5 

5040  LOIR  PI/2 

5045  LORO  6 

5050  Ter.-Ticc2in-l/12*Totai-l.  C/24*Total-.  l/24*Total*Tiinemin/Timeniax 
5055  MOVE  Tea, Ymax/2 

5060  LABEL  USING  "8";*SEAL  TEMPERATURE (F ) * 

5065  FOR  r«l  TO  5 
5070  CALL  Tint 
5075  NEXT  I 

5080  INPUT  "MAX  LEAKAGE-?" ,Yaax 
5085  Lkaax-Yaax 
5090  Yain-0 

5095  SCALE  Q ,Xaax , 0 , Yaax 

5100  INPUT  "LEAKAGE  0IV.-?",Yt 

5105  AXES  xaax ,Yt,Xaax,0, 1,1 

5110  FOR  1-1  TO  5 

5115  CALL  Tint 

5120  NEXT  I 

5125  CSIZE  3,  .5 

5130  LOIR  0 

5135  LORG  4 

5140  INPUT  "LEAKAGE  STEP-?" ,E 

5145  FOR  Yl-Ymin  TO  Yaax  STEP  E 

5150  MOVE  Tiaeoax+.3/24*Total, Y1+.3/ 50*Ynax 

5155  LABEL  USING  "MOOO";Yl 

5160  NEXT  Y1 

5165  FOR  1-1  TO  5 

SI 70  CALL  Tint 

5175  NEXT  I 

5180  CSIZE  3,. 5 

5185  LOIR  PI/2 

5190  LORG  6 

5195  MOVE  Tiaaaax+1. 3/24*Total,Yaax/2 
5200  LABEL  USING  "K" ; "LEAKAGE (ML/MIN .) " 

5205  FOR  I -l  TO  5 


5220  MOVE  Tiaecax.Yaax 
32  25  2?.-'.:  Ticsc  in,V.7S:< 

3  230  ?e::o? 

5233  LCC'-TE  10*120/26*2,130,10,97 
5240  SCBEXIT 
5245  1 
5250  ( 

5255  1 

5260  SOB  15940 

5265  Oil  ERROR  GOSU3  5295 

5270  OPTION  BASE  1 

S275  COM  Zeco, lnlbov,Relay,flaq$,Ph2o,Rpn,TqBax,Tomax,Tf:3ax ,Tbraax,Twmax 

5280  COM  T  ,Tinete£  ,3<iOFT  thou t .Taout ,Twout , REAL  LXold , Tlkold ,Tlk ( 200) ,Thte£ ,Xl k 

5235  COM  160, I60cn, 112, Ngo, SHORT  TiiB€,Tfout,Tqout,LKOut,Zll,Z12,Z13,Zl4,Vibset,N 

vib.Tvib 

3290  GOTO  3305 

5295  PRINT  ERRMS 

5300  RETURN 

5305  OUTPUT  5  WHS  USING  *  * ,X* ;-3936 

5310  WRITE  IO  5,4:28672 

5315  READ  IO  5  , 4  ;  X 

5320  X-X+32763 

5325  X«4u95-X 

5330  Y*4 

5335  IP  X-Y> . 5  THEN  5360 
5340  Y-2 

5345  IP  X-YJ-.5  THEN  Cil2 
3350  Y-l 

3353  If  X-Yi-.S  THEN  Chi 

3360  ?CINT  ’EVENT  SENSE  CARO  IS  RETURNING  AN  INCORRECT  SIGNAL, Y-* ; Y 
5365  CTEEXIT 

**  •  i 

5375  ! 

5380  Chi:  Nlk-Nlk+1 
5385  CALL  Tiaa (Thr s ,TioaS) 

5390  Tlk(Nlk) -Thcs-Tiaeref 
5395  WAIT  50 
5400  C-FNP.o£f (4) 

5405  'WAIT  10 
5410  C»FNRon ( 4) 

5413  SL’SEUIT 
3420  1 
5425  ! 

5430  Ch 2: i 

5435  IP  V ibsct»G  THEN  Cycle 

5440  ilv ih»”vih+l 

5445  CALL  Tiae  (Thes,Tiae$) 

3450  Tr.ew«Thrs-TirccsJ 

2455  IP  Tnc’.v-Tvic>  10/3600  THEN  A1 

5<«60  GOTO  Resec 

5463  Alt  I?  ■:vib<S  TIEN  Reset 

5470  PRINT  ’VIErVTi:;:  ALA9N  OUETTCNN* 

3475  PRINT  *Nvlt»  ’:Nvib 

3400  print  *Tnew«  *;Tnew 

5483  PRINT  •T7ib»  *;Tvih 

5  4T0  N;e»7 

5i0  2  CCTC  P.esot 

55C0  Cycle:  Svib-i 

3503  CALL  Tiae  (Thrs,  Ti.r.s?] 

3210  TviL*If.rs-Tltracc£ 

2515  Vibset«l 
5320  GOTO  Reset 
5323  Reset:1. 

5330  SUBEND 
5535  ! 


5540  l 

5545  SUB  Tiase  (Thrs.TiseS) 

5550  OPTION  2AS”  i 

5555  COM  2ero,Inlb«v,Aalay,Flag$,Ph2o,apn!,?qn!ax,Toffax,Tf:nax,ThPax,Twpax 
5560  COM  T,Tir.ecef  ,3UCRT  Thout,Toout,Twout,REAL  Lkold,Tlkold,Tllc(200)  ,Thref  /.Ilk 
5565  COM  160, I6Gon, 112, Hqo, SHORT  T  iaa  ,  Of  out ,  Tqout  ,Lkout,  2 11, 212,  Z13  ,2  14  ,V  idsac  ,N 
vib  ,Tvib ,0ays ( 12) 

5570  l  THIS  IS  TIME  FROM  BEGINING  OF  TEAR 

5575  I  FSB.  IS  TAKEN  AT  28  OATS 

5530  ON  ERROR  GOSUB  3b 

5S3S  Er-0 

5590  Aa:  DISABLE 

5595  OUTPUT  luj*R* 

5600  ENTER  10;>:cnth , Day , Hour  ,Minuta  ,S«cond 

5605  ENABLE 

5610  IF  Er-1  THEN  Aa 

5615  Er-0 

5620  GOTO  Ce 

5625  Bb:  PRINT  ERRMS 

5620  PRINT  Month, Oay , Hour ,Kinuta .Second 
5635  IF  ERRNO160  THEN  Od 
5640  Es-EC+1 
5645  Od:  RETURN 

5650  Cc:  Thrs«(Days  (M-onth)  +Oay-l)  *24+Hour+(Minute+S«cond/60)  /60 
5655  MS-VALS (Month) 

5660  OS-VALS(Oay) 

5665  BS-VALS (Ho  ur ) 

5670  Min$-VAL $  (Ml nuta ) 

5675  TlcaS-MSs*/“iSSi "/80  •iHSi*:" 

5680  IF  Minuta<10  TEEN  E« 

5635  TimaS-TlaaStMinS 
5690  3U6EXIT 

5695  Ee:TlH«S-TiaeS-"0" SHinS 
5700  SUBEMO 


152 


Appendix  B 
Test  Results 
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Appendix  C 
Performance  Curves 
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Figure  C-2,  cont. 
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Figure  C-3.  Test  #40,  1800  RPM,  PH  0  =  3.45  MPa, 
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Figure  C-3,  cont. 


Figure  C-3,  cont. 
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Figure  C-4.  Test  #41,  1800  RPM,  PH  Q  =  3.45  MPa ,  p  =  3.45  MPa- 
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Figure  C-4,  cont. 
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Figure  C-4,  cont. 
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Figure  C-5,  cont. 
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Figure  C-5,  cont. 
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Figure  C-5,  cont 
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Figure  C-6,  cont. 


Figure  C-6,  cont. 
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Figure  C-6,cont. 
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Figure  C-7.  Test  #44,  1800  RPM,  pH  n  =  3.45  MPa,  p  =  8.3  MPa. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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Figure  C-8,  cont. 
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